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SECTION  1 


INTRODUCTION 


1.1  Scope 


The  present  report  gives  an  account  of  the  progress  made  during  the  re¬ 
porting  period  in  four  principal  subject  areas: 

(1)  Large  Space  Structure  (LSS)  system  identification. 

(2)  The  effect  of  actuator  dynamics  on  LSS  control. 

(3)  Modal  dashpot/modal  spring  control  design. 

(4)  A  controlled  experiment  for  ACOSS  design. 

1.2  LSS  Identification 

In  Section  2,  the  basic  approaches  to  LSS  identification  and  the  major 
difficulties  encountered  in  the  identification  are  discussed.  Subsequently, 
attention  is  focused  on  the  high-resolution  identification  problem  typifying 
LSS.  A  novel  procedure  for  identifying  lightly  damped,  closely  spaced  modes 
is  derived.  The  approach  is  rooted  in  the  notion  of  numerical  rank  of  matrices 
which  has  bearing  upon  the  spectral  estimation  problem.  Optimizing  the  sampl¬ 
ing  rate  and  relaxing  the  Nyquist  rate  to  construct  a  "zoom-in"  procedure  are 
utilized  in  deriving  efficient  and  numerically  stable  high-resolution  spectral 
analysis  algorithms. 

Although  the  procedures  reported  here  are  obtained  in  the  context  of  LSS 
identification,  they  can  also  be  applied  to  arbitrary  high-resolution  spectral 
analysis  problems  such  as  radar  signal  processing. 

1.3  Actuator  Dynamics 

The  common  assumption  in  work  on  LSS  control  is  the  availability  of  ideal 
sensors  and  actuators.  Namely,  the  bandwidth  and  the  power  of  the  devices  are 
assumed  unlimited  and  structure/actuator  interaction  is  ignored.  Section  3  re¬ 
ports  our  observations  .<a  the  effects  of  incorporating  the  actuator  dynamics  in 
the  model.  The  proper  way  to  model  conventional  actuators  is  presented  first. 
Subsequently,  two  approaches  to  the  incorporation  of  the  actuator  dynamics  are 
explored. 

(1)  The  structural  model  is  augmented  by  the  actuator  model,  and  the 
controller  is  designed  for  the  augmented  plant.  It  is  seen  that 
the  augmented  plant  has  the  same  dynamic  representation  as  struc¬ 
tures  with  optimal  actuators.  Thus,  the  conventional  LSS  control 
design  approaches  are  applicable  to  the  augmented  system. 
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(2)  A  controller  is  designed  assuming  ideal  actuators,  and  a  compensator 
is  designed  around  the  actuator.  The  design  Issues  for  actuator 
compensators  are  addressed  and  some  solutions  are  proposed. 

1.4  Modal  Spring  and  Modal  Dashpot  Design 

Section  4  is  complementary  to  the  second  ACOSS  11  semiannual  report 
[1-1],  where  various  designs  to  accommodate  a  broadband  disturbance  are  con¬ 
sidered.  In  the  present  section,  a  procedure  for  a  combined  velocity  and  posi¬ 
tion  feedback  (dashpot  and  spring,  respectively)  are  to  be  designed  to  augment 
damping  and  move  frequencies  of  the  closed-loop  plant.  A  root  locus  procedure 
for  assessing  the  impact  on  the  full  order  closed-loop  eigenvalue  behaviour  is 
also  described. 

1.5  A  Controlled  Experiment  for  ACOSS  Design 

Our  previous  efforts  to  design  a  control  system  to  accommodate  a  broad¬ 
band  disturbance  and  achieve  stringent  line-of-sight  (LOS)  error  requirements 
( [ 1—1 J  and  Section  3  of  this  report)  have  achieved  only  partial  success. 


Either  the  control  design  has  met  the  LOS  error  requirement  for  the  design 
model,  but  proved  unstable  when  applied  to  the  full-order  LSS  (adverse  spill¬ 
over  effect);  or  spillover  effect  has  been  accommodated,  only  to  fail  to  meet 
LOS  error  requirements.  We  attribute  lack  of  success  to  certain  mutual  incom¬ 
patibilities  in  the  distinct  processes  of  reduced-order  model  selection, 
sensor/actuator  placement,  and  reduced-order  controller  design. 

In  Section  5,  an  account  is  given  of  a  procedure  we  have  devised  to  com¬ 
bine  the  various  selection  and  design  steps  to  simultaneously  achieve  spillover 
accommodation  and  meet  LOS  error  requirements.  We  have  not  reported  here  on 
the  progress  made  on  the  design  itself.  This  will  be  discussed  when  the  itera¬ 
tive  procedure  has  converged  into  an  acceptable  design. 

1.6  Perspective 

An  attempt  is  made  to  document  major  achievements  and  observations,  even 
if  they  are  merely  stepping  stones  towards  the  design  and  implementation  of 
identification  and  active  control  of  LSS.  Such  stepping  stones  might  be  dis¬ 
carded  in  the  study  if  they  are  found  to  have  little  merit. 

Some  work  in  progress  has  not  been  discussed  in  this  report  since  it  has 
not  matured  enough  for  documentation.  In  the  context  of  identification  for  ex¬ 
ample,  the  procedures  for  recovering  the  mode  amplitudes  at  the  sensors,  how 
amplitude  information  would  be  used  to  construct  mode  shapes,  and  our  work  in 
input-output  characterization  from  measurement  and  the  subject  of  input  design 
have  not  been  reported. 

The  discussion  on  actuator-structure  dynamic  interaction  is  the  first,  to 
our  knowledge,  in  the  ACOSS  literature.  It  is  the  first  step  towards  the  rig¬ 
orous  analytic  incorporation  of  actuators  into  models.  More  work  is  being  done 
to  evaluate  its  impact  on  ACOSS  design. 


SECTION  2 


STRUCTURE  IDENTIFICATION— THE  RESOLUTION  PROBLEM 
AND  ITS  SOLUTION 


2.1  Introduction 


The  organization  of  Section  2  is  as  follows.  Section  2.1  discusses  the 
background  and  difficulties  encountered  in  structure  identification  problems 
and  their  solutions.  Section  2.2  gives  the  background  of  the  role  of  the  sam 
pling  rate  in  system  identification  and  spectral  estimation.  In  Section  2.3, 
we  discuss  the  effect  of  sampling  rate  on  the  fast  Fourier  transform  (FFT). 
This  is  a  rather  straightforward  discussion.  Section  2.4  establishes  the 
foundation  for  parametric  spectral  analysis  for  a  sampled  data  stochastic  pro 
cess.  The  role  of  the  sampling  rate  within  the  classical  Nyquist  constraints 
is  discussed  in  Section  2.5.  A  zoom-in  procedure  is  described  in  ction 
2.6.  Generic  examples  are  brought  in  Section  2.7.  The  foregoing  :ctions 
were  devoted  to  narrowband  signals.  Section  2.8  describes  a  deci  ion  tech¬ 
nique  to  enable  the  analysis  of  broadband  signals. 


In  applying  identification  to  Draper  model  simulated  data,  ve  so 

far  demonstrated  limited  success  in  the  sense  that  we  have  succes  r  ,.y  iden¬ 
tified  reduced-order  models,  where  the  reduction  in  order  corresponds  to  a 
failure  to  resolve  closely  spaced  modes.  The  decimation  procedure  and  the 
zoom-in  algorithm  are  presently  being  tested  against  Draper  model  structural 
data.  Thus,  the  report  of  the  experiment  results  is  deferred  for  future 
documentation. 


2.1.1  Background 

One  of  the  first  tasks  to  be  undertaken  in  the  quest  for  large  space 
structure  (LSS)  deployment  and  application  is  the  establishment  of  procedures 
for  modelling  such  structures.  Analytical  modelling,  although  an  essential 
step,  does  not  provide  sufficient  Information  for  practical  utilization  of 
LSS  for  a  few  reasons. 

(1)  Some  phenomena,  like  the  damping  mechanism,  high-frequency  behav¬ 
ior,  and  joint  dynamics,  are  not  well  understood,  and  thus  cannot 
be  reliably  incorporated  into  analytical  models. 

(2)  LSSs  are  properly  modelled  by  distributed  parameter  systems  (par¬ 
tial  differential  equations).  However,  such  models  are  often  too 
difficult  to  construct,  and  are  usually  not  amenable  to  control 
design.  Furthermore,  the  state-of-the-art  of  identifying  distri¬ 
buted  parameter  systems  has  not  matured  enough  for  routine  appli¬ 
cations.  Hence,  the  usual  models  for  LSS  are  of  the  lumped- 
parameter,  finite-dimensional  type.  Such  models  can  be  derived 
using  an  approximation  of  the  distributed  parameter  model,  e.g., 
using  finite-element  techniques,  or  can  be  obtained  directly  from 
the  physical  structure,  say,  by  representing  it  as  a  mesh  of  rigid 


bodies  subject  to  massless  elastic  constraints.  The  final  model  is 
often  represented  via  the  vibration  equation 

Mx  +  Dx  +  Kx  =  f  (2--1) 

The  evaluation  of  such  an  approximation  and  the  estimation  of  the  unknown 
entities  in  Eq.  (2-1)  (e.g.,  the  D  matrix)  will  require  the  implementation  of 
a  parameter  estimation  technique,  henceforth  referred  to  as  identification. 

The  problem  of  identification  can  be  philosophically  addressed  from  two 
points  of  view. 

(1)  Model  Validation — Construct  an  experiment  and  employ  a  procedure  to 
estimate  the  parameters  of  a  model  constructed  from  physical  con¬ 
siderations.  For  example,  estimate  the  mass  and  stiffness  matrices 
in  Eq.  (2-1).  On  a  more  fundamental  level,  estimate  the  elastic 
properties  of  the  structure,  such  as  Young's  moduli. 

(2)  Engineering  Modelling — Estimate  input/output  characteristics  of  the 
"plant",  i.e.,  the  LSS  response  to  actuator  or  disturbance  forces 
as  observed  via  a  given  set  of  sensors.  Both  actuators,  sensors, 
and  their  dynamic  characteristics  and  placement  on  the  structure, 
might  be  determined  based  on  control  considerations  using  an  ana¬ 
lytical  approximation  of  the  plant. 

Obviously,  the  two  approaches  are  not  disjoint.  Information  extracted 
from  either  approach  might  be  utilized  beneficially  to  construct  the  comple¬ 
mentary  model. 

This  section  concentrates  on  the  extraction  of  f requency/damping  infor¬ 
mation  from  LSS  measured  data.  Such  information  ij  the  basis  for  constructing 
engineering  models  as  discussed  in  Option  2  above. 

2.1.2  Difficulties 

There  are  a  few  features  of  LSS  which  hamper  the  extraction  of  fre¬ 
quency/damping  information  from  experimental  data. 

(1)  The  "Curse  of  Dimensionality" — As  already  discussed,  LSSs  are 
essentially  infinite  dimensional.  Even  if  attention  is  restricted 
to  the  bandwidth  of  the  sensors  and  actuators,  the  number  of  modes 
to  be  identified  is  still  large.  For  example,  Figure  2-1  gives  the 
frequency  response  of  Draper  Model  if 2  to  a  15-Hz  disturbance.  (See 
Reference  2-28  for  more  information.)  The  identification  of  a 
high-dimensional  system  has  so  far  not  been  successfully  demon¬ 
strated  in  the  literature. 

(2)  Closely  Spaced  Modes — The  signal-processing  literature  has  devoted 
a  considerable  amount  of  attention  to  constructing  algorithms  which 
are  capable  of  estimating  sinusoids  with  small  frequency  separa¬ 
tions  (closely  spaced  modes).  As  is  obvious  from  Figure  2-1,  this 
is  a  fundamental  problem  for  LSS  identification. 


(3)  Small  Damping — The  ’  dentif ication  of  poles  close  the  the  ju  axis 
has  also  been  a  challenge  to  the  signal  processor.  The  difficulty 
in  estimating  small  damping  is  nicely  explained  in  Reference  2-29 
using  statistical  arguments.  The  damping  ratios  commonly  assumed 
in  the  context  of  LSS  are  of  the  order  of  0.1  to  1  percent. 

(4)  High  Process  and  Measurement  Noise — The  present  state-of-the-art  in 
sensor  technology  in  the  very  low  frequency  range  for  LSS  applica¬ 
tion  suggests  that  the  quality  of  measurements  to  be  used  in  the 
identification  process  might  require  an  estimation  technique  tai¬ 
lored  for  a  high-noise  environment. 

One  or  more  of  the  above-mentioned  problems  appear  in  a  variety  of  iden¬ 
tification  applications.  It  is  the  culmination  of  all  four  difficulties  which 
makes  the  identification  of  LSS  especially  challenging. 

Difficulties  (2)  and  (3)  can  be  classified  as  spectral  resolution  prob¬ 
lems.  In  this  section,  we  report  a  novel  approach  to  the  high  resolution 
spectral  analysis  problem  developed  at  CSDL.  The  algorithm  developed  is  espe¬ 
cially  suited  for  handling  difficulties  (2)  and  (3).  However,  it  is  also 
shown  to  be  noise  resilient.  The  modification  of  the  algorithm  to  handle  high 
dimensional  models  is  also  discussed. 

Although  we  have  developed  our  high  resolution  spectral  analysis  tech¬ 
nique  in  the  context  of  LSS  identification,  it  should  be  noted  that  the  pro¬ 
cedure  is  valuable  in  any  signal-processing  application,  where  high  resolution 
is  desirable,  as  in  radar  signal  processing  or  geophysical  data  analysis. 

2.2  Spectral  Estimation  and  Sampling 


High  resolution  spectral  estimation  has  been  the  subject  matter  of 
numerous  recent  publications  in  the  signal-processing  literature.  The  wide 
range  of  applications  and  the  theoretical  interest  in  the  subject  are  reflect¬ 
ed  in  the  special  issue  on  spectral  estimation  of  the  IEEE  Proceedings  [2-1]. 

In  surveying  the  literature  devoted  to  spectral  estimation,  we  have 
found  that  the  premise  of  most  of  the  work  is  the  availability  of  a  sequence 
of  uniformly  sampled  data  points  of  a  continuous  (stochastic)  process.  The 
analysis  and  the  subsequently  constructed  algorithms  treat  the  discrete  time 
process  while  ignoring  the  sampling  rate  issue. 

It  is  the  purpose  of  this  section  to  demonstrate  the  crucial  dependence 
of  spectral  resolution  on  the  sampling  rate.  It  is  shown  that  by  examining 
the  sampling  rate  issue  we  obtain,  using  a  computationally  simple  algorithm, 
spectral  resolution  as  good  as  with  other  more  complex  procedures. 

To  be  more  specific  in  the  discussion  following,  let 
[y(t) :  0  <  t  <  tf]  ; 

1  I 


t^  finite  or  ® 


be  the  (stochastic)  process  to  be  analyzed.  Digital  spectral  analysis  is  per 
formed  on  the  data 

[y(kT):  k  =  0,1,..  .,N]  ;  N  finite  or  » 

Intuitively,  the  information  contained  in  the  data  is  inversely  proportional 
to  T.  Furthermore,  if  the  process  is  band-limited  (tf  -  »  in  this  case) 
by  the  Shanon  interpolation  formula,  the  data  to  represent  the  process  if 
T  <  T  *  ir/BW,  where  Nyq  abbreviates  Nyquist,  and  BW  =  bandwidth  (in  rad/s) 

Since  tf  is  usually  finite,  aliasing  can  be  reduced  only  by  choosing  T  to  be 
as  small  as  possible.  We  believe  that  these  mathematically  sound  and  intui¬ 
tively  appealing  statements  have  hindered  the  use  of  T  as  a  parameter  in 
improving  spectral  resolution.  Specifically,  most  of  the  work  on  spectral 
estimation  is  based  on  the  availability  of  the  sequence  [y(k)j,  i.e.,  T — the 
sampling  interval  is  usually  suppressed  and  subsequently  ignored  (e.g..  Refer 
ence  2-1). 

The  sampling  rate  plays  a  role  in  nonparametric  spectral  estimation, 
though  not  as  dramatic  a  role  as  in  the  parametric  analysis.  In  Section  2.3, 
the  role  of  T  when  spectral  analysis  is  performed  using  the  conventional  dis¬ 
crete  Fourier  transform  (DFT)  will  be  discussed.  The  discussion  is  rather 
straightforward  and  probably  would  not  surprise  practitioners  in  the  field. 

We  bring  up  this  result  only  for  the  sake  of  completeness. 

The  dramatic  role  of  the  choice  of  T  is  in  the  domain  of  parametric 
spectral  estimation.  Assuming  tf  large  enough  to  enable  resolution,  say,  of 
two  sinusoids  in  white  noise,  then  the  problem  of  resolution  is  fundamentally 
a  numerical  one.  Namely,  regardless  of  the  sophistication  and  statistical 
soundness  of  an  algorithm  employed,  if  the  frequency  separation  is  small 
enough,  the  algorithm  will  fail.  To  demonstrate  this  point  consider  the  fol¬ 
lowing: 

Example  1 :  We  have  applied  covariance  least  squares*  to  128  samples 
of 

y(t)  =  sin  (ujjt)  +  sin  (o^t)  +  v(t) 

with  1-percent  frequency  separation,  i.e.,  e  •  (u i  -  w2)/u)f  =  0.01.  For 
v(t)  =  0,  the  frequencies  wf  and  W2  were  successfully  Identified.  How¬ 
ever,  with  v(t)  white  Gaussian  noise  with  Var[v(t)J  =  10“^,  the  algorithm 
failed. 

The  reason  for  this  failure  is  explained  in  Section  2.5.  In  Section 
2.4,  the  problem  of  parametric  spectral  estimation  for  sampled  data  is  formu¬ 
lated  and  some  modelling  questions  are  discussed.  Section  2.5  is  devoted  to 

*This  algorithm  has  been  proposed  by  Beex  and  Scharf  [2-2)  and  by  Cadzow 
[2-3].  It  will  be  further  discussed  in  the  following. 


the  question  of  spectral  resolution  as  a  function  of  the  sampling  rate  for 
T  <  Tfjyq.  In  Section  2.6,  a  zoom-in  procedure  is  described.  The  procedure 
allows  for  a  choice,  T  >  Tfjyq,  for  the  purpose  of  resolving  closely  spaced 
frequencies.  Our  observations  and  algorithms  are  demonstrated  via  simulations 
in  Sections  2.7  and  2.8.  Section  2.9  contains  concluding  remarks. 

2.2.1  Historical  Perspective 

Although  we  have  failed  to  find  articles  addressing  the  impact  of  sam¬ 
pling  rate  on  spectral  resolution,  the  general  problem  of  sampling  rate  in 
parameter  estimation  (system  identification)  has  received  attention  in  the 
literature. 

We  have  found  the  most  extensive  treatment  of  the  problem  in  the  econo¬ 
metric  literature.  A  compendium  of  results  on  the  subject  can  be  found  in 
Bergstrom's  book  [2—4 1 .  The  main  concern  in  the  econometric  work  is  the  ap¬ 
proximation,  or  the  exact  representation  of  continuous  stochastic  differential 
equations  by  difference  equations.  Consequently,  an  effort  is  made  to  assess 
the  quality  of  the  identification  of  the  continuous  equation  as  inferred  from 
the  estimation  of  the  discrete  approximation.  This  evaluation  is  performed 
via  bias  and  variance  evaluation  for  T  +  0.  However,  this  work  does  not  pro¬ 
vide  us  with  useful  tools  in  the  context  of  high  resolution  spectral  analysis. 

Some  work  on  the  sampling  rate  issue  has  been  done  and  some  practical 
experience  gained  by  researchers  in  the  system  identification  field.  Notable 
is  the  work  by  Astrom  [2-5]  Goodwin,  Pain,  and  Zarrop  [2-7].  Ast’rom  has  anal¬ 
yzed  a  first-order  stochastic  process,  the  parameters  of  which  are  identified 
using  maximum  likelihood  (ML)  estimation  employing  the  sampled  data.  He  has 
shown  that  the  error  variance  is  minimized  for  the  choice  T  *  a  =  0.797  where 
a  is  the  time  constant  of  the  process.  Goodwin,  et  al.,  have  addressed  the 
problem  of  the  experiment  design,  including  sampling  rate  determination,  data 
filtering,  and  optimal  input  design  using  measures  of  the  information  matrix. 
They  have  shown  that  with  an  anti-aliasing  filter  and  a  finite  number  of  data 
points,  the  Fisher  information  matrix  is  inversely  proportional  to  T,  and 
thus,  is  maximized  by  the  choice  T  -  (1  -  e)T}jyq,  e-  being  an  arbitrarily 
small  value.  This  result  coincides  with  our  observations  in  Section  2.3. 
However,  it  does  noc  agree  with  our  observations  on  parametric  spectral  esti¬ 
mation  (Section  2.5  and  2.6)  nor  with  practical  considerations  for  spectral 
resolution. 

We  should  note  that  practitioners  of  spectral  estimation  have  often  used 
good  "rules  of  thumb"  for  obtaining  high  resolution,  although  we  could  not 
find  the  justification  for  this  practice.  For  example,  Bendat  and  Pearsol 
suggest  ([2-8],  page  288)  T  «  0.5  Tfjyq  if  "the  correlation  function  has 
frequencies  near  l/(2Tfjyq),"  and  T  *  0.8  Tfjyq  "if  power  spectra  measure¬ 
ments  are  of  prime  consideration"  (compare  these  to  Sections  2.5  and  2.3, 
respectively).  In  the  context  of  identification,  Iserman  [2-9]  recommends 
T  ■  (5...15)*Tg5  where  T95  is  the  95-percent  settling  time  of  the  tran¬ 
sient  function.  This  we  found  to  be  impractical  for  very  low  damping  ratios. 


Finally,  an  important  result,  noteworthy  and  relevant  to  our  problem, 
was  reported  recently  by  Kay  [2-10]:  In  estimating  the  autocorrelation 
function,  the  sampling  rate  should  be  at  least  twice  the  Nyquist  rate 
(T  <  0. 5  Ttfyq),  but  statistically,  there  is  insignificant  gain  in  sam¬ 
pling  faster  than  that  rate.  This  result  coincides  with  the  Bendat-Pearsol 
recommendation  cited  previously  and  nicely  complements  our  observation  in 
Section  2.5. 

2.3  Nonparametric  Spectral  Estimation 


The  most  widely  used  tool  in  spectral  analysis  is  the  discrete  Fourier 
transform  implemented  via  the  FFT  algorithm  on  N  samples  of  a  continuous  sig¬ 
nal.  Typical  commercial  spectrum  analyzers  [2-11  to  2-13]  are  based  on 
N  *  2*®  =  1024  point  FFT.  In  analyzing  the  effect  of  the  sampling  interval 
T  on  spectral  resolution  when  using  the  FFT,  the  following  facts  should  be 
noted. 

(1)  The  Fourier  transform  is  represented  by  N/2  +  1  points  spaced  at 
frequency  intervals  Aw  *  2tt/NT.  Thus,  the  relative  resolution 
around  a  frequency  wq  is 


where  n  is  the  number  of  cycles  of  a  sinusoid  with  frequency  wq. 

(2)  The  width  of  the  main  lobe  associated  with  N  points  sampled  with 
interval  T  is 

<5o )  *  4n  /NT 

Consider  the  signal  of  Example  (1).  The  two  frequencies  are  resolvable 
in  light  of  the  above  facts  if 


T  >  T, 


N  •  (u>1  -  u>2) 


Thus  for  N  *  1024  and  (wj  -  W2)/wi  *  0.01,  frequency  resolution  is 
attainable  for  T  >  0.391  T{jyq.  Note,  however,  that  aliasing  the  main  lobe 

N  -  2 

with  itself  is  prevented  if  T  <  T  •  — - —  which  is  0.2  percent  faster  than 
the  Nyquist  rate.  ^ 


In  the  above  discussion,  we  have  completely  ignored  the  effects  of  leak 
age  which  will  require  even  tighter  bounds  than  presented.  From  the  above 
discussion,  0.688Tjjyq  is  the  sampling  rate  required  to  resolve  W2  from 


Wj  if  a m  9/11  and  N  *  32. 
0.975  are  given  in  Figure  2-2. 


The  FFTs  for  T/T. 


0.3,  0.35,  0.45,  and 


The  preceding  discussion  is  admittedly  rather  elementary.  We  have  pre¬ 
sented  it  to  illustrate  the  often  overlooked  upper  bound  on  the  sampling 
rate.  Furthermore,  it  demonstrates  in  the  frequency  domain  our  main  results 
to  be  derived  below  in  the  context  of  parametric  spectral  estimation. 

It  should  be  noted  that  the  preceding  results  stem  from  the  practical 
constraint  on  finite  N.  If  N  is  not  limited,  T  can  be  chosen  arbitrarily 
small.  However,  even  for  N  -*•  ®  in  the  context  of  parametric  spectral  estima¬ 
tion,  T  should  be  carefully  chosen. 

2.4  Sampled  Data  Spectral  Analysis 


In  this  section,  we  shall  establish  the  relationship  between  a  contin¬ 
uous  time  stochastic  process  and  the  discrete  time  representation  of  its  sam¬ 
ples,  which  is  appropriate  for  spectral  estimation.  Subsequently,  the  param¬ 
etric  spectral  estimation  problem  is  stated.  This  section  is  the  basis  for 
assessing  the  impact  of  the  sampling  rate  on  spectral  analysis  as  presented  in 
Sections  2.5  and  2.6. 

2.4.1  Continuous  and  Discrete  Models 


System  identification  in  general,  and  spectral  estimation  in  particular, 
are  usually  based  on  the  second-order  statistics  of  the  data.  Thus,  in  deter¬ 
mining  a  suitable  discrete  time  representation  of  a  continuous  time  process,  a 
modelling  criterion  that  Perl  and  Scharf  call  "covariance  invariance"  [2-14] 
plays  a  fundamental  role.  Covariance  invariant  models  are  discussed  first. 

It  is  assumed  throughout  this  section  that  the  discussed  processes  are 
wide  sense  stationary  (WSS). 

We  shall  use  the  subscripts  "c"  and  "d”  to  denote  continuous  and  dis¬ 
crete  time  entities  respectively.  Continuous  time  is  denoted  by  t,  xeR+ 

(the  positive  reals),  and  the  discrete  time  by  k,  kel+  (the  positive  inte¬ 
gers).  k  actually  represents  the  time  instant  kT  where  T  is  the  sampling 
interval  being  suppressed  if  it  is  not  discussed  explicitly. 

A  stochastic  (continuous  or  discrete  time)  process  with  a  rational  spec¬ 
trum  can  be  formally  represented  in  the  state  space  via  the  state  equation 

D*  x(t)  *  A  x(t)  +  w(t)  ;  xeRn,  u>eRn  (2-2) 


the  output  equation 

y(t)  -  C  x(t)  ;  yeR*  (2-3) 

and  the  measurement  equation 

z(t)  -  y(t)  +  v(t)  (2-4) 

i 

Or,  in  terms  of  the  impulse  response  convolved  with  the  input  o>(t) 


y(t) 


h(t)  *  w(t)  ;  z(t) 


y(t)  +  v ( t ) 


(2-5) 


which  corresponds  to  the  transfer  function  representation 


Y(q) 


H(q)  W(q)  ;  Z(q)  -  Y(q)  +  V(q) 

-1 


(2-6] 


H(q)  -  C(ql  -  A) 

The  various  entries  in  Eq.  (2-2)  to  (2-6)  are  interpreted  as  follows: 


Term 


Continuous  Time 


Discrete  Time 


time  t 

D 


teR 


differentiation  -r- 
dt 


U),  v 


Wiener  processes* 


kel 

forward  shift  Z 
uncorrelated  white  noise 


t  >  0 


h  (t)  «  Ce 
c 


At 


I  ct(i) 


X(i)r 


h  (k)  -  CAk  -  l  a(i)X(i)k 


h(t) 


t  <  0 


h  (t)  -  0 
c 


hd(k)  -  0 


s  »  the  Laplace  transfer 
variable 


the  z  transform  variable 


The  input  noise  w  and  the  measurement  noise  u  are  assumed  to  be  white 
and  mutually  uncorrelated  with  covariance  matrices 


E[w(t1)w(t2) ' ]  *  W6(t i  -  t2) 

E[u(t1) u(t2) ' ]  -  VCit!  -  t2) 


where  6(t)  is  the  Dirac  delta  function  in  the  continuous  time  case,  and  a 
Kronecke  delta  in  the  discrete  time  case. 


Of  special  interest  is  the  output  covariance  function  associated  with 


Eq.  (2-2)  and  (2-3),  i.e.,  Ryy( t)  =  E[y(t1)y(t1  +  t ) ’ J 


Ce  AtV  C'  t  <  0 

Rc(t)  =  a't 

CV  e  C’  t  >  0 
c 


Rd(k) 


CA'^V.C'  k  <  0 
a 


CVd(A')kC'  k  >  0 


(2-7) 


(Note  that  in  (2-7)  we  have  replaced  the  subscript  yy  by  d  or  c  to  indicate 
discrete  or  continuous  time.) 


*  We  are  using  here  formally  and  somewhat  loosely  constructed,  the  representa¬ 
tion  of  the  continuous  stochastic  process.  More  precise  treatment  can  be 
found  in  Reference  2-7. 


V  «  E[x(t)x(t)'j  is  the  state  covariance  matrix,  satisfying  the  Liapunov 
equation 


A  V  +  V  A'  +  R  -0 

C  C  C  C  0)10 


V.  -  A.V.A'  +  R  (2-8) 

d  d  d  d  oxo 


or  equivalently 


«°At  A' t  <x> 

V  *  /  e  C  R  e  C  dx  V .  -  l  A7R  (A' )k 

c  J  wo)  d  L  d  (0(0 v  d ; 

o  o 


(2-9) 


In  terms  of  the  impulse  response  and  the  transfer  function  we  have 


R  (t  )  -  /  h  ( t )  R  h  (t  +  t)'  dt  R.(k)  -  J]  h.(£)R  h  (£  +  k) ' 

c  coxoc  d  d  <oo>  d 


£=■0 


(2-10) 


n  A  (i)x 

l  0c(i>e  C 

i-1 


l  8d(i)Ad(i) 


i-1 


(2-11) 


~  /H  (-s)R  H  (s)'e8T  ds 

Zn  j  '  c  UMO  C 


Sj/W*"1)11 . H(*),*k'1  d: 


0)0) 


(2-12) 


where  the  coefficients  8(i)  satisfy 


Mi>  *  'I  ac(i>Ru«oac(j),/l^(i)  +  Uj>) 


J-l 


(2-13) 


6d(i) 


l  ad(i)R0)0)°d(j),/l1  '  XH<i)^<J>l 


j“l 


d'~#  <f 


Definition :  (Rephrasing  [2-14]).  A  WSS  discrete  time  stochastic  pro¬ 
cess  is  T-covariance  invariant  (TCI)  with  its  continuous  counterpart  if 


R  (kT)  =*  R.(k)  -V-  k 

c  d 


(2-14) 


In  this  section  we  shall  not  delve  into  the  theory  of  covariance  invar¬ 
iance,  but  rather  state  a  few  facts  (without  proof)  which  can  be  derived  from 
Eq.  (2-3)  to  (2-14). 


Let  E.  -  (A.,  C, ,  B. )  and  E  *  (A  ,  C  ,  B  )  be  the  triplets  characterizing  a 
cl  a  a  a  c  ccc 

discrete  time  and  a  continuous  time  stochastic  process  as  defined  above. 

Here,  Bj  and  Bc  represent  the  square  root  of  -  the  input  covariance 

matrix  for  the  discrete  time  and  continuous  time  cases,  respectively.  We 

ignore  Rvv  since  its  inclusion  is  straightforward. 

Fact  A  Ed  is  T  covariance  invariant  with  Ec  if  the  following  are 

satisfied 


A  T 
•  c 

A^  ■  e  (necessary  condition) 

(2-15) 

C,  *  C 
d  c 

(2-16) 

C.V.  =  C  V 
d  d  c  c 

(2-17) 

where  *  means  equality  up  to  a  similarity  transformation.  Equations  (2-16) 
and  (2-17)  give  sufficient  conditions.  It  is  not  clear  that  Eq.  (2-16)  and 
(2-17)  are  necessary. 

Fact  B  Let  Eq.  (2-15)  through  (2-17)  be  satisfied.  A  necessary  condi¬ 
tion  for  the  existence  of  Ed,  which  is  covariance  invariant  with  Ec>  is 

p(B.)  -  p(C  )  (2-18) 

d  c 


where  p(A)  denotes  the  rank  of  A. 

Fact  B  follows  from  stochastic  realization  theory  and  from  the  spectral 
factorization  theorem  [2-7].  This  statement  imposes  a  condition  on  the  mini- 
I  mum  number  of  uncorrelated  inputs  which  are  required  to  realize  a  discrete 

time  covariance  invariant  representation  of  a  continuous  WSS  stochastic  pro¬ 
cess.  Note,  in  particular,  that  replacing  Eq.  (2-17)  by  Vc  =  Vd  [2-14] 
implies  p(Bd)  *  n.  Fact  B  has  ramifications  beyond  the  scope  of  this  sec¬ 
tion  and  will  be  discussed  in  a  separate  publication. 

Fact  C  Let  p(Cc''  ■*  p(Bd)  ■  1  (Single  Output  and  thus  Single  Input 
by  Fact  B).  The  relative  degree  of  Hd(z-*),  which  is  T  covariance  invari¬ 
ant  with  H^(s ) ,  is  almost  always  one. 

Fact  C  is  a  crucial  observation  which  is  too  often  overlooked  in  the 
system  identification/spectral  estimation  literature.  It  implies  simply  that 
autoregressive  (AR)  models,  and  hence,  the  Yule-Walker  equation  and  the  maxi¬ 
mum  entropy  method  as  in  [2-15],  are  not  appropriate  techniques  for  spectral 
estimation  of  sampled  data.  Any  failure  of  these  methods  in  spectral  estima¬ 
tion  of  "two  sinusoids  in  white  noise"  should  first  of  all  be  attributed  to 
the  wrong  formulation  of  the  estimation  problem  rather  than  to  a  deficiency  in 

I 
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these  algorithms.  Note,  however,  that  an  autoregressive  moving  average 
ARMA  (n,m)  can  be  approximated  by  an  AR(N)  process  with  N  »  n,  and  thus, 
a  methodology  for  spectral  estimation  developed  for  AR  models  can  be  applied 
to  ARMA  models  by  considering  high  order  models. 

2.4.2  Problem  Formulation 

With  the  above  background,  we  are  ready  to  state  the  parametric  spectral 
estimation  problem. 

Main  Problem:  Assume  the  sequence  [z(kT):  k  =  0,1,..., N]  represents  or 

can  be  closely  approximated  by  Eq.  (2-2)  to  (2-3).  (See  Reference  2-2.) 

Estimate  Ec. 

The  above  problem  can  be  broken  into  two 

Subproblem  1 ;  Suppress  T.  Estimate  E^.  Denote  the  estimate  by 

Subproblem  2:  Given  (E,(T)),  evaluate  a  TCI  E  . 

-  a  c 

Our  main  concern  is  usually  in  estimating  the  following 

(1)  The  system  eigenvalues 

Ac(i)  -  -c(i)m(i)  +  j(-l)ioi(i)  ;  u>(2i)  =  u(2i  -  1)  :  i  »  1...2&J 

Ac(i)  *  -o(i)  :  i  *  2£1  +  l...*2  ;  2^  +  *2  *  n 

n  being  the  assumed  or  estimated  order. 

(2)  The  system  input  output  matrices  Bc,  Cc  where  BCBC'  *  R^.. 

Note  that  by  a  proper  transformation  to  canonical  form,  C  can  be  assumed 
known.  Furthermore,  recall  that  there  exists  B  such  tha£  p(B  )  -  p(C  ), 
which  solves  the  problem  at  hand.  c 

In  the  following  sections,  we  shall  concentrate  on  Subproblem  1  and  the 
role  of  T  in  its  solution,  and  on  Subproblem  2(1).  The  problem  of  estimating 
the  structure  of  Ec  is  beyond  the  scope  of  this  section.  We  shall  only 
consider  here  the  problem  of  spectral  estimation  under  the  assumption  of  a 
single  input /single  output  (SISO)  system  driven  by  white  noise.  In  this  con¬ 
text,  the  estimation  of  closely  spaced  modes  and  the  dependence  of  the  estima¬ 
tion  on  T  is  our  main  concern.  This  yields  the  formulation  of  the  problem  via 
Eq.  (2-11),  i.e.,  we  assume  (like  Beexe  and  Scharf  [2-2])  that  the  covariance 
function  of  the  WSS  process  measured  can  be  represented,  or  can  be  closely 
approximated,  by  a  finite  linear  combination  of  stable  exponential  functions. 
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2.5  The  Role  of  Sampling  Rate 

Assuming  the  covariance  sequence  to  be  as  defined  in  Eq.  (2-11)  for  the 
S1S0  case,  i.e. 

**  I  ip  i 

Rd(k)  "  l  6d(1)Xd(1)  k  =  0;  ±1;  ±2  (2-19) 

i-1 

we  have  the  recursive  relation 

n 

Rd(0  =  l  a(i)RdU  -  i)  £  >  n  (2-20) 

Rewriting  Eq.  (2-20)  L  times,  i.e.,  £  *  n,  n  +  1. . .n  +  L  -  1,  where 
L  ^  n,  we  obtain  in  matrix  form 

R  (O.n.L)  •  a(n)  -  R(n,L)  (2-21) 

where  R  (j,n,L)  is  the  n  x  L  Toeplitz  matrix  defined  as 


R  ( j ,n,L)  * 


R(j  +  n  -  1) 
R(j  +  n) 

R(j  +  n  +  L  -  2) 


R(j  +  n  -  2) 
R(j  +  n  -  1) 


R(j ) 


R(j  +  L  -  1) 


(2-22) 

and  the  vectors  a(n)  and  R(n,L)  are 

a(n)  -  [a(l)  ...  a(n)j'  (2-23) 

R(n,L)  -  (R(n)  ...  R(n  +  L  -  1)'  (2-24) 

Equations  (2-20)  to  (2-24)  are  the  basis  for  many  parametric  spectral 
estimation  techniques  [2-2,  2-3,  2-16  through  2-18],  Since  Eq.  (2-19)  is 
often  an  approximate  relation,  and  since  the  covariance  sequence  is  often 
obtained  from  the  measurements  of  a  sample  pass  of  the  WSS  ergodic  stochastic 
process,  a  reasonable  approach  to  the  spectral  estimation  problem  is  that  of 
solving  Eq.  (2-20)  to  (2-24)  using  a  large  L.  By  solving  the  equations,  we 
mean  solving  in  the  least-squares  sense,  where  one  allows  for  weighting  on  the 
equations  using  the  weighting  matrix  W,  i.e. 


a 


R  (0,n,L)*  •  R(n,L) 


(2-25) 


( 


R 

w 


(  R  WR  )  R  W 


(2-26) 


The  following  observations  can  be  made. 

(1)  If  the  process  is  an  AR  process  (an  unlikely  event  in  the  case  of 
sampled  data),  Eq.  (2-20)  to  (2-24)  are  valid  for  all  Jt.  Writing 
the  equations  for  l  *  l...n,  and  taking  advantage  of  the  symmetry 
R(i)  ■  R(-i),  yield  the  celebrated  Yule-Walker  equations  and  esti¬ 
mation  procedures  which  are  recursive  in  order. 

(2)  Assuming  an  ARMA  (n,m)  model,  i.e. 

n  m 

y(k)  -  l  a(i)y(k  -  i)  +  \  b(i)  «(k  -  i)  (2-27) 

i=l  i*0 

we  have  Eq.  (2-20)  to  (2-24)  satisfied  for  l  >  m.  Using  the  in¬ 
strumental  variable  (IV)  approach  to  system  identification  [2-19] 
with  [y(k  -  £ ) ;  £  =  m  +  l  ...  m  +  n]  as  the  '‘instrument",  and  as¬ 
suming  the  approximation 

1  N 

RflO  -  h  l  y(i)  yO-  +  k>  ;  y(j>  *  o  v  ]  <  o  (2-28) 

i»i 

we  have  another  way  of  justifying  Eq.  (2-20)  through  (2-26).  In 
particular,  for  the  AR  case  (m  =  0)  the  least-squares  procedure 
is  equivalent  to  the  IV  procedure,  which  is  equivalent  (with  Eq. 
(2-28))  to  the  Yule-Walker  equations. 

(3)  To  accommodate  the  inaccuracies  in  the  values  of  R(* )  used  in  Eq. 
(2-20)  through  (2-24),  one  has  to  solve  a  structured  total  least- 
squares  problem.  This  subject  is  addressed  in  Reference  2-20. 

(4)  If  { z( ) }  ,  rather  than  { y ( ) } ,  is  used  in  estimating  { R( ) } ,  R  (l,n,L) 
and  R(n  +  1,L)  should  be  used  in  Eq.  (2-25)  and  (2-26). 

Since  our  main  concern  is  using  the  sampling  rate  to  improve  reso¬ 
lution,  we  consider  here  the  case 

A  (i)T 

Ad(i)  =  e  ;  i  -  l...n  (2-29) 


where 


A  (2i  -  1) 


A  ( 2i ) 


-C.o)  +  jw 


;  i  «  1 . . .n/2 


X  denotes  the  complex  conjugate  of  X,  Ci  being  the  damping  ratio  and 
the  ic^  frequency.* 

Assuming 

T  <  TNyq  ,  (T  •  w±  <  it  V  i)  (2-30) 

we  proposed  the  following  criterion  for  choosing  T 

Definition :  The  optimal  sampling  period  T*  (T*  <  T^yq)  satisfies 

<[  R  (0,n,L)(T*)]  <  k(  R  (0,n,L)(T)]  V  T  <  T  (2-31) 

where  k [A]  is  the  condition  number  (CN)  of  the  matrix  A.  (See  References  2-21 
through  2-24  for  definition  and  properties  of  ic.) 

Note  the  explicit  dependence  of  R  on  T  induced  by  the  substitution  of 
Eq.  (2-29)  into  (2-19),  and  subsequently  into  R  as  defined  in  Eq.  (2-22). 

Rationale:  Although  the  choice  of  the  optimal  sampling  rate,  T*,  as 

defined  above,  and  the  modification  to  be  discussed  in  Section  2.6  are  based 
on  a  concept  from  numerical  analysis — the  condition  number — it  should  be 
stressed  that  the  reason  for  defining  T*  that  way  lies  in  the  domain  of  esti¬ 
mation  theory.  If  the  covariance  sequence  were  known  exactly,  if  Eq.  (2-20) 
through  (2-24)  were  satisfied  exactly,  and  if  the  computer  used  to  solve  the 
set  of  equations,  say  via  Eq.  (2-25),  was  of  infinite  precision;  the  solution, 
i.e.,  the  vector  £  and  the  corresponding  set  of  eigenvalues  [X^(i)J,  could 
be  determined  (estimated)  with  no  error.  Unfortunately,  none  of  the  above 
suppositions  are  valid,  the  computer  accuracy  being  the  lesser  evil. 

The  inaccuracy  in  the  representation  (Eq.  (2-20)  through  (2-24))  is  han¬ 
dled  rather  routinely  by  applying  the  least-squares  concept  to  a  large  number 
of  equations,  i.e.,  taking  L  large  in  Eq.  (2-20)  through  (2-24).  The  most 
severe  obstacle  to  the  spectral  estimation  is  induced  by  the  inaccuracy  of  the 
entries  to  the  matrix  R  (0,n,L)  in  Eq.  (2-22).  This  inaccuracy  stems  from 
using  the  estimated  covariance  sequence.  In  solving  a  set  of  linear  equations 
like  Eq.  (2-22),  perturbation  of  the  coefficient  matrix  (  R  (0,n,L))  has  a 
severe  consequence  on  tie  solution,  if  and  only  if  this  perturbation  is  of  the 
order  of  the  inverse  condition  number  (ICN)  of  the  matrix  [2-23].  Closely 
spaced  modes  are  shown  in  the  following  to  cause  the  matrix  to  be  numerically 
rank  degenerate  [2-23]  (i.e.,  to  have  a  small  ICN).  Hence,  our  objective  as 
stated  above  is  to  minimize  this  negative  effect. 


*  In  the  engineering  context,  the  convention  is  to  consider  the  second-order 

2  2  /  2 
equation,  X  +  2^wX  +  uj  =0,  which  gives  X  =  ±  jw/l  -  ?  ,  as  opposed  to 

the  present  convention.  For  £  «  1,  the  above  is  a  close  approximation. 


The  actual  procedure  for  solving  Eq.  (2-22),  say,  via  Eq.  (2-25), 
(2-26),  or  by  using  a  numerical  procedure  based  on  the  singular  value  decom¬ 
position  (SVD),  or  the  QR  factorization  (2-22),  is  of  secondary  Importance. 

To  demonstrate  the  preceding  discussion  and  to  establish  rules  for 
choosing  T*,  we  analyze  the  two  sinusoid  cases,  i.e.,  let 

R(t)  =  cos  WjT  +  B2  cos  i»)2t  ;  0^  +  02«1,  Bj>0 

with  the  corresponding  sampled  covariance  sequence 


R(k)  =  0^  cos  (u^T  k)  +  02  cos  (o^T  k) 


VW1 


We  thus  obtain  (n  =  4) 

R  (0,4,4)  *  R  (-3,4,4) 

where,  taking  advantage  of  the  symmetry  R(-j)  =  R(j) 

R(0)  R( 1 )  R(2) 

R  (-3  4  4)  =  R(1)  R(0)  R(1) 


R(0) 

R(  1 ) 

R(2) 

R(3) 

R(l) 

R(0) 

R(l) 

R(2) 

R(2) 

R(l) 

R(0) 

R(  1 ) 

R(3) 

R(2) 

R(l) 

R(0) 

and  Aj  is  the  state  transition  matrix  of  the  process  in  companion  form 


|a(4) . . ,a(l^j 

it  follows  from  the  definition  of  the  condition  number  that 

1  <  max  <  <(AB)  <  <(A)  •  <(B) 

Hence,  if  A  is  a  well-conditioned  matrix  (x(A)  is  ‘'small"),  AB  is  ill-condi¬ 
tioned  if  and  only  if  B  is  ill-conditioned.  Since  is  a  well-conditioned 
matrix,  it  suffices  to  consider  <[  R  (-3,4,4)].  We  prefer  to  discuss 
R  (-3,4,4),  since  we  can  derive  analytical  expressions  for  its  condition 
number. 

Some  algebra,  taking  advantage  of  the  symmetric  Toeplitz  structure  of 
R  (-3,4,4),  yields  the  1CN 


20 


1/k[  R  (-3,4,4) 


min  (X j  ,X2) 


=  1  +  |  { [  R(  1 )  -  R(3)  ]  2  +  4 [  R(  1 )  +  R(2) ]  2}  1/2 

L  -  Mil  +  .R(3).  -  *.  { [  R(  1)  -  R(3)  ] 2  +  4[R(1)  -  R(2)]2}1/2 


Assuming  T  <  ,  and  closely  spaced  modes  (u>2  -  w^)/u>2  «  1,  the  ICN  can  be 

shown  (using  standard  approximation  techniques)  to  be  of  the  order  of 
2 

[(oj2  -  (Oj )T]  .  Furthermore,  if  R  (-3,4,4)  is  rank  degeneration  [2-23]  for 

T  «  T*,  it  is  seen  that  the  numerical  rank  of  R  is  two,  i.e.,  only  one  sinusoid 
describes  the  data.  As  T  +  0,  the  numerical  rank  approaches  1.  Note  that  for 

T  ®  0,  the  perfect  predictor  is  y[(k  4-  1 ) T ]  *  y(kT)  which  corresponds  to  rank 
one  covariance  matrix.  Plots  of  the  ICN  for  few  values  of  u2>  g  ,  and 

32  are  given  in  Figures  2-3  through  2-11.  Figure  2-12  gives  the  value  of  T* 

(normalized  by  )  as  function  of  the  frequency  ratio.  We  note  that  as 

w^/w-  approaches  one,  the  dependence  of  T*  on  g  and  $2  diminishes.  Further¬ 
ed  ^  +  ui2  ^ 

more,  it  is  seen  that  - 2 -  T*  «  j  for  high  Wj/u>2.  We  propose  the  following 

"rule  of  thumb"  for  selecting  the  sampling  rate 


T  -  ir/toj  +  o>2)  (2-32) 


instead  of  T*.  It  follows  that  T  satisfies  the  following  properties  (see 
Tables  2-1  trough  2-4) 

(1)  T  is  a  good  approximation  of  T*,  except  for  <di/u>2  <  1/3. 

(2)  k[  R  (!••)]  -  k[  R  (T ) ] ,  except  for  small  frequency  ratio  and  large 
gl  (i.e.,  large  frequency  separation  and  the  signal  energy 
concentrated  at  the  lower  frequency). 

(3)  For  <  1/3,  T  poorly  approximates  T*,  but  k~^[  R  (T)]  is 

sufficiently  large  to  enable  good  estimation  of  coi,  u>2  via  Eq. 
(2-25)  and  (2-26). 
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0.068 

0.132 

0.159 

0.185 

0.221 

0.2 

0.224 

0.414 

0.502 

0.575 

0.383 

0.4 

0.358 

0.678 

0.778 

0.731 

0.387 

0.6 

0.162 

0.270 

0.293 

0.288 

0.191 

0.8 

0.037 

0.057 

0.060 

0.059 

0.040 

0.9 

0.009 

0.013 

0.014 

0.013 

0.009 

Table  2-2.  Optimal  sampling  rate,  T*. 


\  2 

0*1  /oj^\1 

0.2 

0.4 

0.5 

0.6 

0.8 

0.1 

0.954 

0.925 

0.909 

0.889 

0.817 

0.2 

0.918 

0.866 

0.833 

0.789 

0.740 

0.4 

0.675 

0.694 

0.714 

0.734 

0.743 

0.6 

0.596 

0.614 

0.625 

0.636 

0.656 

0.8 

0.544 

0.551 

0.556 

0.560 

0.568 

0.9 

0.521 

0.524 

0.526 

0.528 

0.532 

0.2  0.4  0.5  0.6  0.8 


0.947 


0.927 


0.927 


0.907 


0.934 


0.942 


0.968 


0.985 


0.957 


0.962 


0.979 


0.992 


01.00 


01.00 


01.00 


01.00 


01.00 


01.00 


0.693 


0.712 


0.885 


0.901 


0.952 


0.986 


0.292 


0.544 


0.855 


0.769 


0.8/ J 


0.918 


Table  2-4. 
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Left-hand  s 

Lope/right-har 

- 1 

id  slope  at  T* 

—4 

• 

\b? 

W1 'W2  x1 

0.2 

0.4 

0.5 

0.6 

0.8 

0.1 

0.105 

0.170 

0.187 

0.187 

0.088 

-2.83 

-3.32 

-3.25 

-3.02 

-1.89 

0.2 

0.218 

0.252 

0.179 

0.002 

1.60 

-4.70 

-4.99 

-4.53 

-3.62 

-1.59 

2.41 

5.57 

5.84 

4.18 

1.91 

0.4 

-u.681 

-1.44 

-2.17 

-3.20 

-1.86 

0.6 

1.85 

2.87 

2.93 

2.68 

i.53 

-0.510 

-0.971 

-1.15 

-1.25 

-0.991 

0.8 

0.468 

0.693 

0.721 

0.683 

0.439 

-0.159 

-0.264 

-0.289 

-0.291 

-0.212 

0.9 

0.112 

0.166 

0.173 

0.165 

0.109 

-0.042 

-0.066 

-0.070 

-0.069 

-0.048 

T*  -  T  =  */(«.  +  u,)  =  0 . 5it / ( to  )  -  0.5T  (2-33) 

12  2  Nyq 

It  ts  appropriate  at  this  point  to  recall  Kay's  result  for  estimating  a 
covariance  function  [2-10].  Kay  shows  that  T  =  T^yq/2  is  a  necessary  and 
sufficient  sampling  interval  for  obtaining  estimates  of  the  covariance 
function  with  small  error  variance.  This  result  is  a  nice  complement  to 
Eq.  (2-33). 

2 .6  Zooming  In 

2.6.1  Phasor  Representation 


Returning  now  to  Eq. 


RdOO 


(2-11),  we  may  rewrite  Eq.  (2-19) 


=  R  (kT) 
c 


n 


l 

i  =  l 


(i)T 

Bc(i)(e 


as 


i=l 


edU)  xd(i) 


(2-34) 


Equation  (2-34)  is  valid  for  any  T,  namely. 


it  is  not  restricted  to  T  <  T  , 

Nyq 


Thus,  one  may  employ  Eq.  (2-20)  through  (2-24)  to  estimate  a(T)  using  Eqs. 
(2-25)  and  (2-26),  and  then  solve  uniquely  for  {\,j(i)}.  In  this  section  we 
shall  explore  the  benefits  attainable  by  abolishing  the  constraint  T  <  Tpjyq 
and  assess  the  impact  of  this  operation  on  the  estimation  of  {xc(i)}. 


We  assume  here  that  the  covariance  sequence  |Rd(k)  =  R^(kT);  k  =  0,  1...} 
is  available  or  has  been  estimated.  Note  that  in  the  latter  case,  we  abide  by 
the  Nyquist  rate  constraint  for  sampling  the  process,  or  preferably,  abide  by 
Kay's  recommendation  for  obtaining  good  estimates  of  the  covariance  sequence, 
i.e.,  the  process  is  sampled  at  least  twice  the  Nyquist  rate. 


A  figurative  representation  of  Eq.  (2-34)  is  via  the  concept  of  phasors 
as  follows.  Define  the  it'1  phasor 

X  (i)T  (-U(i)w(i)|  +  jw (i )T 

F(i)  =  e  C  =  e  (2-35) 

where  1  n,  with  n  even  and  u>(2i)  =  -w(2i-l)  >  0,  i  =  l,2,...,n/2.  Namely, 
the  it‘1  phasor  rotates  (as  T  varies  over  R+)  at  an  angular  velocity,  w(i), 

around  a  spiral  which  collapses  at  a  rate  of  e  per  revoi_ution,  The  as¬ 

sumption  of  w(2i)  =  -w(2i-l)  is  practical,  but  not  essential  for  the  following 
discussion.  Actually,  the  clockwise/counterclockwise  rotation  adds  complexity 
to  the  solution  described  below.  The  phasor  interpretation  of  Eq.  (2-34) 
suggests  the  choice  of  T  such  that  the  n  phasors  will  be  evenly  distributed 
angle-wise,  as  long  as  the  damping  effect  does  not  dominate,  i.e.,  F(i)  is 
sufficiently  large. 
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2.6.2  Procedure 


Since  our  objective  is  the  estimation  [Xc(i)j,  we  can  formally  write 


X  (i)  -  ^inlF(i)] 

C  T 


(2-36) 


where  Eq.  (2-35)  should  be  interpreted  as 

Ta31e{|2nk  ±  /F(i )  |;  k  -  0,1,2,  ...  [T/TNyqI }  (2-37) 

'i  *  -T^tn  (|F(1)|>  <2-38) 

| F C i ) j  and  /F(i)  denote  the  magnitude  and  angle  of  F(i),  respectively, 
and  [x]  is  the  integer  value  of  x.  The  ambiguity  (nonuniqueness)  in  Eq. 
(2-36),  as  reflected  in  Eq.  (2-37),  and  thus  Eq.  (2-38),  in  the  case  that 
T  >  Tfjyq,  is  another  way  of  stating  the  Nyquist  rate  constraint.  Neverthe¬ 
less,  we  propose  to  use  T  >  T^yq  to  achieve  high  spectral  resolution.  The 
procedure  that  we  propose  is  outlined  in  the  following  steps,  which  we  call 
"zooming  in." 

Assume  the  process  consists  of  a  single  cluster  of  n/2  closely  spaced 
lightly  damped  modes  (i.e.,  ^  «  1  w(2i)=-ui(2i  -  1)  i  =  l...n/2).  The 

set  of  modes  is  resolved  as  follows: 

Step  1 :  Identify  the  center  frequency  of  the  cluster.  This  step  can  be 
accomplished  by  one  of  the  following  two  procedures: 

(1)  Plot  the  FFT  of  either  the  data  sampled  with  period  T^  <  Tpjyq 

or  of  the  sequence  [R, (k)  =  R  (kT,  )]. 

a  cl 

(2)  Solve  Eq.  (2-25)  and  (2-26)  in  conjunction  with  Eq.  (2-35)  using 
the  sequence  [Rc(kTj)j  and  n^  such  that  R  (0,n^,L)(T^)  is 

well  conditioned. 

Preferable  choice  of  Tj  is  such  that 

Ti  “c  -  f  <2*39> 

where  u >c  is  roughly  the  center  frequency  of  the  cluster. 

If  n^  =  n,  namely  the  matrix  R  (0,n,L)  is  well-conditioned,  then  the 
n/2  modes  are  spaced  sufficiently  apart  to  enable  their  resolution,  and  thus, 
Step  1  completes  the  spectral  analysis.  If  nj  <  n,  or  no  prior  knowledge  is 
available  regarding  the  number  of  modes  which  characterize  the  signal,  proceed 
to  Step  2. 
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Step  2:  Choose  T2  =  £  •  T*,  where  £  is  an  odd  integer  larger  than 
one,  and  T*  satisfies  Eq.  (2-39).  Repeat  Step  1  using  the  sequence  {R  (kT»)} 
If  the  parametric  method  is  used,  i.e.,  Eq.  (2-25),  (2-26),  and  (2-36)Cto 
(2-38)  are  solved,  repeat  this  step,  gradually  increasing  £  and  determining 
n  =  n(£)  to  be  the  largest  integer  for  which  R  (0,n,L)(T2)  has  full 
numerical  rank,  (i.e.,  is  well  conditioned). 


Stop  when  n(£)  does  not  increase  for  £  £j,  or  if  n(£j)/2  is  the 

expected  number  of  modes. 

Step  3 :  Resolve  the  ambiguity  introduced  by  Eq .  (2-36)  and  (2-37)  by 
selecting  the  set  of  modes  residing  near  the  center  frequency  identified  In 
Step  1.  Alternatively,  use  a  conventional  identification  technique  for 
T  <  Tfjyq  where  the  set  of  identified  parameters  is  the  finite  set  corre¬ 
sponding  to  the  different  solutions  of  these  equations.  (See  also  Reference 
2-27  for  a  selection  procedure.) 

We  have  earlier  stated  that  R  (0,n,L)(T)  is  least  singular  if  the  pha- 
sors  are  evenly  distributed.  This  statement,  although  correct,  should  not 
mislead  the  reader  to  believe  that  T2  should  be  so  chosen.  To  achieve  even 


distribution  of  the  phasors  we  need  T„/T„ 

r  2  Nyq 


2(0^  _  w2) 


when  we  consider  two 


sinusoids.  Thus  for  a  1%  frequency  separation,  a  choice  of  T2  =  100  T^yq 
gives  least  singularity,  but  at  the  same  time,  one  has  to  select  among  200 
possible  frequencies  in  Step  3  of  the  procedure. 


A  practical  solution  to  this  dilemma  is  to  choose  T2  as  small  as  pos¬ 
sible,  while  rendering  R  (0,n,L)(T2)  well-conditioned.  This  recommendation 
is  further  justified  when  considering  "Intermingling"  as  discussed  in  the  fol¬ 
lowing. 

The  behavior  of  the  ICN  for  R  (0,n,n)  and  R  (-n+l,n,n)  (symmetric  covar¬ 
iance  matrix)  for  two  sinusoids  as  function  of  T  is  illustrated  in  Figures 
2-13  through  2-16. 

The  above  discussion  focused  on  the  resolution  of  a  single  cluster  of 
modes.  In  the  case  of  multiclusters,  these  ideas  can  be  applied  using  one  of 
the  following  options. 

Option  A:  Increase  T2  while  preventing  cluster  "intermingling."  The 
problem  of  intermingling  is  most  simply  explained  by  considering  the  phasors 
rotation  as  in  Figure  2-17.  The  shaded  sections  A  and  B  represent  the  two 
clusters  of  modes  for  a  given  T  <  T^yq.  Thi  '  iasors  in  the  sectors  labelled 
1  rotate  clockwise  and  those  labelled  2  rotate  counterclockwise  with  T. 
Intermingling  occurs  when  two  shaded  sections  overlap.  Using  the  Information 
from  Step  1,  picking  T2,  which  gives  an  intermingling-f ree  phasor  distribu¬ 
tion,  is  a  simple  algebraic  exercise.  Note  that  for  a  large  enough  T2, 
intermingling  can  also  occur  with  a  single  cluster;  aliasing  being  an  example 
of  such  an  occurrence. 


Option  B:  Use  frequency  decimation  to  construct  two  signals,  each  rep¬ 
resenting  one  cluster  [2-13,  2-15].  Spectrum  analyze  each  signal  separately 
and  recombine. 
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Figure  2-17.  Illustration  of  two-cluster  distribution. 

2.6.3  Over  Parameterization 

In  the  quest  for  a  high  spectral  resolution  algorithm,  Cadzow  [2-3]  and 
Beex  and  Scharf  [2-2]  have  proposed  over  parameterization  as  one  tool  for 
improving  algorithms;  namely,  they  have  proposed  the  use  of  large  n  in  the 
representation  Eq.  (2-19)  to  (2-24).  For  example,  in  Reference  2-2,  n  *  7  is 
used  for  two  sinusoids  in  white  noise  (fourth-order  system).  In  Reference 
2-3,  reasonable  success  is  reported  in  resolving  two  sinusoids  with  a  0  dB 
signal-to-noise  ratio  and  7-percent  frequency  separation  using  n  «  20. 
Actually,  the  use  of  over  parameterization  is  quite  common  in  the  literature 
(e.g.,  [2-18],  [2-27],  and  references  therein);  however,  it  has  not  received 
sufficient  justification. 

Let  Eq.  (2-34)  be  a  correct  representation  of  the  covariance  sequence 
and  consider  the  over  parameterization  via  the  recursion 


R  (kT) 
c 


l  a(i)R  [(k  -  i>T] 
i-1  c 


k  >  L 


(2-40) 


where  T  <  Tfjyq,  and  L  is  an  integer  larger  than  n,  e.g.,  L  =  5  x  n.  Then, 

the  solution  for  the  vector  a  -  [a(l),  ...,  a(L)]'  via  Eq.  (2-25)  and  (2-26) 
involves  the  matrix 


R(5n  -  1) 


R  [0,5n, 5n] (T) 


(5n  x  5n) 


R 


On  the  other  hand,  considering  Step  2  in  our  zooming-in  procedure  with 
T2  *  5T,  we  have  to  solve  n  linear  equations  with  the  coefficient  matrix 


R  [0,n,n] (5T) 


R[5(n  -  1)1...  R(5) 


R(0) 


R[10(n  -  1)]...  R[5(n  +  1)]  R[5(n  -  1) ) 


(n  x  n) 


that  is,  R  [0,n,n](5T)  is  obtained  from  R  [0,5n,5n](T)  by  retaining  only  every 
fifth  column  and  fifth  row.  If  the  representation  (Eq.  (2-34))  is  merely  an 
approximation  [2-2],  it  would  be  advisable  to  retain  every  fifth  column 
of  R  [0,5n,5nJ  and  use  the  covaraince  least  squares  as  in  Eq.  (2-25)  and 
(2-26).  Now,  since  the  process  has  order  n  by  assumption,  the  numerical  rank 
of  R  [0,5n,5n](T)  is  n,  provided  the  numerical  rank  of  R  [0,n,n](5T)  is  n. 
Thus,  in  effect,  over  parameterization  is  equivalent  to  spreading  the  phasors 
to  increase  the  numerical  rank  of  the  covariance  matrix.  Using  sound 
numerical  procedure  in  solving  ji(5n)  yields  a  solution  which  corresponds  to  n 
dominant  modes  and  (5-1 )n  superfluous  modes.  This  is  exactly  Beex's  and 
Tuft's  observation  [2-2,  2-27]. 

The  advantages  of  our  zooming-in  procedure  are  rather  obvious.  We  avoid 
solving  a  very  high-dimensional  ill-conditioned  set  of  equations,  thus  elimi¬ 
nating  the  need  to  use  a  procedure  such  as  the  singular  value  decomposition 
( SVD)  used  by  Tufts  [2-18]  and  Cadzow  [2-3].  More  important,  we  have  to  solve 
only  an  nth  order  well-conditioned  polynomial  for  the  eigenvalues  compared 
to  an  Lth  order  polynomial.  The  penalty  paid  for  these  advantages  is  Step  3 
of  our  procedure. 

Finally,  we  note  that  using  l  +  1  parallel  processors,  Step  1  and  (Step 
2)  l  times  can  be  performed  simultaneously,  Step  3,  recombining  the  results  of 
the  parallel  processors.  The  architecture  of  such  a  parallel  numerical  anal¬ 
yzer  is  beyond  the  scope  of  this  section. 

2.7  Simulation  Results 


The  effect  of  sampling  within  and  beyond  the  Nyquist  rate  has  been 
extensively  tested  in  simulation.  A  few  examples  of  the  simulation  results 
are  given  here. 

The  signals  assessed  are  of  the  type 
n/2 

y(k)  -  l  A(i)  exp[-^(i)w(i)T.  (k)  ]  cos  [u(i)T  k  +  ( i)  ]  +  u(k)  (2-41) 

j  _ 1  *  *■ 


The  frequencies  are  -chosen  to  generate  a  closely  spaced  cluster.  We  specify 
for  each  cluster  the  frequency  separation  (FS)  (in  percent)  which  is  defined 

as 


(2-42) 

(u(k)}  is  a  zero  mean  white  Gaussian  noise  with  covariance,  V,  as  specified 
for  each  example.  We  use  the  equation  (see  Eq.  (2-25)  and  (2-26)) 

*  (0,;,n, )(T2)  •  a(T2)  =  R(n,n)(T2)  (2-43) 

to  estimate  a(T2) .  Subsequently,  the  characteristic  polynomial 
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(2-44) 


is  solved  to  find  (x<j(i)(T2);  i»l...n}  and  Eq.  (2-36)  to  (2-38)  are  used 
to  estimate  {c(i),  w(i)}.  As  entries  to  the  matrices  R  and  R,  we  use  the 
estimated  covariances 


RUV  “  r  l  y(kT2)y[(k  -U)T2]  ;  £  -  0...2n-l  (2-45) 

^  k=l 

where  L  is  fixed  and  T2  is  chosen  such  that  T2/Ti  is  an  integer.  The 
values  of  L,  Tj ,  and  T2  are  specified  for  each  run.  Note  that  Eq.  (2-45) 
gives  the  biased  estimate  of  the  covariance  (see  Reference  2-26  for  justifica¬ 
tion).  In  the  tabulated  results,  we  give  the  estimated  frequencies,  to(i).  If 

iXd(i)|  _<  1»  we  also  give  the  estimated  damping  ratios;  otherwise,  we  speci- 
y  the  estimate  of  | X^C i ) | •  A  *  in  the  tables  denotes  a  failure  in  estimat¬ 
ing  the  corresponding  value.  If  the  estimation  failed  for  T2  <  TN  ,  the 
optimal  order  and  the  corresponding  estimated  system  is  given  in  tlie  discus¬ 
sion. 


Example  1:  Two-sinusoid  case. 


u)(D  - 

1.5  rad/s 

w(  2)  = 

c(D  - 

0 

4(2)  = 

A(  1)  - 

1.0 

A(2)  = 

1.57  rad/s  FS  =  4.56% 

0 

1.0 
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(1)  V  -  10-/  ;  T2  <  TNyq 

The  estimated  values  are  given  in  Table  2-5.  This  example  demon¬ 
strates  our  observations  in  Section  2.2.  The  estimates  are  reason¬ 
ably  good  for  the  range  0.5Tfjyq  <  T2  <  0.9T^yq.  The  optimal 
sampling  rate  is  (see  Figure  2-11)  T*  =  0.55  T  =  1.1. 


T2  <  TNyq 


The  results  of  the  simulations  are  tabulated  in  Table  2-6. 

Here  reasonable  estimates  are  obtained  for  the  range 

0.5TN^  <  T2  <  0.8TN^.  In  case  (1),  the  optimality  of  T2  =  T* 

is  obvious.  In  case  (2),  to(l)  is  best  estimated  for  T2  *  T*;  how¬ 
ever,  a  somewhat  better  estimate  of  o>(2)  is  obtained  for  T2  =  1.3. 


Both  examples  consisted  of  undamped  sinusoids,  and  the  error  in  the 
estimates  of  | i )  |  did  not  exceed  6  percent,  and  was  usually 

considerably  smaller.  However,  attempts  to  estimate  the  damping 
ratios  for  z,  >  0  have  failed  for  T^  <  T^  1°  following  two 

examples,  we  extend  the  range  of  T2  and  demonstrate  success  in 
estimating  the  damping  ratios  and  considerable  improvement  in  the 
estimation  of  the  frequencies  compared  to  the  case  T2  <  T^yq. 

The  data  is  the  same  as  in  case  (2)  except  for  the  damping  ratios 


C(l)  *  0.01 


£ (2)  =  0.01 


(Note  that  1  percent  is  the  order  of  magnitude  of  material  damping — 
an  entity  of  interest  in  elasticity  and  structure  engineering.)  The 
results  of  the  simulations  are  summarized  in  Table  2-7. 


(4)  Same  as  (3)  with  damping  ratios 
C ( 1 )  =  0.05  i 


C(2)  =  0.05 


Results  are  given  in  Table  2-8. 

Example  1  demonstrates  the  value  of  selecting  the  sampling  rate  judi¬ 
ciously  within  the  Nyquist  rate  as  well  as  the  improvement  that  might  be  gained 
by  zooming  in.  However,  the  frequency  separation  and  the  noise  level  in 
Example  1  are  not  extreme  enough  to  cause  failure  of  the  estimation  except  for 
the  very  poor  choice  of  T2.  In  the  following  we  consider  more  severe  cases. 


Example  2 

u>(  1)  =  1  .5543  rad/s 


j ( 2)  =  1.57  rad/s 


FS  =  1  % 


Table  2-9  summarizes  the  simulation  results.  Only  one  frequency  has  been 

Identified  successfully  with  T,  <  T,  ,  where  the  optimal  order  of  the  recur- 

2  Nyq 

sion  is  two.  The  estimation  of  a  second-order  recursion  yields  the  estimates 
u)  =  1.562  and  £  =  0.12E-2.  We  see  •  increasing  benefits  with  larger  T2.  For 
T2  =  35T*  perfect  estimates  of  the  frequencies  are  observed. 

Example  3:  Three  sinusoids. 

w(l)  =  1.5543  m( 2)  =  1.57  w(3)  =  1.5857  FS  =  1% 

c  ( 1 )  =  0.001  r,(2)  =  0.001  c(3)  =  0.001 

A ( 1 )  =  1  A(2)  =  1  A( 3)  =  1 

V  =  10-4 

Results  are  summarized  in  Table  2-10.  We  have  observed  failure  for 
^2  <  3TNyq* 

The  optimal  order  of  the  recursion  is  n  =  2,  with  the  estimates  w  =  1.57 
and  r,  =  0.6E-2  when  using  T^  =  =  1.  For  T  =  25T*,  success  is  registered 

in  estimating  the  damping  ratios  and  the  frequencies. 

We  have  conducted  experiments  similar  to  the  previous  examples  for  clus¬ 
ters  of  up  to  five  frequencies  with  1-percent  frequency  separation  and  similar 
s ignal-to-noise  ratio.  The  results  were  qualitatively  identical  to  those  dis¬ 
cussed. 

We  bring  as  a  final  example  the  case  studied  by  Cadzow  (2-3]  and  others. 
This  example  culminates  a  few  difficulties.  Although  the  frequency  separation 
(FS  =  6.5  percent)  is  larger  than  considered  in  the  previous  examples,  estima¬ 
tion  is  hampered  here  by  the  disparity  in  the  modal  amplitudes  and  by  a  high 
noise  level  (0-dB  s ignal-to-noise  ratio). 

Example  4 

u(l)  =  0.4ir  rad/s  to ( 2 )  =  0.426tt  rad/s  FS  =  6.5% 

C(l>  =  0  ?( 2)  =  0 

A ( 1 )  =  4.47  A(2)  =  1.414 

V  =  1.0 

The  results  of  the  simulations  are  given  in  Table  2-11.  For  T^  =  1  the 
optimal  order  is  n  =  3  with  estimates  w  =  1.266,  £  =  0.3E-2,  and  the  discrete 
eigenvalue,  =  -0.09,  which  does  not  have  a  continuous  interpretation. 
Examination  of  the  table  reveals  the  Importance  of  choosing  T2  such  that 
T2  5  ( 2i  +  1)T*  for  some  £  (see  Step  2).  Note  that  using  the  estimate, 
u>  =  1.266,  T*  can  be  approximated  as  ti/2u  =  1.24.  The  choice  (Tj  =  1.0; 

T2  =  5)  does  not  satisfy  the  recommended  values  and  yields  poor  estimates. 

T^  =  i.2  with  T2  =  6  or  T2  =  18  approximates  this  requirement  and  gives 
correspondingly  good  estimates.  The  choice  T^  =  1  and  T2  =  16  =  13  x  1.24 
aiso  satisfies  the  recommendation  of  Step  2  of  the  zoom-in  procedure. 


Table  2-9 


2.8  General  Comments 


The  benefits  of  a  judicial  choice  of  the  sampling  rate,  as  well  as  our 
proposed  zoom-in  procedure,  ha'’e  clearly  been  demonstrated.  Zooming  in  can  be 
mechanized  in  various  ways.  Our  approach  was  to  select  the  group  of  frequen¬ 
cies  closest  to  the  center  frequency  identified  with  T2  <  T^yq. 

To  demonstrate  the  selection  procedure  for  resolving  the  ambiguity  in  the 
case  'i’2  >  Tfjyq  (Step  3  of  the  zoom-in  procedure),  consider  Example  2  with 
T2  =  35T*.  The  estimation  for  T2  <  Tfjyq  ^as  failed,  but  estimating  the 
fourth-order  data  using  a  second-order  model  gave  the  following  estimates: 
frequency  estimate  =  1.562,  damping  ratio  estimate  =  0.26E-2.  In  Table  2-12, 
a  subset  of  the  70  frequencies  around  1.562  and  the  damping  ratios  generated 
in  Step  3  for  T2  =  35Tj ,  Tj  =  1  are  given.  The  frequencies  marked  with  * 
are  the  closest  to  1.562  and  are  the  obvious  choice.  Note  that  the  estimated 
frequencies  in  this  case  are  the  exact  values  (see  Table  2-9).  This  example  is 
the  most  ambiguous  one.  In  other  examples  tested,  the  selection  among  the  dif¬ 
ferent  solutions  of  Eq.  (2-35)  was  even  clearer. 

2.9  Concluding  Remarks 

We  have  established  analytically  and  experimentally  in  this  section  the 
role  of  the  sampling  rate  in  high  resolution  spectral  analysis.  Presenting  the 
problem  of  resolution  as  a  numerical  one  enabled  us  to  establish  desirable  sam¬ 
pling  rates  for  high  resolution.  Moreover,  via  the  phasor  interpretation  of 
the  covariance  function  we  derived  a  zoom-in  procedure,  which  is  a  powerful 
spectral  analysis  tool.  Although  we  have  concentrated  on  two  procedures  for 
spectral  analysis,  the  FFT  and  the  parametric  regression  on  the  covariance 
sequence,  the  results  of  the  studies  have  ramifications  to  other  techniques. 

The  symmetric  covariance  matrix  can  be  interpreted  under  certain  conditions  as 
the  Fisher  information  matrix  [2-6].  The  inverse  condition  number  of  this  ma¬ 
trix  is  thus  an  information  measure  which  we  have  selected  for  numerical  rea¬ 
sons.  Maximizing  the  ICN  as  a  function  of  the  sampling  rate  thus  has  ramifica¬ 
tions  for  any  statistical  estimator. 

Although  the  emphasis  of  this  study  was  the  sampling  rate  for  continu¬ 
ous  signals,  the  results  subsume  the  case  of  purely  discrete  time  data.  In 
the  discrete  case,  T^  should  be  suppressed,  say  in  Eq.  (2-41);  and  thus, 
w(i)  •  T should  be  replaced  by  u)(i)e(0,ir).  Correspondingly,  T2  should  be 
replaced  in  all  the  expressions  by  the  integer,  T2/T^,  and  the  correspond¬ 
ing  results  follow  trivially. 

We  believe  that  the  algorithm  we  propose  is  numerically  superior  and  com¬ 
putationally  more  efficient  than  other  algorithms  proposed  recently  for  im¬ 
proved  spectral  resolution.  Our  success  can  be  attributed  to  the  fact  that  we 
have  addressed  the  source  of  the  difficulty.  Our  analysis  also  extends  to 
other  algorithms  as  discussed  in  Section  2.6.3. 

The  simulation  results  reported  in  Section  2.7  support  our  claims. 

It  should  be  noted  that  the  zoom-in  procedure  proposed  in  Section  2.6  is 
only  applicable  for  narrow-band  signals.  If  a  signal  spectrum  consists  of  a 
set  of  narrow  band  signals,  frequency  decimation  and  modulation  can  be  employed 


before  applying  the  proposed  procedure.  Wide  band  signals  with  no  distinct 
clusters  of  energy  cannot  be  analyzed  via  our  zoom-in  procedure,  though  our 
recommendations  in  Section  2.5  are  still  valid. 

Finally,  we  would  like  to  give  an  intuitive,  though  somewhat  ad  hoc  ex¬ 
planation,  of  the  main  features  of  the  high  resolution  spectral  analysis  pre¬ 
sented  in  this  section.  Replacing  a  very  narrow  band  signal  with  center  fre¬ 
quency  wc  by  sin(u>ct),  we  observe  that  sampling  the  signal  at  intervals 

T*,  satisfying  Eq.  (2-39),  generates  the  sequence  {sin(y  •  k)}.  This  is  the 
sequence  with  the  maximum  variation  that  can  be  generated  from  the  sinusoidal 
function.  Since  spectral  estimation  and  particularly  parametric  spectral  esti¬ 
mation  is  concerned  with  the  dynamics  of  the  signal,  observing  the  maximum 
variation  is  most  beneficial  for  its  estimation.  The  prism  is  a  classical  de¬ 
vice  for  splitting  closely  spaced  colors.  Desampling  the  covariance  function 
while  maintaining  maximum  variation  of  the  signal  as  in  the  zoom-in  procedure 
con  be  viewed  as  a  digital  prism  in  light  of  the  phasor  representation. 


Table  2-12 


Frequency 

Damping 

Frequency 

Damping 

1.0315 

0.28E-4 

1.0158 

0.36E-4 

1.1123 

0.28E-4 

1.1385 

0.36E-4 

1.2110 

0.24E-4 

1.1195 

0.31E-4 

1.3023 

0.24E-4 

1.3180 

0.31E-4 

1.3905 

0.21E-4 

1.3748 

0.27E-4 

1.4818 

0.21E-4 

1.4975 

0.27E-4 

1.5700 

0.19E-4 

1.5543 

0.24E-4 

1.6613 

0.19E-4 

1.6770 

0.24E-4 

1.7495 

0.17E-4 

1.7338 

0.22E-4 

1.8409 

0.17E-4 

1.8566 

0.22E-4 

1.9291 

0.15E-4 

1.9134 

0.20E-4 

Subset  of  estimated  frequencies  and  damping  ratios 
0)C  -  1.562,  T2  -  35,  Tj  -  1. 
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SECTION  3 


CONTROL  OF  LARGE  SPACE  STRUCTURES 
USING  ELECTROMECHANICAL  ACTUATORS 


3 . 1  Int  roduct ion 


The  published  literature  on  vibration  control  of  large  space  structures 
is  mainly  concerned  with  the  theoretical  issues  associated  with  the  design  of 
control  laws  or  control  strategies.  Numerous  suggestions  have  been  made  re¬ 
garding  the  form  of  suitable  control  laws  and  the  selection  of  control  law 
parameters.  Such  investigations  have  generally  based  the  design  and  analysis 
of  control  laws  solely  on  the  properties  of  the  flexible  structure  to  be  con¬ 
trolled  and  on  the  vibration  control  objectives.  The  important  effects  of 
physical  transducers,  namely  actuators  and  sensors,  in  achieving  the  vibration 
control  objectives  have  generally  been  ignored.  Of  course,  research  has  con¬ 
tinued  on  the  development  of  transducer  hardware  suitable  for  use  in  large 
space  structure  applications  [3-6,  3-7].  The  objective  of  this  section  is  to 
reconcile  the  work  on  vibration  control  of  large  space  structures  (where 
transducer  effects  have  generally  been  ignored)  with  the  work  on  transducers 
(where  application  to  the  particular  vibration  control  objectives  have  gen¬ 
erally  been  ignored).  Although  our  development  will  focus  solely  on  the  role 
of  actuator  transducers,  a  similar  development  likely  holds  with  regard  to  the 
role  of  sensors.  Thus,  our  main  emphasis  in  this  work  is  consideration  of  the 
role  of  actuators  as  an  important  and  necessary  part  of  the  problem  of  vibra¬ 
tion  control  of  large  space  structures.  Some  recent  work  [3-2,  3-3]  has 
examined  the  question  of  proper  placement  of  actuators,  ignoring  all  effects 
of  actuator  dynamics.  In  this  section,  the  effects  of  actuator  dynamics  are 
examined,  while  the  placement  of  actuators  is  assumed  fixed. 

The  overall  plan  of  Section  3  is  as  follows.  In  Section  3.2,  the 
important  issue  of  modelling  the  actuators  and  large  space  structure  is 
examined.  Section  3.3  takes  a  specific  viewpoint  towards  the  vibration  con¬ 
trol  problem--use  of  the  centralized  controller,  where  both  the  structural 
motion  and  the  actuator  motions  are  used  for  feedback  to  the  actuator  inputs, 
is  suggested.  An  alternative  using  a  decentralized  controller  consisting  of  a 
force  controller  and  actuator  controllers,  one  for  each  actuator,  is  con¬ 
sidered  in  Section  3.4.  Several  examples  are  presented  in  Section  3.5.  An 
indication  of  possible  extensions  and  a  summary  are  given  in  Section  3.6. 
Throughout  this  section  there  is  a  dual  development  for  each  of  two  distinct 
classes  of  electromechanical  actuators  mentioned  in  Section  3.2.  In  each  sec¬ 
tion,  use  of  a  single  actuator  is  initially  studied  in  some  detail,  and  the 
results  are  subsequently  extended  to  the  case  of  multiple  actuators. 

In  order  for  our  conclusions  to  be  as  general  as  possible,  the  detailed 
nature  of  actuator  hardware  is  not  considered.  Rather  our  interest  is  in  the 
dynamic  effect  that  the  actuators  have  on  the  large  space  structure  with  res¬ 
pect  to  the  vibration  control  objectives.  Thus  the  content  of  the  following 
sections  is  primarily  concerned  with  the  mathematical  characterization  of  the 
actuators,  the  large  space  structure,  and  the  feedback  controller.  This 


generality  allows  us  to  develop  qualitative  insight  into  the  role  of  actuators 
in  vibration  control  of  large  space  structures  for  a  wide  variety  of  physical 
actuator  hardware. 

3.2  Modelling 

3.2.1  Actuators 

Actuators  are  considered  to  represent  the  physical  devices  whereby 
forces  and  moments  are  actually  applied  to  a  large  space  structure.  They  can 
also  be  considered  as  energy  transducers  since  they  transform  electrical 
energy  into  mechanical  energy.  Actuators  can  also  be  viewed  as  power  ampli¬ 
fiers  since  they  typically  require  low  (electrical)  power  input  while  they 
generate  higher  (mechanical)  power  output.  Actuator  hardware  is  a  necessary 
ingredient  in  any  actively  controlled  system  [3-8] . 

An  actuator  can  be  characterized  in  terms  of  the  physical  nature  of  its 
input  and  output  signals.  Throughout,  an  actuator  input  will  be  considered  to 
be  a  voltage  signal.  The  nature  of  the  actuator  output  depends  on  the  way  in 
which  the  actuator  is  physically  connected  with  the  space  structure.  Although 
large  space  structures  may  consist  of  plate  and  shell  segments,  many  large 
space  structures  will  be  truss  connections  of  structural  members.  In  such  a 
case,  an  actuator  may  be  configured  to  provide  a  longitudinal  force  at  a  joint 
connection  of  several  members,  an  axial  force  along  a  member,  a  moment  at  a 
joint,  or  a  bending  moment  on  a  member.  Since  the  mathematical  characteriza¬ 
tions  of  all  of  these  different  physical  actuator  configurations  are  equiva¬ 
lent,  it  suffices  in  the  sequel  to  consider  the  output  of  an  actuator  to  be  a 
suitably  defined  generalized  force  applied  to  the  structure. 

All  electromechanical  actuators  consist  of  an  electrical  subsystem  and  a 
mechanical  subsystem  with  some  electromechanical  interaction.  The  electrical 
subsystem  is  assumed  to  be  of  a  simple  resistive  form  so  that  its  associated 
dynamics  can  be  ignored.  The  mechanical  subsystem  is  assumed  to  consist  of  an 
actuator  mass  and  spring  combination  which  are  significant  in  defining  the 
actuator  dynamics.  The  electromechanical  interaction  provides  a  force  on  the 
actuator  mass  of  electrical  origin.  These  features  are  typical  of  electro¬ 
mechanical  actuators  where  the  electromechanical  interaction  is  of  magnetic 
field  origin,  of  electric  field  origin,  or  of  piezoelectric  origin  [3-9]. 

Various  devices  have  been  mentioned  in  the  published  literature  as  feas¬ 
ible  actuators  for  vibration  control  of  large  space  structures:  reaction 
wheels,  proof  mass  actuators,  piezoelectric  actuators,  electronic  dampers, 
voice  coil  actuators,  cable  spool  actuators,  and  tendon  cable  actuators.  Each 
of  these  electromechanical  devices  falls  within  the  framework,  of  our  subse¬ 
quent  development.  Our  development  does  not  explicitly  include  the  case  of 
mass  expulsion  type  actuators,  control  moment  gyros,  or  hydraulic  or  pneumatic 
actuators.  However,  our  development  may,  with  some  modification,  be  extended 
to  such  cases. 


3.2.2  Modelling  Assumptions 

In  this  section,  the  basic  assumptions  which  hold  throughout  Section  3 
are  stated.  In  addition,  two  classes  of  actuators  are  described  for  detailed 
examination  in  the  subsequent  sections. 

Although  a  large  space  structure  can  often  be  described  by  a  distributed 
mass  model  using  partial  differential  equations,  our  development  will  make  use 
of  a  lumped  model  for  a  large  space  structure.  Such  a  model  can  often  be  ob¬ 
tained  by  finite-element  methods,  even  for  complex  large  space  structures. 

Thus  the  second-order  vector  differential  equation 


Mx  +  Dx  +  Kx  =  Bf  (3-1) 

is  chosen  to  represent  the  model  for  structural  dynamics.  Here  x  denotes  the 

•  •  • 

n-vector  of  generalized  structural  displacements,  so  that  x  and  x  are  struc¬ 
tural  velocity  and  acceleration  vectors.  The  n  x  n  mass  matrix  M  is  assumed 
to  be  symmetric  and  positive  definite.  The  n  x  n  structural  stiffness  matrix 
K  is  assumed  to  be  symmetric.  Since  our  interest  is  in  vibration  motion  rath¬ 
er  than  rigid  body  motion,  K  is  assumed  to  be  positive  definite.  The  n  x  n 
damping  matrix  D  is  often  rather  arbitrarily  chosen  since  there  is  sparse 
theory  for  structural  damping  available  for  guidance.  In  this  work,  the  damp¬ 
ing  matrix  D  is  assumed  to  be  symmetric  and  nonnegative  definite;  the  common 
choice  D  =  0  is  allowed.  The  right-hand  side  of  Eq.  (3-1)  characterizes  the 
influence  of  the  actuator  forces  on  the  structure;  the  m-vector  f  denotes  the 
generalized  actuator  force  vector,  while  B  is  an  n  x  m  dimensionless  influence 
matrix.  The  specific  form  of  the  actuator  force  vector  is  examined  in  con¬ 
siderable  detail  in  the  following  sections. 

There  are  numerous  mathematical  models  that  could  be  used  to  describe 
electromechanical  actuator  dynamics;  the  literature  on  actuator  dynamics  is 
abundant.  Unfortunately,  there  seems  to  be  no  single,  tractable  model  for 
electromechanical  actuator  dynamics  which  is  general  enough  to  characterize 
all  of  the  physical  actuators  mentioned  in  the  previous  section.  However,  it 
seems  to  be  the  case  that  each  of  the  physical  actuators  mentioned  does  fit 
into  one  of  two  distinct  classes  of  actuators.  Hence,  our  subsequent  theory 
is  developed  along  two  separate  but  parallel  lines,  corresponding  to  the  two 
classes  of  electromechanical  actuators  considered. 

The  two  electromechanical  actuator  classes  are  reaction  types  and  trans¬ 
mission  types.  Physically,  an  actuator  is  a  reaction  type  if  the  resultant 
force  on  the  structure  is  an  inertial  reaction  force;  an  actuator  is  a  trans¬ 
mission  type  if  the  resultant  force  on  the  structure  is  transmitted  through  a 
flexible  connection.  More  formally,  an  actuator  is  either  a  reaction  type  or 
a  transmission  type  according  to  the  specific  form  of  its  mathematical  model. 
In  the  subsequent  sections,  detailed  mathematical  models  are  developed  for 
each  of  these  two  classes  of  actuators. 

Specifically,  reaction  wheels  and  proof  mass  actuators  are  examples  of 
reaction  type  actuators,  while  piezoelectric  actuators,  electronic  dampers, 
voice  coil  actuators,  cable  spool  actuators,  and  tendon  cable  actuators  are 
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examples  of  transmission  type  actuators.  Although  the  literature  is  vast. 
References  3-10  through  3-12  examine  various  examples  of  reaction  type  actu¬ 
ators  for  use  in  control  of  elastic  systems.  References  3-13  through  3-13 
examine  various  examples  of  transmission  type  actuators  for  use  in  control  of 
elastic  systems.  This  classification  of  actuators  seems  natural  and  is  cer¬ 
tainly  convenient;  however,  the  specific  terminology  has  been  chosen  by  the 
author  and  has  not  been  used  elsewhere.  There  is  no  claim  that  the  classifi¬ 
cation  of  electromechanical  actuators  is  exhaustive,  but  many  examples  of 
electromechanical  actuators  known  to  the  author  do  seem  to  fit  nicely  into  one 
of  the  two  classes.  Thus,  the  classification  is  suggested  as  a  practical 
matter  whereby  mathematically  tractable  actuator  models  are  obtained. 

The  structural  dynamics,  as  defined  by  Eq.  (3-1),  are  defined  in  terms 
of  a  linear  differential  equation.  In  the  actuator  models  developed  in  the 
next  section,  the  linearity  assumption  is  maintained.  But  the  present  work 
could  serve  as  a  framework  for  development  of  nonlinear  actuator  models. 

For  clarity  in  the  presentation,  a  mathematical  model  of  a  large  space 
structure  controlled  by  a  single  reaction  type  actuator  is  first  developed. 
Next,  a  mathematical  model  of  a  large  space  structure  controlled  by  a  single 
transmission  type  actuator  is  developed.  Then  the  general  case  where  multiple 
actuators  are  used  to  control  a  large  space  structure  is  easily  handled. 


3.2.3  Single  Reaction  Type  Actuator 


In  this  section,  a  mathematical  model  is  developed  for  a  large  space 
structure,  as  described  by  Eq.  (3-1),  controlled  by  a  single  reaction  type 
actuator. 


The  basis  fur  our  development  in  this  section  is  the  schematic  diagram 
in  Figure  3-1.  The  large  space  structure  is  represented  by  a  lumped  mass 
model  of  which  a  single  mass  element  is  shown.  The  dynamics  of  the  electro¬ 
mechanical  actuator  are  characterized  by  an  actuator  mass,  ma  ,  connected  to 
the  structure  by  a  spring  and  damper  in  parallel.  A  reaction  type  actuator 
has  the  property  that  the  electromechanical  interaction  force,  fa,  appears 
as  a  reaction  pair  as  shown  in  this  figure.  It  is  claimed  that  the  simple 
schematic  represents  the  general  case  of  a  structure  controlled  by  a  single 
reaction  type  actuator. 

Based  on  Figure  3-1,  equations  of  motion  can  be  developed  using  a  La- 
grangian  approach.  Introduce  the  notation 

x  =  n  -  vector  of  generalized  structural  displacements 
T 

B  x  =  scalar  gi  m  rail  zed  displacement  of  the  structure  at  the  actuator 
1  si  at  i  .in 


z  =  sialat  i  :  i  1  1  i  sp  1  ac  cment  of  the  actuator  mass  relative  to 

r  h.  or  •  a, 


Figure  3-1.  Single  reaction  type  actuator. 

The  actuator/structure  Lagrangian  is  given  by 

1  .  T  .  l  t •  1 2  IT  1  1  1 2 

L  =  tt  x  Mx  +  -r-  m  z  +  B  x  -  r  x  Kx  —  —  k  j  z  (3-2) 

2  2  3  I  |  2  2  3  f  | 

where  M  and  K  are  n  x  n  matrices  as  in  Eq.  (3-1),  B  is  actuator  influence 
vector,  and  ma  >  0,  ka  >  0  define  the  actuator  mass  and  stiffness.  From 
Lagrange's  equations  obtain 

T 

Mx  +  m  B  (B  x  +  z)  +  Kx  =  Q  (3-3) 

3  I 

m  (BTx  +  V)  +  k  z  =  Q  (3-4) 

3  a  Z 

The  generalized  forces  are 

Q1  =  -Dx  (3-5) 

and 

Q  =  c  z  +  f  (3-6) 

l  a  a 

where  fa  represents  the  electromechanical  interaction  force  given  by 

f  =  b  u  (3-7) 

a  a 

Here,  u  is  actuator  input  voltage,  and  ba  >  0  is  actuator  servo  constant; 
ca  >  0  represents  actuator  damping  coefficient.  The  choice  of  coordinates 
with  z  defined  as  a  relative  displacement  results  in  Q2  depending  explicitly 
on  the  reaction  force  fa,  while  does  not. 

In  summary,  the  actuator /structure  model  obtained  is 

Mx  +  Dx  +  Kx  =  -B  m  (z  +  B^x)  (3-8) 

a 
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b  u 
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(3-9) 


k--  - 

►*  _  • 

•’  * 

• ' 


* 

r  ■ 

* 


,  •  •  r.  •.  • 

mlz+Rx)+cz+kz  = 
a  a  a 


Those;  equations  art:  inhert-iulv  coupled  through  the-  actuator  inertial  term. 
The  effective  actuator  f •  •  . t  on  the  structure  is  given  by 


f  =  -m  (V  +  B^x*)  (3-10) 

a 

which  can  be  viewed  as  a  force  in  reaction  to  the  acceleration  of  the  actuator 
mass.  It  is  this  view  that  motivates  the  designation  of  the  actuator  as  being 
a  reaction  type.  Note  that  Eq.  (3-8),  (3-9),  and  (3-10)  have  a  transmission 
zero  at  s  =  0.  Thus,  a  reaction  type  actuator  cannot  be  used  to  control  rigid 
body  motion;  the  actuator  cannot  excite  the  rigid  body  motion. 

There  are  two  approximations  that  are  commonly  made,  since  the  actuator/ 
structure  dynamics  are  considerably  simplified.  The  approximation  is  based  on 
the  "no  loading  assumption"  that  the  acceleration  of  the  structure  at  the 
actuator  location  is  negligible  compared  with  the  acceleration  of  the  actuator 
mass ,  i . e . 


j  BT*x|  «  |z*|  (3-11) 

The  result  is 

Approximation  Model  I 

Mx*  +  Dx  +  Kx  =  -B  m  z  (3-12) 

a 

mz  +  cz  +  kz  =  bu  (3-13) 

a  a  a  a 

This  approximation  essentially  decouples  the  actuator  motion  from  the  struc¬ 
tural  motion.  A  further  simplification  that  is  often  made  is  to  assume  that 
"actuator  inertial  terms  are  dominant"  in  that 

Ik  z|  <<  |m  z’l  (3-14) 

la  I  I  a  I 


The  result  is 

Approximation  Model  II 


Mx  +  Dx  +  Kx  =  -B  m  z  (3-15) 

a 


. .  • 

mz  +  cz  =  iu 
a  a  a 


(3-16) 


Approximation  Model  11  is  often  employed  to  characterize  reaction  wheel  dynam¬ 
ics,  since  in  that  case,  k  -  0. 

’  ’  a 
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The  preceding  approximation  models  should  be  recognized  as  specific  as¬ 
sumed  simplifications  of  the  general  model  described  by  Eq.  (3-8)  and  (3-9). 
There  is  often  no  a  priori  reason  for  accepting  the  "no  loading  assumption"; 
thus,  in  the  subsequent  developments,  the  exact  equations  (3-8)  and  (3-9)  are 
used.  Of  course,  all  our  subsequent  results  are  substantially  simplified  if 
one  of  the  above  approximations  is  employed. 

3.2.4  Single  Transmission  Type  Actuator 

A  mathematical  model  for  a  large  space  structure  controlled  by  a  single 
transmission  type  actuator  is  now  developed. 

The  basis  for  ^he  development  in  this  section  is  the  schematic  diagram 
in  Figure  3-2.  As  before,  the  large  space  structure  is  represented  by  a 
lumped  mass  model,  of  which  a  single  mass  element  is  shown.  The  dynamics  of 
the  electromechanical  actuator  are  characterized  by  an  actuator  mass,  ma, 
connected  to  the  structure  by  a  spring  with  a  damper.  A  transmission  type 
actuator  has  the  property  that  the  electromechanical  interaction  force,  fa, 
acts  only  on  the  actuator  mass  as  shown  in  this  figure.  It  is  claimed  that 
the  simple  schematic  of  Figure  3-2  represents  the  general  case  of  a  structure 
controlled  by  a  single  transmission  type  actuator. 


Figure  3-2.  Single  transmission  type  actuator. 

Based  on  Figure  3-2,  equations  of  motion  can  be  developed  using  a  La- 
grangian  approach.  Introduce  the  notation 

x  =  n  -  vector  of  generalized  structural  displacements 

T 

B  x  =  scalar  generalized  displacement  of  the  structure  at  the  actuator 
location 

z  =  scalar  generalized  displacement  of  the  actuator  mass 
The  actuator/structure  Lagrangian  is  given  by 
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where  M  and  K  are  n  x  n  matrices  as  in  Eq.  (3-1),  B  is  an  actuator  Influence 
vector,  and  ma  >  0,  ka  >  0  define  the  actuator  mass  and  stiffness.  From 
Lagrange's  equations  obtain 


Mx  +  Kx-Bk  (z-Bx)  =  Q. 

a  X1 


m  z  +  k  (z  -  B  x)  =  Q„ 
a  a  L 


The  generalized  forces  are 


and 


Qx  =  “Dx 


Q0  =  -c  z  +  f 
2  a  a 


where  fa  represents  the  electromechanical  interaction  force  given  by 


f  =  b  u 
a  a 


(3-18) 

(3-19) 

(3-20) 

(3-21) 

(3-22) 


Here  u  is  actuator  input  voltage,  ba  >  0  is  the  actuator  servo  constant,  and 
ca  >  0  is  the  actuator  damping  coefficient. 


In  summary,  the  actuator/structure  model  obtained  is 

Mx  +  Dx  +  Kx  =  Bk(z-  BTx) 

a 

m*z*+cz+k(z-Bx)  =  bu 
a  a  a  a 


(3-23) 


(3-24) 


These  equations  are  inherently  coupled  through  the  actuator  flexibility  terra. 
The  effective  actuator  force  on  the  structure  is  given  by 


f  =  k  (z  -  B  x) 
a 


(3-23) 


which  can  be  viewed  as  a  force  transmitted  from  the  actuator  mass  to  the 
structure.  It  is  this  view  that  motivates  the  designation  of  the  actuator  as 
being  a  transmission  type. 

Two  approximations  of  the  preceding  equations  are  now  mentioned;  the 
resulting  actuator/structure  dynamics  are  considerably  simplified.  The  first 
approximation  is  based  on  the  "no  loading  assumption,"  that  the  displacement 
of  the  structure  at  the  actuator  location  is  negligible  compared  with  the  dis¬ 
placement  of  the  actuator  mass,  i.e. 
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The  result  is 

Approximation  Model  III 


Mx  +  Dx  +  Kx  -  Bkz 

a 


(3-26) 


(3-27) 


mz  +  cz+kz  =  bu  (3-28) 

a  a  a  a 


This  approximation  essentially  decouples  the  actuator  motion  from  the  struc¬ 
tural  motion.  A  further  approximation  that  is  often  made  is  to  assume  that 
"actuator  flexibility  terms  are  dominant"  so  that 


m  z 

« 

Ik  z 

a  I 

!  a 

(3-29) 


The  result  is 

Approximation  Model  IV 


Mx  +-  Dx  +  Kx  =  Bkz 

a 


(3-30) 


cz+kz  =  bu 
a  a  a 


(3-31) 


The  above  approximation  models  should  be  recognized  as  specific  assumed 
simpli i  > cations  of  the  general  model  described  by  Eq.  (3-23)  and  (3-24). 

There  is  no  a  priori  reason  for  accepting  the  "no  loading  assumption"; 

thus,  in  ,e  ubsequent  developments,  the  exact  equations  (3-23)  and  (3-24) 
are  used.  Ai ’  of  our  subsequent  results  are  substantially  simplified  if  one 
of  the  above  approximations  is  employed. 

3.2.5  Multiple  Actuators 


In  this  section,  a  mathematical  model  for  a  large  space  structure 
controlled  by  multiple  actuators  is  developed.  It  is  assumed  that  there  are  a 
total  of  m  actuators,  ordered  so  that  the  first  p  actuators  are  reaction  type 
and  the  last  m  -  p  actuators  are  transmission  type.  The  notation  used  in  the 
previous  sections  is  maintained,  so  that  the  actuators/structure  Lagrangian  is 
given  by 
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The  parameters  m  >  0,  k  >0  represent  the  i  actuator  mass  and  stiffness 

th 


values,  and  B^  represents  the  i  actuator  Influence  vector,  i  = 
From  Lagrange's  equations  obtain 
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The  generalized  forces  are 
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In  summary  the  actuators/structure  model  obtained  is 

Mx+Dx+Kx  =  -  y  B  m1(z  +  B^x)  +  £  B  ki(z  -  bJx)  (3-3 
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These  coupled  equations  provide  the  basic  model  for  the  actuators/structure. 

It  is  observed  that  Eq.  (3-39)  takes  on  the  form  of  Eq.  (3-1)  where  the 
actuator  force  vector  f  =  (f^,  fm)  is  given  by 

fi  =  ^a  (zi  +  BIX)  *  ieIl 

f  =  kj  (z  -  bTx)  ,  iel  (3-42 

i  a  l  i  2 

and 

B  =  (B.  ,  ...  ,  B  )  (3-43 

1  m 

is  n  x  m  influence  matrix. 

Actuator  approximations,  as  indicated  in  the  previous  sections,  can  be 
made.  The  resulting  simplifications  are  straightforward. 

3.3  Centralized  Control  Design  Viewpoint 

In  Sections  3.3  and  3.4,  two  viewpoints  toward  control  design  for  an 
actuators/structure  interconnection  are  expressed  and  developed.  In  this  sec¬ 
tion,  the  viewpoint  of  a  centralized  controller  is  examined,  where  feedback  of 
both  the  structural  motion  and  the  actuator  motion  is  possible,  with  no  a 
priori  restriction  on  the  controller  form.  A  certain  decentralized  control 
design  viewpoint  is  taken  in  Section  3.4,  where  only  local  feedback  of  the 
actuator  motion  is  allowed.  This  leads  to  a  natural  partition  of  the  control¬ 
ler  into  a  force  controller  which  depends  on  the  feedback  of  the  structural 
motion  and  actuator  servo  controllers,  one  for  each  actuator,  which  depend  on 
feedback  of  the  local  actuator  motion  only.  This  latter  control  design  view¬ 
point  is  natural  in  the  case  where  a  force  controller  is  developed  by  com¬ 
pletely  ignoring  all  actuator  effects,  and  internal  compensation  for  each 
actuator  is  used  to  justify  that  assumption. 

Although  these  two  viewpoints  have  not  been  explicitly  stated  in  the 
published  literature.,  they  have  been  implicitly  followed  in  two  cases  where 
actuator  dynamics  have  been  examined.  In  Reference  3-16,  the  actuator  dynam¬ 
ics  were  included  as  part  of  the  plant  for  purposes  of  controller  design;  this 
is  characteristic  of  the  centralized  control  design  viewpoint.  In  Reference 
3-17  a  controller  is  developed,  ignoring  actuator  dynamics,  and  the  effects  of 
actuator  dynamics  on  the  closed  loop  is  evaluated;  this  is  characteristic  of 
the  decentralized  control  design  viewpoint. 

Consideration  of  the  two  control  design  viewpoints  does  seem  useful 
since,  depending  on  the  context,  the  role  of  the  actuators  in  the  two  cases 
does  differ.  These  two  viewpoints  are  conceptual.  A  given  controller  for  an 
actuator/structure  cannot  be  classified  as  necessarily  having  been  developed 
using  a  particular  viewpoint;  however,  the  different  viewpoints  do  lead 
naturally  to  certain  controller  configurations  that  are  discussed  in  this 
section  and  in  Section  3.4. 


3.3.1  Approach 


A  large  space  structure  with  associated  actuators  can  be  viewed  as  a 
single  plant.  This  is  a  natural  viewpoint  since  the  mathematical  equations 
which  describe  the  actuators  and  structure,  e.g. ,  Eq.  (3-39),  (3-40),  and 
(3-41),  are  inherently  coupled.  Since  there  is  no  reason  to  distinguish 
between  the  actuators  and  structure  in  defining  the  plant,  a  centralized 
controller  is  based  on  feedback  of  both  the  actuator  motion  and  the  structural 
motion.  A  schematic  of  a  closed  loop,  incorporating  a  centralized  controller, 
is  shown  in  Figure  3-3.  The  actuators  and  structure  are  indicated  separately 
for  conceptual  reasons.  The  centralized  control  design  viewpoint  is  defined 
by  the  feedback  system  in  Figure  3-3. 


Figure  3-3.  Centralized  control  representation. 

3.3.2  Single  Reaction  Type  Actuator 

Centralized  control  using  the  scheme  shown  in  Figure  3-3  is  applied  to  a 
large  space  structure  controlled  by  a  single  reaction  type  actuator.  The 
development  is  based  on  Eq.  (3-8)  and  (3-9)  for  the  actuator /structure 
combination.  They  can  be  written  in  compact  form  as 


Mq  +  Dq  +  Kq  =  Bu 

where  q  =  (x,z)  and  Mt  D,  K,  B  are  partitioned  matrices 
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It  is  important  to  recognize  that  the  mathematical  form  of  the  augmented 
Eq.  (3-44)  is  the  same  as  the  mathematical  form  of  the  large  space  structure 
Eq.  (3-1)  alone,  in  the  sense  that  the  (n  +  1)  x  (n  +  1)  matrices,  M,  D,  K, 
are  symmetric  with  M,  K  positive  definite,  and  D  nonnegative  definite.  The 
practical  implications  of  this  observation  are  great:  all  control  design  pro¬ 
cedures  developed  explicitly  for  large  space  structure  equations  are  formally 
applicable  to  the  actuator/structure  combination.  Thus,  there  is  a  wide  range 
of  design  and  analysis  procedures  which  are  directly  applicable  to  Eq.  (3- 
44).  A  few  simple  results  are  mentioned. 

To  illustrate  the  centralized  control  design  viewpoint  consider  the  use 
of  direct  velocity  feedback. 


In  the  most  general  case,  both  the  actuator  velocity  and  the  structural 
velocity  can  be  used  for  feedback  control.  In  this  case,  the  controller  is 
given  by 


u 


(3-49) 


In  this  case,  the  second  term  in  Eq.  (3-49)  is  assumed  to  have  a  positive 
coefficient  in  order  to  characterize  positive  damping.  The  constant  scalar, 
Cf ,  and  vector,  C^,  represent  the  feedback  gains.  In  this  case,  the 
closed  loop  is  described  by  equations 

M  q  +  D  q  +  K  q  =  0  (3-50) 


where 


D 

c 


(3-51) 


69 


0 

c 

a 


(3-46) 


and 


K 


T 

0 


0 

k 

a 


(3-47) 


(3-48) 


It  is  important  to  recognize  that  the  mathematical  form  of  the  augmented 
Eq.  (3-44)  is  the  same  as  the  mathematical  form  of  the  large  space  structure 
Eq.  (3-1)  alone,  in  the  sense  that  the  (n  +  1)  x  (n  +  1)  matrices,  M,  D,  K, 
are  symmetric  with  M,  K  positive  definite,  and  D  nonnegative  definite.  The 
practical  implications  of  this  observation  are  great:  all  control  design  pro¬ 
cedures  developed  explicitly  for  large  space  structure  equations  are  formally 
applicable  to  the  actuator/structure  combination.  Thus,  there  is  a  wide  range 
of  design  and  analysis  procedures  which  are  directly  applicable  to  Eq.  (3- 
44).  A  few  simple  results  are  mentioned. 

To  illustrate  the  centralized  control  design  viewpoint  consider  the  use 
of  direct  velocity  feedback. 

In  the  most  general  case,  both  the  actuator  velocity  and  the  structural 
velocity  can  be  used  for  feedback  control.  In  this  case  the  controller  is 
given  by 


u  = 


(3-49) 


In  this  case,  the  second  cerm  in  Eq.  (3-49)  is  assumed  to  have  a  positive 
coefficient  in  order  to  characterize  positive  damping.  The  constant  scalar, 
Cj,  and  vector,  C^,  represent  the  feedback  gains.  In  this  case,  the 
closed  loop  is  described  by  equations 
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where 
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with  c  =  c  +  c  .  In  this  case,  the  matrix  D  is  not  symmetric,  in  fact, 
a  f  c 

the  skew-sytnme t r ic  part  of  D  represents  the  usual  gyroscopic  type  terms 

c 

[3-18],  The  closed-loop  system  (Eg,  (3-50))  is  stable  if  and  only  if  each 
eigenvalue  of  Eq.  (. 3  —  50 )  has  negative  real  part.  A  simple  sufficient  condi 
t ion  for  the  stability  of  Eq.  (3-50 ) ,  based  on  resu  It  s  in  Reference  3-18,  i 
the  following.  Any  choice  of  feedback  gains  for  which  the  symmetric  matrix 
T 

D  ♦  D  is  positive  definite  guarantees  that  the  closed  loop  defined  by  E 


(3-50)  is  stable.  Suppose  D  is  positive  definite;  the  closed  loop  can  be 
stabilized  for  any  actuator  velocity  gain  Cf  >  0  if  the  structural  velocity 
gain,  ,  is  sufficiently  small. 

One  special  case  of  Eq.  (3-49)  of  particular  interest  is 


(3-52 


where  only  actuator  velocity  is  used  for  feedback.  In  this  case,  the  closed 
loop  can  be  viewed  as  the  feedback  connection  of  the  positive  real  transfer 
funct ion 
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and  the  feedback  gain,  — j  .  Hence,  based  on  the  results  of  Reference  3-19, 

a 

the  closed-loop  system,  Eq.  (3-50),  with  C  =  0,  is  stable  if  c,  >  0.  More- 


Vi  X- 

over,  as  shown  in  Reference  3-19,  the  closed  loop  has  certain  desirable 
robustness  properties. 

Even  the  simplest  case  of  Eq.  (3-49),  corresponding  to  open-loop  control 

u  =  0  (3-54) 

may  be  a  feasible  control  strategy.  The  closed  loop  system  (Eq.  (3-50))  with 
c,  -  0,  C,  -  0,  is  stable  if  D  is  positive  definite  and  c  >  0 .  This  simple 
control  configuration  corresponds  to  the  case  where  the  reaction  actuator  is 
used  as  a  passive  damper  [3-20]. 

The  above  example  control  strategies  are  only  illustrative  of  the 
centralized  design  approach  applied  to  a  large  space  structure  with  a  single 
reaction  type  actuator. 


3.3.3  Single  Transmission  Type  Actuator 

Consider  Eq.  (3-24)  and  (3-25)  developed  previously  to  describe  the 
dynamics  for  a  large  space  structure  controlled  by  a  single  transmission  type 
actuator.  Those  equations  can  be  written  in  the  augmented  form 
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Mq  +  Dq  +  Kq  =  Bu 

where  q  =  (x,z)  and  M,  D,  K,  B  are  partitioned  matrices 
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In  this  case,  it  is  recognized  that  the  mathematical  form  of  the  augmen¬ 
ted  Eq.  (3-55)  is  the  same  as  the  mathematical  form  of  the  large  space  struc¬ 
ture  Eq.  (3-1)  alone,  in  the  sense  that  the  (n  +  1)  x  (n  +  1)  matrices  M,  D, 

K  are  symmetric  with  M,  K  positive,  and  D  nonnegative  definite.  The  practical 
implication  of  this  observation  is  that  any  control  design  procedure  developed 
for  large  space  structure  equations  is  directly  applicable  to  the  actuator/ 
structure  combination.  Thus,  there  are  many  design  and  analysis  procedures 
that  can  be  applied  to  Eq .  (3-55).  A  few  such  results  are  now  mentioned. 

Direct  velocity  feedback  is  again  used  to  illustrate  the  centralized 
control  design  viewpoint  in  this  case. 

Assuming  both  actuator  velocity  and  structural  velocity  are  used  for 
feedback,  consider  the  controller  given  by 


(3-60) 
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where  Cf  and  Cj  are  feedback  gains.  The  closed  loop  is  described  by 
equation 


M  q  +  D£q  +  K  q  =  0 


(3-61) 


where 


(3-62) 


with  c  =  ca  +  Cf.  Comparing  Eq.  (3-60)  through  (3-62)  to  Eq.  (3-49) 
through  (3-51),  the  discussion  following  Eq.  (3-51)  applies  here  verbatim. 

3.3.4  Multiple  Actuators 

In  this  section,  centralized  control  of  a  large  space  structure  con¬ 
trolled  by  multiple  actuators  is  considered;  the  assumptions  made  in  Section 

3.2.5  are  continued  so  that  the  actuators/structure  are  described  by  Eq. 
(3-39),  (3-40),  and  (3-41).  The  augmented  variables  z  =  (zf,...,zm), 

q  =  (x,z),  and  u  =  (uf,...,um)  are  introduced  so  that  Eq.  (3-40),  (3—4 1 ) , 
and  (3-42)  can  be  written  as 


where 


l  q  - 

6  K  q 

=  B  u 

(3-63) 

(n 

+  m)  x 

(n  +  m)  partitioned  matrices 

lll 

M12 

0  ~ 

[T 

12 

M22 

0 

(3-64) 

0 

M33 

M  + 

l  m 

1b,bJ 

(3-65) 

[m^  B  ,  . . . ,  m^B  ] 
a  1  a  p 


diag  (m  ,  . . .  ,  nr  ) 
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,  P+l  m. 

diag  (ra  ,  . . . ,  m  ) 
a  a 


(3-66) 

(3-67) 

(3-68) 


i 


In  this  general  case,  the  mathematical  form  of  the  augmented  Eq.  (3-63) 
is  the  same  as  the  mathematical  form  of  the  structure  of  Eq.  (3-1)  alone  in 
the  sense  that  the  (n  +  m)  x  (n  +  m)  matrices  M,  D,  K  are  symmetric,  with  M,  K 
positive  definite,  and  D  nonnegative  definite.  Again,  the  implication  of  this 
observation  is  that  all  control  design  procedures  developed  for  large  space 
structure  equations  are  directly  applicable  to  the  actuators/structure 
comb inat ion . 

Obviously,  there  are  innumerable  possible  control  configurations  in  the 
case  of  multiple  actuators.  A  few  cases  are  considered  using  feedback  of  the 
actuator  velocities  and  the  structural  velocities. 

Consider  the  feedback  controller 


u  =  ~Ba  (C£z  +  CHX) 


(3-78) 


where  is  m  x  m  gain  matrix  and  is  m  x  n  gain  matrix. 
The  closed  loop  is  described  by 


Mq+D^q+Kq 
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where 
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The  nonsymmetric  matrix  includes  gyroscopic  terms  as  before.  The  closed-loop 
system,  Eq.  (3-79),  is  stable  if  and  only  if  each  eigenvalue  of  Eq.  (3-79)  has 
a  negative  real  part.  A  simple  sufficient  condition,  based  on  results  in  Ref¬ 
erence  3-18  is  that:  any  choice  of  feedback  gains  for  which  the  symmetric 
ma t r l x.  D  +  D  is  pocitive  definite  guarantees  that  the  closed  loop  defined  by 


by  Eq.  ( 3  —  7 9 )  is  stable.  Suppose  D  is  positive  definite.  The  closed  loop  can 
be  stabilized  for  any  symmetric,  positive  definite  actuator  velocity  gain 
matrix,  Cj,  if  the  structural  velocity  gain  matrix,  ,  is  sufficiently 


sma l 1 . 


One  special  case  of  the  control  Eq.  (3-78)  is 


u  =  -B  C cz 
a  f 


(3-81) 


where  feedback  control  depends  only  on  the  actuator  velocities.  In  this  case, 
the  closed  loop  can  be  viewed  as  the  feedback  connection  of  the  positive  real 
transfer  function  matrix 


3-19,  the  c  losed- loop ,  Eg.  (.3-63.),  with  =  0,  is  stable  if  C  u  symmetric 

and  positive  definite.  Moreover,  as  shown  in  Reference  3-19,  the  closed  loop 
has  certain  desirable  robustness  properties. 

The  case  where  all  actuators  are  used  as  passive  dampers  corresponds  to 

u  =  0 


The  closed  loop,  Eg.  (3-63),  with  C  =  0,  =  0,  is  stable  if  D  is  positive 

definite  and  C1  >  0,  i  =  1,  m. 

- a - - 

These  constant  gain  velocity  feedback  strategies  are  illustrative  of  the 
centralized  design  approach  applied  to  a  large  space  structure  controlled  by 
multiple  electromechanical  actuators. 

!  3.3.5  Comment  s 

The  centralized  control  design  approach  is  to  impose  no  a  priori  con¬ 
straints  on  the  controller  configuration.  This  approach  is  natural  when  the 
ac.uators  and  structure  are  viewed  as  a  "plant"  and  there  is  no  necessary  dis- 
|  t met  ion  between  feedback  of  the  actuators  motion  and  the  structural  motion. 

The  advantage  of  this  control  design  viewpoint  is  its  generality.  Since 
no  a  priori  control  constraints  are  imposed,  selection  of  a  particular  control 
configuration  can  be  made  with  freedom.  Furthermore,  any  particular  control¬ 
ler  design  obtained  on  the  basis  of  the  augmented  model,  e.g.,  Eq.  (3-63),  has 
expLicitLy  obtained  elosed-Loop  properties.  In  other  words,  Eq.  (3-63)  repre¬ 
sents  a  complete  model  of  the  actuators  and  structure. 

The  disadvantage  of  designing  a  controlle’-  on  the  basis  of  Eq.  (3-63) 
lies  in  the  fact  that  it  contains  n  +  m  modes,  i.e.,  2(n  +  m)  state  vari¬ 
ables.  This  increase  in  model  order,  due  to  inclusion  of  the  actuator  dynam¬ 
ics,  may  increase  the  computational  difficulty  in  achieving  a  good  control 
design.  Model  order  reduction  methods  have  Deen  commonly  employed  as  a  way  of 
achieving  a  reduced-order  Large  space  structure  model  of  tractable  order  [ 3 — , 
3-5].  Our  suggestion  is  that  it  is  preferable  to  perform  any  model  order  re¬ 
duction  on  the  basis  of  Eq.  (3-63),  which  incLudes  the  actuator  dynamics, 
rather  than  performing  order  reduction  on  the  large  space  structure  model 
alone  and  subsequently  incorporat l ng  the  actuator  dynamics. 
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3.4.1  Approach 


A  natural  approach  to  the  design  of  a  controller  is  to  develop  a  con¬ 
troller  for  the  large  space  structure  alone,  ignoring  all  actuator  dynamics. 

It  is  then  necessary  to  provide  local  compensation  for  each  actuator  as  a 
means  of  justifying  the  original  assumption.  The  controller  is  thus  separated 
into  two  parts:  a  force  controller,  typically  based  on  the  structural  dynam¬ 
ics  only,  and  an  actuator  servo  controller  for  each  actuator.  This  a  prion 
assumed  controller  configuration  is  referred  to  here  as  a  decentralized  con¬ 
troller.  A  schematic  of  a  closed  loop,  incorporating  a  decentralized  control¬ 
ler,  is  shown  in  Figu.„  3-4. 


Figure  3-4.  Decentralized  control  representation. 

Our  objective  is  to  consider  a  general  class  of  procedures  for  designing 
such  decentralized  control  schemes.  Since  much  attention  has  been  focused  on 
the  design  of  force  controllers,  this  issue  is  generally  ignored  here.  For 
purposes  of  illustration,  a  force  controller  which  depends  on  structural  vel¬ 
ocity  feedback  only  is  subsequently  employed.  Our  main  attention  is  focused 
on  design  of  the  actuator  controllers. 

This  decentralized  control  design  approach  is  conceptually  appealing. 
However,  the  equations  which  describe  the  large  space  structure  and  actuators, 
e.g.,  Eq .  (3-40),  (3-41),  and  (3-42),  are  inherently  coupled.  This  inability 
to  completely  separate  structural  dynamics  from  actuator  dynamics  makes  the 
design  of  a  decentralized  controller  with  guaranteed  closed-loop  properties 
difficult.  Hence,  some  attention  is  given  to  methods  for  evaluating  closed- 
loop  stability,  where  the  loop  is  closed  using  a  decentralized  controller. 


In  this  section,  a  decentralized  control  scheme  as  shown  in  Figure  3-4 
is  developed  for  the  case  of  a  large  space  structure  controlled  by  a  single 
reaction  type  actuator. 

Let  f^  denote  the  output  force  of  the  force  controller.  The  special 
force  controller  form 

f,  =  -c,BTx  (3-83) 

d  d 

T. 

is  assumed,  where  B  x  represents  the  structural  velocity  feedback  and  c<j  ls 
the  feedback  gain.  Development  of  the  actuator  controller  characteristics  do 
not  strongly  depend  on  the  force  controller,  so  subsequent  results  can  easily 
be  generalized  if  a  different  force  controller  form  is  assumed. 

Our  objective  now  is  to  develop  insight  into  methods  for  choosing  a 
suitable  actuator  controller.  Since  the  actuator  dynamics  cannot  be  separated 
from  the  structural  dynamics,  care  must  be  taken.  The  first  task  is  to 
determine  the  transfer  function  from  u  to  f,  based  on  Eq.  (3-8)  and  (3-9), 
where  the  actuator  force,  f,  on  the  structure  is  given  by  Eq.  (3-10).  To  this 
end,  define  the  scalar  transfer  function 

T  2  - 1 

G  (s)  =  B  (Ms  +  Ds  +  K)  B  (3-84) 

T 

which  represents  the  influence  of  the  force,  f,  on  the  displacement,  B  x,  of 
the  structure  at  the  location  of  the  actuator.  It  can  be  shown  that 

-b 

G  (s)  =  - - -  (3-83) 

a  c  s  +  k 

1  +  — — r -  +  G  (c  s  +  k  ) 

2  s  a  a 

m  s 
a 

Thus  Gg(s)  represents  the  effective  actuator  transfer  function,  where  the 
structural  loading  effects  on  the  actuator  are  taken  into  account.  An  actua¬ 
tor  controller  can  be  developed  using  this  actuator  transfer  function,  Eq. 
(3-85). 

The  actuator  controller  is  supposed  to  servo  the  actuator  so  that  the 
actuator  force,  f,  tracks  the  desired  force,  f^ .  This  is  achieved  through 
actuator  velocity  feedback  plus  a  suitable  dynamic  interface  between  f^  and 
the  actuator  input  voltage,  u.  Formally,  the  actuator  controller  should  be 
realizable  and,  as  closely  as  is  possible,  represent  the  system  inverse  of 
Gg(s).  It  is  also  desirable  for  the  actuator  controller  not  to  depend  ex¬ 
plicitly  on  properties  of  the  structure.  Such  constraints  are  severe. 


F  (  s  ) 
U  (  s  ) 


One  form  for  the  actuator  controller,  which  represents  a  realizable  ap¬ 
proximation  to  the  actuator  inverse,  is  given  by 


U(s)  =  -  i-  (1  +  —  )F  (s)  -  s  Z(s) 

b  v  m  s '  d  b 

a  a  a 


(3-86) 


where  c  =  ca  +  Cf .  This  actuator  controller  depends  on  actuator  velocity 
feedback  plus  dynamic  (integral)  forward  compensation.  The  resulting  transfer 
function  from  fj  to  f  is 


Fd(s) 


G^  (s  ) 


(3-87) 


k  +  m  s  (cs  +  k  )G 
a  a  a  s 

ms  +  cs 
a 


where  Gf(s)  represents  the  effective  actuator  controller  transfer  function. 
This  dimensionless  transfer  function  characterizes  the  total  effect  of  the 
actuator  and  actuator  controller  dynamics  on  the  structure. 

The  actuator  controller,  Eq.  (3-86),  is  a  reasonable  choice,  but  other 
actuator  controller  forms  could  be  selected.  In  Eq.  (3-86),  the  actuator  con¬ 
troller  output  voltage,  u,  depends  on  feedback  of  the  velocity  of  the  actuator 
mass;  it  is  also  proportional  to  fj  and  to  its  integral.  Clearly,  the  actu¬ 
ator  controller  described  by  Eq.  (3-86)  is  realizable.  In  the  case  that  f^ 
is  given  by  Eq.  (3-83),  the  actuator  controller  output,  u,  can  be  expressed 
directly  in  terras  of  the  structural  displacement  and  velocity,  thereby  avoid¬ 
ing  need  for  explicit  integration. 

The  effective  actuator  transfer  function,  Gf(s),  can  be  represented  as 
a  feedback  connection  of 

m  s(cs  +•  k  ) 
a  a 

2 

ms  +  cs  +  k 
a  a 


Each  of  these  transfer  functions  is  positive  real,  with  the  former  strictly 
positive  real  if  c  >  0.  Using  the  result  in  Reference  3-19,  the  actuator/ 
actuator  controller  system  defined  by  Eq.  (3-87)  is  stable  if  c  >  0.  Thus ,  a 
reaction  type  actuator  cannot  be  destabilized  when  connected  to  any  elastic 
structure  in  an  open-loop  configuration. 

The  closed-loop  configuration  of  Figure  3-4  is  now  examined,  where  the 
force  controller  is  given  by  Eq.  (3-83),  the  actuator  controller  is  given  by 
Eq.  (3-86),  and  the  actuator/ large  space  structure  is  described  by  Eq .  (3-8) 
and  (3-9).  The  closed-loop  system  is  defined  by 
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Mq  +  Dq  +  Kq 


0 


(3-88) 


M 


u 


where  q  =  (x,z)  and  M,  D,  K  are  defined  by  the  partitioned  matrices 


M 


T 

M  +  m  BB  m  B 
a  a 


T 

m  B  m 

a  a 


(3-89) 


D  = 


0 

c 


(3-90) 


K  = 


O' 

k 


(3-91) 


This  set  of  cLosed-loop  equations  does  not  have  the  symmetric  properties  of 
the  structural  Eq.  (3-1);  the  matrices  D  and  K  are  not  generally  symmetric. 

But  Eq.  (3-88)  could  form  the  basis  for  analysis  of  the  closed  loop. 

An  alternative  description  of  the  closed  loop  can  be  given  in  terms  of 
transfer  functions.  The  return  difference  function  is  given  by 

1  +  G  ( s)  G  (s)  (3-92) 

£  P 

where  the  con t rol ler /s t rue t ure  transfer  function 

Gp(s)  =  CdsBT(Ms2  +  Ds  +  K)_1B  (3-93) 

Stability  of  the  closed  Loop  is  criticaL.  The  following  conditions  are 
well  known .  The  closed  loop,  defined  by  Eq.  (3-88),  is  stable  if  and  only  if 
each  eigenvalue  of  Eq.  (3-88)  has  a  negative  real  part.  Equivalently,  the 
closed  loop  is  stable  if  and  only  if  each  zero  of  Eq.  (3-92)  has  a  negative 
real  part.  Such  necessary  and  sufficient  conditions  for  closed-loop  stability 
give  little  insight  into  the  effect  of  the  actuator  dynamics. 

From  the  viewpoint  of  decentralized  control,  the  key  issue  is  whether 
the  effective  actuator  dynamics,  defined  by  Eq.  (3-87),  destabilize  the  closed 
loop.  Recently  developed  tests  [3-17,  3-22]  are  particularly  suitable  in  this 
case  since  they  characterize  the  robustness  of  the  closed  loop  explicitly  in 
terms  of  the  effective  actuator  dynamics.  The  conditions  are  sufficient  for 
stability  of  the  closed  loop.  In  the  present  context  the  following  is  ob- 
ta  ined . 
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M 


jose  that  each  zero  of 


has  a  negative  real  part 


1  +  G  (s) 
P 


(3-94) 


G  (jw)  -  1  <  1  +  G  1  ( jco) 

i  P 


(3-93) 


holds  for  all  to  >  0 ,  the  closed  loop  is  stable. 

The  hypothesis  of  Eq.  (3-94)  is  that  the  closed  loop,  ignoring  actuator 
dynamics,  is  stable;  : s  hypothesis  is  satisfied  in  the  particular  case,  Eq. 
(3-93),  if  c^  >  0  [ 3— ' 3 ,  3-19].  The  frequency  condition,  Eq.  (3-95),  has  a 
simple  interpretation:  the  left-hand  side  represents  the  frequency  reponse 
magnitude  of  the  actuator  error,  while  the  right-hand  side  is  the  frequency 
response  magnitude  of  the  nominal  inverse  return  difference.  The  inequality 
Eq.  (3-95)  can  easily  be  checked  using  graphical  procedures.  Use  of  this 
robustness  test  is  illustrated  in  a  subsequent  example  in  Section  3.5. 

In  this  case,  as  a  partial  justification  for  the  decentralized  control 
design  approach,  consider  the  following  property.  Suppose  a  force  controller 
is  chosen  so  that  the  closed  loop,  ignoring  actuator  dynamics,  is  stable. 


There  are  actuator  parameters,  m  >  0,  c  >  0,  k  >0,  such  that  the  closed 


a  a 

loop,  including  actuator  dynamics,  is  stable.  In  other  words,  it  is  always 
possible  to  determine  an  actuator  and  actuator  controller  of  the  form  of  Eq. 
(3-86)  so  that  the  closed  loop  is  stable. 

A  few  general  remarks  can  be  made  regarding  selection  of  the  actuator 
parameters  for  a  reaction  type  actuator.  The  actuator  mass  parameter,  ma, 
should  be  selected  so  that  the  inertial  effects  of  the  actuator  are  dominant 
over  the  frequency  range  of  interest.  The  actuator  stiffness  parameter,  ka, 
should  be  chosen  as  small  as  possible. 

3.4.3  Single  Transmission  Type  Actuator 

Decentralized  control  of  a  large  space  structure  as  shown  in  Figure  3-4 
is  now  examined  where  the  structure  is  controlled  by  a  single  transmission 
type  actuator. 

A  force  controller  given  by 


T* 

-c  ,  B  x 
d 


(3-96) 


is  once  again  assumed  for  simplicity,  where  c^  is  the  feedback  gain  con¬ 
stant  . 


A  basis  for  choosing  an  actuator  controller  is  now  developed.  The 
transfer  function  for  the  actuator  relates  the  actuator  voltage,  u,  to  the 
actuator  force,  f.  Using  the  structural  transfer  function,  Gg ,  defined  by 


G  =  BT(Ms2  +  Ds  +  K)  *B 

s  ^  ‘ 


(3-97) 


and  Eq.  (3-23)  and  (3-24)  for  a  transmission  type  actuator,  the  actuator  input 
voltage,  u,  and  actuator  force,  f,  defined  by  Eq.  (3-25)  can  be  related  by 


F(s)  _  (  , 

UUJ  '  t’a(s) 


(3-98) 


m  s  +  c  s 

i  .  a  a  _  ,  2  . 

1  +  - - -  +  G  (m  s  +  c  s) 

k  s  a  a 
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Here,  Gg(s)  represents  the  effective  actuator  transfer  function,  where  the 
structural  loading  effects  on  the  actuator  are  taken  into  account.  An  actua¬ 
tor  controller  can  be  developed  using  this  actuator  transfer  function. 

As  mentioned  previously,  the  actuator  controller  should  be  chosen  to 
represent  a  realizable  approximation  to  the  actuator  inverse.  One  suitable 
form  for  the  actuator  controller,  in  this  case,  is  given  by 


"<■>  -  (i  *  f)  v») 


(3-99) 


where  c  =  ca  +  Cj.  This  actuator  controller  depends  on  actuator  velocity 
feedback  plus  dynamic  (derivative)  forward  compensation.  The  resulting  trans¬ 
fer  function  from  to  f  is 


F  ,(s) 
d 


2  2 

ms  +  (m  s  +  cs)k  G 
a _ a _ a  s 

cs  +  k 


(3-100) 


where  Gj(s)  represents  the  effective  actuator  and  actuator  controller  trans¬ 
fer  function. 

The  actuator  controller  form  Eq.  (3-99)  is  one  choice,  where  the  actua¬ 
tor  controller  output,  u,  depends  on  feedback  of  the  velocity  of  the  actuator 
mass  and  it  is  proportional  to  f^  and  its  derivative.  In  the  case  that  f^ 
is  given  by  Eq .  (3-96),  the  actuator  controller  output,  u,  can  be  expressed 
directly  in  terms  of  the  structural  velocity  and  acceleration,  thereby 
avoiding  need  for  explicit  differentiation. 


The  effective  actuator  transfer  function,  Gjis),  can  be  represented  as 
a  feedback  connection  of 


and 


k  (ra  s  +  c) 
a  a 

- 2 - 

ms  +  cs  +  k 
a  a 


sG  (  s  ) 
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Each  of  these  transfer  functions  is  positive  real,  with  the  former  strictly 
positive  real  if  c  >  0.  Using  results  from  Reference  3-19,  the  actuator/ 
actuator  controller  system  defined  by  Eq.  (3-100)  is  stable  if  c  >  0.  Conse¬ 
quently,  a  transmission  type  actuator  can  not  be  destabilized  when  connected 
to  any  elastic  structure  in  an  open-loop  configuration. 

The  closed-loop  configuration  of  Figure  3-4  is  now  examined,  where  the 
force  controller  is  given  by  Eq.  (3-96),  the  actuator  controller  is  given  by 
Eq.  (3-99),  and  the  actuator / large  space  structure  is  given  by  Eq.  (3-23)  and 
(3-24).  The  closed-loop  system  can  be  described  by 

Mq  +  Dq  +  Kq  =  0  (3-101) 

where  q  =  (x,z)  and  M,  D,  K  are  defined  by  the  partitioned  matrices 


(3-102) 
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(3-104) 


The  set  of  closed-loop  equations  (Eq.  (3-101))  does  not  have  the  symmetry 
properties  of  the  structural  equat ion  (Eq.  (3-1));  the  matrices  M,  D  are  not 
generally  symmetric.  But  Eq.  (3-101)  could  form  the  basis  for  analysis  of  the 
closed  loop. 

Another  description  of  the  closed  loop  can  be  given  using  transfer 
functions.  The  return  difference  functioq  is  given  by 

1  +  GAs)  G  (s)  (3-105) 
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.  (3-101)  has  a  negative  real  part.  Equivalently,  the  closed  loop  is 


stable  if  and  only  if  each  zero  of  Eq.  (3-105)  has  a  negative  real  part. 
Simpler  conditions  which  give  explicit  insight  into  the  effect  of  actuator 
dynamics  are  also  available  using  the  robustness  tests  of  References  3-17  and 
3-22. 

Suppose  that  each  zero  of 


1  +  G  (s) 
P 


(3-107) 


has  a  negative  real  part.  If 


G  ( jui)  -  1  <  1  +  C  1  (j w) 

f  P 


(3-108) 


holds  for  all  _>  0,  the  closed  loop  is  stable. 

As  stated  previously,  the  general  hypothesis  of  Eq .  (3-107)  is  that  the 
closed  loop,  ignoring  the  actuator  dynamics,  is  stable;  the  hypothesis  is  sat¬ 
isfied  in  the  special  case,  Eq.  (3-106),  if  c^  >  0  [3~18,  3— 1 9 J .  The  fre¬ 
quency  condition,  Eq.  (3-108),  has  a  simple  interpretation  and  can  be  easily 
checked  graphically.  Use  of  these  robustness  tests  is  illustrated  by  a  subse¬ 
quent  example  in  Section  3.5. 

As  a  partial  justification  for  the  decentralized  control  design  ap¬ 
proach,  in  this  case,  consider  the  following  property.  Suppose  a  force  con- 
troller  is  chosen  so  that  the  closed  loop,  ignoring  actuator  dynamics,  is 


stable.  There  are  actuator  parameters,  m  >  0,  c  >  0,  k  >0,  such  that  the 


closed  loop,  including  actuator  dynamics,  is  stable.  In  other  words,  it  is 
always  possible  to  determine  an  actuator  and  actuator  controller  of  the  form 
of  Eq.  (3-99)  so  tti<»t  the  closed  loop  is  stable. 

A  few  remarks  can  oe  made  regarding  selection  of  the  actuator  parameters 
for  a  transmission  type  actuator.  The  actuator  stiffness  parameter,  k  , 

Cl 

should  be  selected  so  that  flexibility  effects  of  the  actuator  are  dominant 
over  the  frequency  range  of  interest.  The  actuator  mass  parameter,  ma, 
should  be  chosen  as  small  as  possible. 

3.4.4  Multiple  Actuators 


In  this  section,  decentralized  control  of  a  large  space  structure  con¬ 
trolled  by  multiple  actuators  is  considered;  the  assumptions  made  in  Section 
3.2.5  are  continued  so  that  the  actuators/structure  are  described  by  Eq. 
(3-40),  (3-41),  and  (3-42). 


A  force  controller  is  assumed  to  be  of  the  form 


fd  =  ‘CdB  X 


(3-109) 


where 


£d  ■  (£i-  •■■■  0  <3-II0) 

B  =  (B, .....  B  ) 

1  m 

and  is  m  x  m  feedback  gain  matrix. 

Actuator  controllers  for  each  actuator  are  used  as  in  the  previous 
sections.  On  the  basis  of  the  previous  developments,  the  actuator  controllers 
are  assumed  to  be  given  by 


U.(s)  =  -  -J-r  fl  +  — )  Fj(s) - -  sZ.(s),  i el 

1  ,  i  v  l  '  d  ,ii  1 

b  ms  b 

a  a  a 


(3-111) 


U  (s)  =  ~  (l  +  — )  F*(s) - -  sZ.(s),  i el„ 

l  ,iv  ,  i  ;  d  ,ii  2 

b  k  b 

a  a  a 


(3-112) 


where  =  ca  +  °£i  1  =  . . . ,  m.  Each  actuator  controller  depends  on  actua¬ 

tor  velocity  feedback  plus  dynamic  forward  compensation.  The  resulting  m  x  m 
effective  actuator  and  actuator  controller  transfer  function  matrix  G^(s) 
defined  by 


is  given  by 


F(s)  =  G  (s)  F  ( s ) 
f  d 


(3-113) 


G  (s)  =  (l  +  q  (s)  +  Q  (s)  G  ,(s)  ) 


(3-114) 


where 


0,  (s)  =  diae  f— r- 


p+ 1  2 

m  s 
a 


(3-115) 


p+1 ,  p+1  2 

k  (ra  s  +  c  . s) 
a _ a _  p+1 

,  p+1 

c  xi s  +  k 

p+1  a 


(3-116) 


q2(») 


=  d  lag 


1  2,  1. 
ms  (c .  s  +  k  ) 
r  a  1  a 

"  1  2 

ms  +  c  ,  s 
a  1 


and 


T  2  -1 

Gs(s)  =  B  (Ms  +  Ds  +  K)  B  (3-117) 

This  actuator/a  ’'uator  controller  transfer  function,  G^  ( s ) ,  is  not  in 
general  diagonal  due  to  the  presence  of  the  general  nondiagonal  structural 
transfer  function,  Gs(s),  in  Eq .  (3-114).  In  the  case  where  multiple 
actuators  are  used  to  control  a  large  space  structure,  the  effective  actuator 
dynamics  involve  crossfeed  between  the  actuation  channels,  where  the  cross 
feed  effects  are  explicitly  due  to  structural  loading  of  the  actuators. 

The  effective  actuator  transfer  function,  Gj(s),  can  be  represented  as 
a  feedback  connection  of 


ll  +  Qj(s)J 


Q2(s> 


and 


sG  ( s ) 
s 

Each  of  these  transfer  function  matrices  are  positive  real,  with  the  former 
strictly  positive  real  if  c^  >  0,  i  =  1,  ...,  m.  From  the  results  in  Refer¬ 
ences  3-18  and  3-19,  the  actuator/actuator  controller  system  defined  by  Eq . 
(3-114)  is  stable  if  c  >  0,  i  =  1,  ...,  m.  Consequently,  multiple  actuators 

l  ■  ‘  "  — 

of  the  reaction  and  transmission  types  cannot  be  destabilized  when  connected 
to  any  elastic  structure  in  an  open-loop  configuration. 

The  closed-loop  configuration  of  Figure  3-4  is  now  examined,  where  the 
force  controller  is  given  by  Eq .  (3-109),  the  actuator  controllers  are  given 
by  Eq .  (3-111)  and  (3-11"  and  the  actuator/iarge  space  structure  is 
described  by  Eq .  (3-40),  (3-41),  and  (3-42).  The  closed-loop  system  can  be 
described  by 

Mq+Dq+Kq=0  (3-118) 

where  q  =  (x,  Zj  ,  ...,  z  )  and  M,  D,  K  are  defined  by  (n  +  m)  x  (n  +  m) 

partitioned  matrices 


l-BCd,  ....  -BCJ) 


(3-130) 


21 


and 


where 


II 

,  p+1 

(BCd  . 

m. 

Bcd) 

K11 

o  k13 

(3-131) 
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*22  ° 

(3-132) 
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r  P+ 1 
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m  •, 

.  .  .  ,  -k  B 
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(3-136) 

CT  = 

21 

4  B<v  ■ 

-^bc5] 

(3-137) 

This  set  of  closed-loop  equations  does  not  have  the  symmetry  properties  of  the 
structural  equation  (Eq.  (3-1)),  but  Eq.  (3-118)  could  form  the  basis  for 
analysis  of  the  closed  loop. 


An 

t  rans  f e  r 


alternative  .^nption  of  the  closed  loop 
functions.  The  return  difference  function 


can  be  given  in  terms  of 
is  given  by 


det  [I  +  G(s)  G  (s)] 
f  P 

where  the  cont r o  l  ler /s t rue t ure  transfer  function 

T  2  - 1 

G  (s)  =  CdsB  (Ms  +  Ds  +  k)  B 


(3-138) 


(3-139) 


j 

•  i 


'/ 


Stability  of  the  closed  loop  is  an  essential  propeLty  and  can  be  charac¬ 
terized  as  follows.  The  closed  loop,  defined  by  Eg.  (3-118),  is  stable  if  and 
only  if  each  eigenvalue  of  Eq.  (3-118)  has  a  negative  real  part.  Equivalent¬ 
ly,  the  closed  loop  is  stable  if  and  only  if  each  zero  of  Eq.  (3-138)  has  a 
negative  real  part.  Such  necessary  and  sufficient  conditions  for  closed-loop 
stability  are  so  complicated  that  little  insight  into  the  effect  of  the  actua¬ 
tor  dynamics  is  obtain’d. 

Simpler  conditions,  which  do  give  explicit  insight  into  the  effect  of 
actuator  dynamics  on  c losed- loop  stability,  can  be  developed  using  the  results 
of  References  3-17  and  3-22.  Suppose  that  each  zero  of 

de t [ 1  +  G  (s)]  (3-140) 

P 

has  a  negative  reai  part.  If 

o[G,(ju)  -  1]  <  o[l  +  G_1  ( j  to )  )  (3-141) 

t  —  p 

holds  for  all  w  >  0  the  closed  loop  is  stable.  Here  o[G]  and  o[G]  denote  the 
maximum  and  minimum  singular  values  of  a  matrix  G,  respectively. 

The  general  hypothesis  of  Eq .  (3-140)  is  that  the  closed  loop,  ignoring 
all  actuator  dynamics,  is  stable;  the  hypothesis  is  satisfied  in  the  special 
case  in  Eq.  (3-139)  if  is  symmetric  and  positive  definite  (Eq.  [3-18]). 

The  frequency  condition,  Eq.  (3-141),  has  a  simple  interpretation  and  can  be 
checked  graphically. 

In  the  case  where  multiple  actuators  are  used  to  control  a  large  space 
structure,  it  can  be  shown  that  it  is  always  possible  to  determine  actuators 
and  actuator  controllers  of  the  form  of  Eq.  (3-111)  and  ( 3— 1 i 2 )  so  that  the 
closed  loop  is  stable.  Suppose  a  force  controller  is  chosen  so  that  the 
closed  loop,  ignoring  actuator  dynamics,  is  stable.  There  are  actuator  param- 

l  i 

eters,  >  0 ,  c  ^  >  0,  k^  >  0,  l  =  1,  ...,  m,  such  that  the  closed  loop,  in¬ 
cluding  actuator  dynamics,  is  stable.  Explicit  guidelines  for  selection  of 
the  actuator  parameters  are  difficult  to  obtain.  As  in  the  previous  sections, 

i  i 

it  appears  desirable  to  chouse  large  in  >0  and  small  k  >0,  i  =  1 . .  for 

3  3 

i  l 

the  reaction  actuator  parameters  and  small  m  >0  and  large  k  >  0,  i  -  p  +  1 , 

3  3 

...,  m  fur  the  transmission  actuator  parameters.  Substantial  trial  and  error 
may  be  required  to  select  the  actuator  parameters  effectively  in  a  particular 
case . 


3.4.5  Comment  s 

The  decentralized  control  design  approach  is  to  impose  the  a  priori  con 
strain!  that  the  controller  consists  of  a  force  controller  and  a  set  of  actua 
tor  controllers,  one  for  each  actuator,  as  shown  in  Figure  3-4.  This  appi  >a. 
is  natural  when  the  structure  is  viewed  as  the  "plant"  and  the  aituatoi 


dynamics  are  ignored;  a  controller  obtained  on  the  basis  of  this  assumption  is 
essentially  a  force  controller.  The  actuator  controllers  can  be  subsequently 
developed. 


The  main  advantage  of  this  control  design  viewpoint  is  that,  since  all 
actuator  dynamics  are  ignored  in  developing  a  force  controller,  there  is  a 
reduction  in  dimensionality  and  hence  complexity  of  the  control  design  prob¬ 
lem.  Such  an  order  of  reduction  is  often  desirable  as  a  means  of  obtaining  a 
computationally  tractable  control  design  problem.  Specifically,  this  decen¬ 
tralized  approach  depends  on  obtaining  a  force  controller  based  on  n  struc¬ 
tural  modes,  plus  m  single-loop  actuator  controllers,  as  opposed  to  the  use  of 
n  +  m  modes,  i.e.,  2(n  +  m)  state  variables  if  the  centralized  control  design 
approach  is  taken. 

The  major  disadvantage  with  this  approach  is  that  there  can  be  no  guar¬ 
antee  of  specific  closed-loop  properties,  e.g.,  stability.  As  shown  in  Refer¬ 
ence  3-23  and  elsewhere,  there  is  the  possibility  that  a  closed-loop  system 
with  controllers  chosen  to  stabilize  the  closed  loop,  ignoring  actuator  dyn¬ 
amics,  may,  in  fact,  be  destabilized  by  actuator  dynamics.  Thus  closed-loop 
characteristics,  where  the  controller  is  obtained  using  the  decentralized 
design  approach,  should  always  be  carefully  analyzed.  As  indicated,  a  com¬ 
plete  eigenvalue  analysis  of  the  closed  loop  can  be  performed  using,  for 
example,  Eq.  (3-118).  An  alternative  that  often  gives  insight  into  the  spe¬ 
cific  importance  of  the  actuator  parameters  is  to  make  use  of  the  robustness 
test  as  discussed. 

3.5  Examples 


3.5.1  Control  of  Single-Mode  Structure  Using  Single  Reaction  Type  Actuator 

Our  objective  here  is  to  illustrate  the  theory  developed  for  reaction 
type  actuators  by  examining  in  some  detail  the  simplest  case  of  a  single-mode 
structure  controlled  by  a  single  reaction  type  actuator. 

From  Eq.  (3-8)  and  (3-9),  the  mathematical  model  for  the  structure/ 
actuator  is 


Mx  +  Dx  +  Kx  =  -m  (x  +  V)  (3-142) 

a 

m  (z*  +  x)  +  c  z  +  k  z  =  b  u  (3-143) 

a  a  a  a 

where  M  >  0,  D  >  0,  and  K  >  0  are  structural  mass,  damping,  and  stiffness 
parameters,  and  ma  >  0,  ca  >  0,  ka  >  0,  and  ba  >  0  are  actuator  mass, 
damping,  stiffness,  and  input  parameters;  x  denotes  the  scalar  displacement 
of  the  structure. 


It  is  convenient  to  introduce  the  parameters 


P 


K 

M* 


2  £  w 
s  s 


(3-144) 


a 

in  the  analysis  to  follow.  Several  different  control  laws  are  now  examined. 

First,  consider  feedback  control  using  only  actuator  velocity  feedback 
as  given  by 


If  the  definition 


2;  to 
a  a 


c  + 
a 


m 

a 


(3-145) 


(3-146) 


is  used,  the  closed-loop  characteristic  equation  can  be  written  as 


2^  <jj  s  [(1  +  p)s  +  lx,  u)  s  +  u 

l  *  — a— - 2 - 2 - V - FT  =  0  0-147) 

(s  +  2c  0)  s  +  oj  )(s  +  a)  )  +  ps  to 

It  is  easily  shown  that  the  closed-loop  system  defined  by  Eq.  (3-142), 

(3-143),  and  (3-144),  is  always  stable  for  any  Ca  ^  Cg  —  Thus*  the 

structure  can  be  stabilized  by  using  only  actuator  velocity  feedback  or  by 
using  the  actuator  as  a  passive  damper.  A  root  locus  plot  for  Eq.  (3-147)  for 
the  indicated  parameter  values  is  shown  in  Figure  3-5;  dependence  of  the 
closed-loop  poles  on  the  actuator  damping  parameter,  ca»  indicated.  Al¬ 
though  the  structure  is  stabilized  by  the  actuator,  there  is  a  modest  amount 
of  closed-loop  damping  that  can  be  achieved  using  the  feedback  control  Eq. 
(3-145).  The  maximum  damping  ratio  for  the  dominant  pole  pair  is  0.12,  cor¬ 
responding  to  the  selection  of  ca  =  0.68. 


Consequently,  the  centralized  control  viewpoint  suggests  use  of  the 
controller 


u  =  -^z+r^x  (3-148) 

0  0 
a  a 

where  both  the  actuator  velocity  and  the  structural  velocity  are  fed  back.  If 
the  definition 


lx,  u> 
c  s 


c 

c 


M 


is  used,  the  closed-loop  characteristic  equation  can  be  written  as 


(3-149) 


(2  5  0)  s)  82 

1  +  - - - - £-5 - - - - - —  -  o 

(s  +  2cswgs  +  oj|)(s  +  2cauas  +  u>a)  +  p  s  (25awas  +  </) 

(3-150) 

The  closed-loop  system,  defined  by  Eq.  (3-142),  (3-143),  and  (3-148),  is 
stable  for  any  ca  >  0  and  for  cc,  not  exceeding  some  bound.  A  root  locus 
plot  for  Eq.  (3-150)  for  the  particular  value  5a  =*  0.68  obtained  in  Figure 
3-5  is  shown  in  Figure  3-6;  dependence  of  the  closed-loop  poles  on  the  struc¬ 
tural  damping  parameter,  £c,  is  indicated.  In  this  case,  significant  ad¬ 
ditional  damping  can  oe  obtained  by  using  the  feedback  control  Eq.  (3-148). 
Note  that  Figure  3-o  indicates  that  there  is  a  limit  to  the  closed-loop  damp¬ 
ing  that  can  be  achieved;  the  maximum  damping  ratio  for  the  dominant  pole  pair 
is  0.25  corresponding  to  the  choice  *;c  *•  0.19.  If  the  structural  velocity 
feedback  gain  is  too  large  so  that  Cc  >  1.0,  the  closed  loop  is  unstable. 

The  decentralized  control  viewpoint  suggests  use  of  the  force  controller 


fd  -  -cdi  (3-151) 

together  with  the  actuator  controller  defined  by 

u(s)  -  -  ~  (l  +  ~  F  8Z(8)  (3"152) 

a  a  a 

where  c  -  Cf .<>  As  mentioned  previously,  these  expressions  can  be  com¬ 

bined  to  obtain  the  expression 


Cf  CA  CCA 

»  ■  T’-r^bT"  »-153> 

a  a  a 

which  can  be  compared  with  Eq.  (3-145)  and  (3-148);  an  extra  displacement 
feedback  term  appears  in  Eq.  (3-153).  Hence,  the  actuator  controller  Eq. 
(3-152)  in  conjunction  w*th  structural  velocity  feedback  for  the  force  con¬ 
troller  is  equivalent  to  the  controller  Eq.  (3-153).  The  effective  actuator 
transfer  function  using  Eq.  (3-87)  is 


Gf(s) 


(s2  +  2;au)as)(s2  +  2<;su)ss  +  oj2) 

(s2  +  2cgwgs  +  u>g)(s2  +  2?au>as  +  tu2)  +  ps2(2cau)as  +  u>2) 


(3-154) 
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Figure  3-6.  Closed-  oop  root  locus. 

Type:  reaction  actuator 

Feedback:  centralized 

Parameter:  c 


Using  the  control  Eq.  (3-153),  if  the  definitions 


2cdus  "  iT 


(3-155) 


8 

G  (s)  -  -3 - — - j 

P  2  ,  _  .2 

8  +  tt  8  +  W 

s  s  s 


(3-156) 


are  introduced,  the  ^losed-loop  characteristic  equation  can  be  written  as 


1  +  GAa)  G  (s) 
f  P 


(3-157) 


The  closed-loop  system  defined  by  Eq.  (3-142),  (3-143),  (3-151),  and  (3-152) 
is  stable  for  all  values  z,a  >  0  and  ^  >  0.  A  root  locus  plot  of  Eq. 

3-157),  for  the  particular  value  of  £a  “  0.68  obtained  in  Figure  3-5,  is 
shown  in  Figure  3-7.  Dependence  of  the  closed-loop  poles  on  the  structural 
damping  parameter,  is  indicated.  In  this  case,  significant  additional 
damping  can  be  obtained  by  using  the  force  controller  Eq.  (3-151)  and  the 
actuator  controller  Eq.  (3-152).  Figure  3-7  indicates  that  there  is  a  limit 
to  the  closed-loop  damping  that  can  be  achieved.  For  example,  if  q  *  0.25, 
the  damping  ratios  for  the  two  pole  pairs  are  0.47  and  0.48. 

As  an  alternative  to  a  direct  analysis  of  the  closed-loop  poles,  the 
robustness  test  of  Eq.  (3-95)  can  be  used.  One  procedure  would  be  to  plot  the 
frequency  response  of  the  actuator  error,  |Gf(jw)  -  1  ,  which  does  not  de¬ 
pend  on  the  force  controller.  Then  the  force  controller  could  be  adjusted  to 
guarantee  satisfaction  of  the  robustness  inequality  Eq.  (3-95).  In  this  par¬ 
ticular  example,  the  frequency  response  for  the  actuator  error  function, 
|Gf(ju>)  -l|,  is  indicated  in  Figure  3-8;  and  for  any  Cj  >  0*  the  robust¬ 
ness  inequality  Eq.  (3-95)  is  satisfied.  To  illustrate  the  procedure,  the 

frequency  response  for  the  inverse  return  difference  function,  |l  +  G  *(jaj)|. 


is  indicated  in  Figure  3-8  for  the  value  “  0.25.  Certain  conclusions  can 
be  drawn  from  the  figure  based  on  the  robustness  test,  Eq.  (3-95).  First,  it 
is  clear  that  for  5 j  ■  P.25,  and  in  fact  for  any  5^  >  0,  the  actuator  dy¬ 
namics  do  not  destabilize  the  closed  loop  using  the  force  controller  Eq. 

(3-151).  Second,  since  the  separation  between  |l  -+•  G^  ^  ( j  u>)  |  and  |  ( j  uj)  -  lj 

is  an  Indication  of  the  stability  margin.  Figure  3-8  can  be  used  to  assess 
qualitative  characteristics  of  the  closed  loop.  In  particular,  the  closed 
loop  would  be  most  sensitive  to  disturbances  at  a  frequency  near  2.0  rad/s, 
where  the  "loading  peak"  for  the  actuator  error  occurs.  Further,  the  high 
frequency  margins  are  large  since  the  actuator  error  is  small.  At  low  fre¬ 
quencies,  the  actuator  errors  are  relatively  large,  but  good  margin  is  main¬ 
tained  due  to  the  force  controller  characteristic.  These  qualitative  conclu¬ 
sions  should  be  typical  of  a  large  space  structure  controlled  by  a  single 
reaction  type  actuator,  using  a  properly  chosen  force  controller. 
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Experience  has  indicated  that  the  magnitude  of  the  loading  peak  for  the  actu¬ 
ator  error  is  strongly  affected  by  the  assumption  regarding  open-loop  struc¬ 
tural  damping.  In  particular,  the  robustness  inequality  of  Eq.  (3-95)  may  not 
be  satisfied  if  the  open-loop  structural  damping  is  assumed  too  low. 

3.5.2  Control  of  a  Single-Mode  Structure  Using  Single  Transmission  Type 
Actuator 

The  previously  developed  theory  is  now  illustrated  by  examining  in  some 
detail  the  simplest  case  of  a  single-mode  structure  controlled  by  a  single 
transmission  type  actuator. 

The  mathematical  model  for  the  actuator/structure  is  based  on  Eq. 

(3-23)  and  (3-24),  wl  cn  in  this  case  can  be  written  as 


Mx  +  Dx  +  Kx  -  k  (z  -  x)  (3-158) 

a 

m  z*  +  c  z  +  k  (z  -  x)  =  bu  (3-159) 

a  a  a  a 

where  M  >  0,  D  >  0,  and  K  >  0  are  structural  mass,  damping,  and  stiffness 
parameters  and  ma  >  0,  cg  >  0,  ka  >  0,  and  bfl  >  0  are  actuator  mass, 
damping,  stiffness,  and  input  parameters;  and  x  denotes  the  scalar  displace¬ 
ment  of  the  structure. 


It  is  convenient  to  introduce  the  parameters 
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(3-160) 


in  this  case.  Several  different  control  laws,  each  based  on  some  form  of 
velocity  feedback,  are  now  examined. 

First,  consider  feedback  control  using  only  actuator  velocity  feedback 
as  given  by 


u  =  -~z 

a 


(3-161) 


Define  thu  actuator  damping  parameter  c  from 


2c  to 
a  a 


c  +  c 

a  f 

m 


(3-162) 


The  closed-loop  characteristic  equation  can  be  written  as 


98 


0 


(3-163) 


2?aV  t®2  +  2W  +  (1  +  p)tofl 

1  +  (  2  7  7  2w  2  ^  ITT  2  2 

Is  +  2W  “s^8  +  “aJ  +  P  8  “fl 

Using  this  characteristic  equation,  it  is  easily  shown  that  the  closed-loop 
system  defined  by  Eq.  (3-158),  (3-159),  and  (3-161)  is  always  stable  for  any 
5a  >  0.  The  structure  can  always  be  stabilized  by  using  only  actuator  vel¬ 
ocity  feedback  or  by  using  the  actuator  as  a  passive  damper.  A  root  locus 
plot  for  Eq.  (3-163)  for  a  particular  case  is  shown  in  Figure  3-9;  dependence 
of  the  closed-loop  poles  on  the  actuator  damping  parameter,  ^a,  is  indicat¬ 
ed.  The  structure  can  be  stabilized  in  this  manner,  but  only  modest  damping 
of  the  closed  loop  ran  be  achieved  using  feedback  of  the  actuator  velocity 
only.  The  maximum  damping  ratio  for  the  dominant  pole  pair  is  0.12,  corres¬ 
ponding  to  the  selection  of  £a  =  0.78. 

In  order  to  achieve  additional  structural  damping,  the  centralized  con¬ 
trol  viewpoint  suggests  feedback  of  both  actuator  and  structural  velocity  as 
given  by  the  controller 


(3-164) 


Define  the  structural  damping  parameter,  £c,  from 


2?  ID 

s 


(3-165) 


the  closed-loop  characteristic  equation  can  be  written  as 


(2^cU,s)  “a 

2  2  2  2  2  2 
fs  +  2?  w  s  +  a)  )fs  +  2s  w  s  +  id  1  +  po)  (s  +  2c  us) 
^  ^ss  s  •' v  a  a 1  sv  a  a  ' 


(3-166) 


The  closed-loop  system,  defined  by  Eq.  (3-158),  (3-159),  and  (3-164),  is 
stable  for  any  £a  >  0  ?  .d  for  cc  not  exceeding  some  bound.  A  root  locus 
plot  for  Eq.  (3-166)  for  the  particular  value  ?a  *  0.78  obtained  in  Figure 
3-9  is  shown  in  Figure  3-10;  dependence  of  the  closed-loop  poles  on  the  struc¬ 
tural  damping  parameter  5C  is  Indicated.  It  is  clear  from  Figure  3-10  that 
substantial  additional  closed-loop  damping  can  be  achieved  by  feedback  of  both 
actuator  and  structural  velocities.  In  particular,  a  maximum  damping  ratio  of 
0.21  can  be  achieved  for  the  dominant  pole  pair,  corresponding  to  £;c  =  0.21. 

If  the  structural  velocity  feedback  gain  is  too  large  so  that  £c  >  1.3,  the 
closed  loop  is  unstable. 


Closed-loop  root  locus 
Type:  transmission  actuator 

Feedback:  actuator  only 

Parameter:  r 


IMAGINARY 


w 


The  decenCraltzed  control  viewpoint  suggests  use  of  the  force  controller 


fd  =  "cdx 


together  with  the  actuator  controller  defined  by 


(3-167) 


UCs)  -  A  •  ,1  +?•)  F.(.)  -  rp  sZCs) 


(3-168) 


where  c  =  ca  +  Cf.  The  force  controller  and  the  actuator  controller  can 
be  combined  to  obtain  •’he  equivalent  control  expression 


cf  •  cd  •  C  Cc 

u  =  _  b~  z  ~  b-  x  ir- 

a  a  a 


(3-169) 


which  can  be  compared  with  Eq.  (3-161)  and  (3-164);  this  form  for  the  control¬ 
ler  involves  an  extra  acceleration  feedback  term.  The  effective  actuator 
transfer  function  based  on  Eq.  (3-100)  is  given  by 


G^(s)  - 


_ (2W  +  Ma^s2  *  2csV  *  -s2) _ 

(s2  +  24swss  +  “g  )(s2  +  2^awas  +  %)  +  pw2(s2  +  2cau>as) 


(3-170) 


For  the  control  Eq.  (3-169),  define  the  structural  damping  parameter,  by 


2;dWs  “  M 


(3-171) 


and  define 


2c.u>  s 

Of.)  -  ~2 - ^ - y 

‘  s  +  2;;  u  s  +  iii 

s  s  s 


(3-172) 


The  closed-loop  characteristic  equation  is 


1  +  G,(s)  G  (s)  =  0 

f  P 


(3-173) 


The  closed-loop  system,  defined  by  Eq.  (3-158),  (3-159),  (3-167),  and  (3-168), 
is  stable  for  any  ca  >  0  and  any  £<j  0.  A  root  locus  plot  of  Eq.  (3-173) 

for  the  particular  value  of  =  0.78  obtained  in  Figure  3-9  is  shown  in 
Figure  3—11;  dependence  of  the  closed-loop  poles  on  the  structural  damping 
parameter  is  indicated.  It  is  clear  that  significant  damping  can  be  ob¬ 
tained  using  the  force  control  Eq.  (3-167)  and  the  actuator  control  Eq.  (3-168). 
As  seen  from  Figure  3-11,  there  is  a  limit  to  the  closed-loop  damping  that  can 


be  achieved.  If  *  6.32,  the  damping  ratios  for  the  two  pole  pairs  are 
0.50  and  0.58. 


Instead  of  a  direct  analysis  of  the  closed-loop  poles,  an  analysis  of 
the  closed  loop  using  the  robustness  test,  Eq.  (3-108),  can  be  used.  In  this 
particular  example,  the  frequency  response  for  the  actuator  error  function. 


|Gf(j  w)  -  1 1 ,  is  indicated  in  Figure  3-12;  for  any  Cd  ^  ®  the  robustness 


inequality  of  Eq.  (3-108)  is  satisfied. 


As  an  illustration,  the  frequency 


response  for  the  inverse  return  difference  function,  J 1  —  G  1(jw)j,  is  also 

indicated  in  Figure  3-12  for  the  value  Cd  *  0.32.  Certain  conclusions  can 
be  drawn  from  Figure  3-12,  based  on  the  robustness  test  in  Eq.  (3-108). 

First,  for  the  specify  case  where  =  0.32  as  shown  in  Figure  3-12,  and, 

in  fact,  for  any  ^  ^  0,  the  actuator  dynamics  do  not  destabilize  the  closed 
loop  using  the  force  controller  Eq.  (3-167).  Also,  the  closed  loop  is  most 
sensitive  to  disturbances  with  a  frequency  near  1.0  rad/s,  where  the  "loading 
peak"  for  the  actuator  error  occurs.  The  low-frequency  margins  are  large 
since  the  actuator  error  is  small.  At  high  frequencies,  the  actuator  errors 
are  relatively  large  but  a  good  margin  is  maintained  due  to  the  force  control¬ 
ler  characteristic.  These  qualitative  conclusions  should  be  typical  of  a 
large  space  structure  controlled  by  a  single  transmission  type  actuator  using 
a  properly  chosen  force  controller.  Experience  has  indicated  that  the  magni¬ 
tude  of  the  loading  peak  for  the  actuator  error  is  strongly  affected  by  the 
assumption  regarding  open-loop  structural  damping;  the  robustness  inequality 
of  Eq.  (3-108)  may  not  be  satisfied  if  the  open-loop  structural  damping  is 
assumed  too  low. 


3.5.3  Comments 


The  two  simplest  possible  examples  have  been  considered,  where  the  large 
space  structure  can  be  described  by  a  single  vibration  mode.  Caution  must  be 
exercised  in  making  general  conclusions  on  the  basis  of  two  examples  but  the 
following  comments  are  suggested  by  the  examples. 

It  is  possible  to  add  damping  to  a  structure  by  operating  an  electro¬ 
mechanical  actuator  as  a  passive  damper,  or  by  using  actuator  velocity  feed¬ 
back  only.  It  appears  that  this  control  approach  can  add  only  a  modest  amount 
of  damping  to  a  structure,  however.  In  order  to  increase  the  closed-loop 
damping,  controls  developed  using  the  centralized  and  decentralized  viewpoint 
have  been  developed  and  ar^lyzed.  A  centralized  controller,  where  both  actu¬ 
ator  and  structural  velo  Ities  are  fed  back,  can  substantially  increase  the 
system  damping.  The  feedback  gains  must  be  carefully  selected  so  that  the 
closed  loop  is  stable  with  good  response  characteristics.  A  decentralized 
controller  with  the  force  controller  based  on  structural  velocity  feedback  can 
also  substantially  increase  the  system  damping. 

Our  analysis  of  the  closed-loop  characteristics  has  been  based  on  direct 
evaluation  of  the  closed-loop  poles.  Such  a  direct  approach  is  generally  dif¬ 
ficult  to  use  as  a  method  for  design  or  modification  of  a  force  controller. 

In  the  case  of  decentralized  control,  use  of  a  robustness  test  has  been  demon¬ 
strated  as  a  way  of  evaluating  the  effects  of  the  actuator  dynamics  on  the 
closed  loop.  This  robustness  condition  may  prove  effective  as  a  basis  for 
design  or  modification  of  a  force  controller. 
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Although  the  two  examples  do  not  demonstrate  it,  in  more  complicated 
cases,  there  is  the  possibility  for  closed-loop  instability  unless  the  force 
controller  and  actuator  controller  are  carefully  selected.  The  robustness 
demonstrated  in  the  two  examples  is  a  consequence  of  the  use  of  a  force  con¬ 
troller  defined  in  terms  of  constant  gain  structural  velocity  feedback. 

3.6  Conclusions 

3.6.1  Extensions 

A  rather  detailed  theory  has  been  developed  for  large  space  structures 
controlled  by  certain  classes  of  electromechanical  actuators.  This  theory  ha 
been  based  on  rather  c  'ecific  assumptions  about  the  actuator  characteristics 
and  the  assumed  controller  forms.  The  specifics  of  these  assumptions  are  not 
critical  and  extensions  in  several  directions  can  be  indicated. 

Our  development  has  been  based  on  specific  assumptions  about  reaction 
type  actuators  and  transmission  type  actuators,  as  characterized  by  the  schem 
atics  in  Figures  3-1  and  3-2.  Although  it  appears  that  many  physical  actua¬ 
tors  do  fall  within  one  of  these  two  classes,  there  is  no  doubt  that  a  paral¬ 
lel  theory  could  be  developed  for  electromechanical  actuators  where  slightly 
different  assumptions  are  required.  It  should  be  clear  that  the  qualitative 
system  theoretic  characteristics  of  reaction  type  actuators  and  transmission 
type  actuators  are  similar;  their  characteristics  differ  only  in  detail. 

Hence  our  development  can  be  viewed  as  indicating  qualitative  system  charac¬ 
teristics  of  any  electromechanical  actuator. 

Development  of  the  material  in  Section  3.3  on  centralized  control  was 
illustrated  solely  through  the  use  of  a  constant  gain  actuator  and  structural 
velocity  feedback.  Such  a  simple  control  form  is  natural  where  the  objective 
is  augmentation  of  system  damping.  However,  in  many  cases  other  variables, 
e.g.,  displacements  and  accelerations,  may  be  available  for  feedback.  Fur¬ 
ther,  control  laws  based  on  use  of  optimal  linear-quadratic-Gaussian  theory 
may  be  desired,  where  a  filter  or  observer  is  incorporated  as  part  of  the  con 
troller.  Each  of  these  cases  represents  a  viable  procedure  for  developing  a 
centralized  control  law.  The  key  idea  is  that,  no  matter  what  control  design 
procedures  are  used,  the  design  model  is  assumed  to  include  both  structural 
and  actuator  dynamics. 

Extension  of  the  material  in  Section  3.4  on  decentralized  control  is 
possible,  in  the  sense  t  ac  modifications  to  the  developed  forms  for  both 
force  controller  and  actuator  controllers  can  be  made.  The  development  was 
illustrated  by  considering  a  force  controller  defined  in  terms  of  constant 
gain  structural  velocity  feedback.  The  form  of  the  actuator  controllers  sug¬ 
gested  is  consistent  with  that  class  of  force  controllers.  However,  other 
classes  of  force  controllers  could  be  used,  including  other  output  feedback 
forms  or  control  based  on  optimal  linear-quadratic-Gaussian  theory.  In  prin¬ 
ciple,  the  actuator  controllers  specified  in  Section  3.4  could  be  used  in  con 
junction  with  any  force  controller;  the  only  limitation  arises  from  the  fact 
that  the  force  controller  and  actuator  controllers  must  be  realizable.  Thus 
there  are  many  extensions  to  the  developments  in  Section  3.4  that  could  be 


made.  The  key  idea  is  that  the  force  controller  be  chosen  to  suitably  control 
the  structure,  ignoring  actuator  dynamics,  and  the  actuator  controllers  be 
chosen  to  suitably  suppress  the  actuator  dynamics. 

3.6.2  Summary 

There  has  bei’n  substantial  research  into  the  problems  of  active  control 
of  large  space  structures.  Many  issues  have  been  carefully  addressed  and  a 
coherent  theory  is  evolving. 

With  few  exceptions,  the  effects  of  actuator  dynamics  as  part  of  the 
closed- loop  control  scheme  have  been  ignored.  It  is  the  premise  of  this  work 
that  actuator  dynam  3  may  play  an  important  role  in  the  feedback  control  of 
large  space  structures.  Specifically,  if  actuator  dynamics  are  ignored  in  the 
control  design  process,  their  presence  in  an  actual  closed- loop  system  may 
tend  to  be  destabilizing  due  to  neglected  actuator  phase  shifts.  Such  an  un¬ 
desirable  possibility  is  thought  to  be  more  likely  precisely  in  the  class  of 
control  problems  considered,  namely  where  a  lightly  damped  elastic  structure 
is  coupled  to  actuators  with  damped  oscillatory  dynamics.  There  is  also  the 
possibility  that  actuator  dynamics  play  a  desirable  role  in  the  closed  loop. 
Their  filtering  effect  provides  gain  roll-off  at  low  frequencies  and  at  high 
frequencies;  thus,  electromechanical  actuators  tend  to  suppress  coupling  with 
low-frequency  rigid  body  modes  and  with  high-frequency  unmodelled  modes.  In 
any  event,  the  presence  of  actuators  in  the  closed  loop  should  be  carefully 
taken  into  account. 


Our  objective  has  been  to  develop  a  framework  for  the  control  design 
process  where  effects  of  actuator  dynamics  are  not  completely  ignored.  Two 
different  viewpoints  have  been  suggested.  It  is  impossible  to  select  either 
the  centralized  design  viewpoint  or  the  decentralized  design  viewpoint  as  pre¬ 
ferable;  each  approach  should  have  sufficient  flexibility  to  allow  development 
of  suitable  control  strategies  for  a  large  space  structure  with  multiple  elec¬ 
tromechanical  actuators. 


A  detailed  examination  has  been  made  of  the  control  problems,  associated 
with  each  of  these  design  viewpoints.  It  is  hoped  that  this  work  will  serve 
to  f ocu8  additional  attention  on  the  role  of  actuators  and  other  instrumenta¬ 
tion  as  a  critical  part  of  closed-loop  control  of  large  space  structures. 
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SECTION  4 


MODAL-SPRING  PLUS  MODAL -DASHPOT  DESIGN 
OF  OUTPUT  FEEDBACK  VIBRATION  CONTROLLERS 


A  modal-spring  plus  modal-dashpot  design  of  displacement  and  velocity 
output  feedback  vibration  controllers  was  performed  for  the  20-mode  VCOSS 
example  model.  Only  the  10  primary  modes  (those  modes  which  contribute  the 
most  to  the  RMS  line-of-sight  errors)  we.e  considered  in  the  design  process, 
but  all  20  modes  were  used  in  the  evaluation.  The  20-mode  closed-loop  system 
is  as  stable  as  expected  (see  Reference  4-1),  and  each  of  the  20  modes  has  an 
Increase  in  frequency  or  damping  ratio  or  both. 

The  purpose  of  the  design  was  to  increase  the  frequency  of  some  primary 
modes  by  "modal  springs"  and  the  damping  ratio  of  some  other  primary  modes  by 
"modal  dashpots".  In  other  words,  it  was  desired  to  have  some  closed-loop 
poles  with  larger  imaginary  parts,  and  some  others  to  have  more  negative  real 
parts.  As  in  the  preceding  sections,  the  nine  collocated  actuator /sensor  pairs 
were  used.  The  rows  of  the  10  x  9  modal  actuator  influence  matrix  (and  equiva¬ 
lently,  the  columns  of  the  9  x  10  modal  sensor  influence  matrices)  were  very 
dependent  on  each  other.  Since  it  had  been  observed  that  high  feedback  gains 
would  result  in  large  spillover,  one  design  criterion  was  to  restrain  them  from 
extremely  high  gains.  Feedback  gains  for  both  displacement  and  velocity  feed¬ 
back  can  be  greatly  reduced  by  requiring  the  rows  of  the  modal  actuator  influ¬ 
ence  matrix  (and  the  columns  of  the  modal  sensor  influence  matrices)  to  have  a 
high  degree  of  linear  independence.  A  sequential  "Gramian  test"  was  then 
developed  and  implemented  in  PL/I  which  would  choose  the  most  linearly  inde¬ 
pendent  rows  (or  columns).  The  Gramian  test  indicated  that  only  5  to  6  of  the 
10  primary  modes  should  be  used  in  the  design,  since  their  modal  actuator/ 
sensor  influences  were  most  independent. 

The  displacement  and  velocity  feedback  gain  matrices,  GD  and  Gy,  re¬ 
spectively,  were  computed  using  the  pseudoinverses  that  were  initially  consid¬ 
ered  by  Canavin  for  modal  dashpots.  Damping  was  added  to  modes  9,  14,  22, 

23,  33,  and  34,  while  frequency  increases  were  applied  to  modes  7,  9,  12,  22, 
and  23.  The  desired  damping  ratio,  5^,  and  the  desired  frequency,  Wi(j, 
for  these  modes  are  shown  in  Table  4-1.  The  resulting  gain  matrices,  Gp  and 
Gy,  are  shown  in  Tables  4-2  and  4-3,  respectively. 

These  two  feedback  gain  matrices  were  simultaneously  put  into  the  20-mode 
model  and  the  closed-loop  poles  were  evaluated.  Table  4-4  lists  the  real  and 
imaginary  parts  of  the  poles,  their  frequencies,  and  damping  ratios.  Recall 
that  all  20  modes  are  assumed  to  have  an  inherent  damping  ratio  =  0.001. 

The  corresponding  open-loop  modes  are  also  indicated  in  the  table.  Linking  the 
closed-loop  modes  to  the  open-loop  modes  was  accomplished  by  tracing  the  root 
locus  when  the  gain  matrices  were  gradually  scaled  down.  Specifically,  the 
gain  matrices,  Gp  and  Gy,  were  first  replaced  by  yGp  and  yGy,  respec¬ 
tively;  the  closed-loop  poles  were  then  computed  for  discrete  samples  of  y 
between  0  and  1.  When  y  ■  l,  the  poles  correspond  to  the  closed-loop  modes, 
and  when  y  ■>  0,  they  correspond  to  the  open-loop  modes.  Figure  4-1  shows  such 
a  portion  of  the  root  locus. 
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Velocity  output  feedback  gain  ma 
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It  should  be  noted  even  though  many  poles  other  than  those  used  in  design 
also  moved  due  to  spillover,  the  closed-loop  system  is  stable.  The  reduced- 
order  design  of  moda 1-spr ing/moda 1-dashpot  type  of  output  feedback  controllers 
can  guarantee  full-order  closed-loop  stability,  as  was  proved  analytically  and 
demonstrated  numerically  by  a  simple  free-free  beam  in  Reference  4-1.  This  re¬ 
sult  is  again  demonstrated  numerically  by  this  rather  complicated  example. 
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SECTION  5 


A  CONTROLLED  EXPERIMENT  FOR  ACOSS  DESIGN 


5.1  Motivation 

Attempts  to  synthesize  an  active  control  strategy  for  ACOSS  Model  No.  2 
that  is  stable,  that  accommodates  the  effects  of  a  broadband  disturbance,  and 
that  meets  stringent  specifications  on  the  error  in  line-of-sight  (LOS)  rota¬ 
tion  have  so  far  been  unsuccessful.  Early  synthesis  efforts  using  several 
different,  but  complementary,  approaches  to  the  controller  design  appeared 
promising  at  the  level  of  a  10-mode  design  model  [5-1].  Subsequent  evaluation 
of  these  designs  on  a  20-mode  evaluation  model  proved  disappointing,  exhibit¬ 
ing  spillover-induced  instability.  Further  attempts  to  modify  the  designs  so 
as  to  obtain  the  desired  stability,  disturbance  accommodation,  and  LOS  per¬ 
formance  were  similarly  unsuccessful.  Except  in  isolated  instances  in  which 
unrealistic  assumptions  were  invoked,  the  various  designs  either  failed  to 
meet  LOS  performance  specifications  at  the  design  model  level,  or  proved  un¬ 
stable  in  evaluation.  It  is  not  our  intent  to  give  the  details  of  these  at¬ 
tempted  designs  here.  Rather,  the  purpose  of  the  present  section  is  to  out¬ 
line  a  new  approach  to  the  process  of  synthesizing  active  controllers  that 
holds  high  promise  for  success  in  the  sense  described  above. 

At  least  two  reasonable  alternatives  come  to  mind  at  this  point.  One  is 
to  try  to  answer  the  question:  "What  went  wrong?"  That  is,  to  focus  on  dis¬ 
covering  why  the  synthesis  attempts  so  far  have  not  succeeded,  employing 
available  analysis  tools  to  suggest  ways  of  improving  and  reapplying  the  ap¬ 
proaches  already  tried.  Another  is  to  try  to  answer  the  question:  "What  is 
going  on?”  That  is,  to  reexamine  carefully  and  systematically  all  the  princi¬ 
pal  steps  in  the  synthesis  process,  not  just  the  controller  design  itself,  and 
allowing  what  is  observed  to  suggest  new  approaches  to  the  synthesis  process. 

There  is  a  growing  conviction  that  the  latter  alternative  is  the  course 
that  should  be  followed,  principally  because  it  is  the  more  fundamental  one. 
Since  it  subsumes  the  former,  it  is  certain  to  take  longer,  but  can  be  ex¬ 
pected  to  provide  a  much  deeper  understanding  of  the  synthesis  process,  and 
suggest  answers  to  questions  that  perhaps  have  not  arisen  yet.  This 
conviction  has  been  reaffirmed  through  a  recent  editorial  by  a  well-known 
expert  in  the  field  of  control,  who  has  proposed  a  major  and  intentional 
effort  in  "experimental  . oncrol  science"  [5-2] . 

5.2  Preliminary  Observations 


In  order  to  set  the  stage  properly  for  the  approach  to  be  presented,  a 
few  preliminary  observations  are  needed. 

There  are  at  least  four  principal  elements  that  make  up  the  complete 
process  of  synthesizing  an  active  control  strategy  for  a  complicated  structure 
such  as  ACOSS  Model  No.  2.  First  is  the  selection  of  a  basic  structural  de¬ 
sign  to  work  with.  In  the  case  of  Model  No.  2,  the  original  design  has  been 
substantially  revised  several  times.  In  the  present  discussion,  the  reasons 


for,  or  the  details  of,  these  revisions  are  not  as  important  as  the  recogni¬ 
tion  that  most  of  them  have  substantial  influence  upon  the  control  synthesis 
process.  Second  is  the  selection  of  two  or  more  distinct  mathematical  models 
of  the  structure,  one  for  design,  and  one  or  more  for  evaluation  of  the  de¬ 
sign.  Essentially,  this  may  reduce  to  a  (nontrivial)  decision  as  to  which 
modes  of  the  original  (finite  element)  structural  model  to  retain  for  the  var¬ 
ious  purposes.  Third  is  a  selection  of  active  devices  (sensors  and  actuators) 
to  implement  the  control  and  observation  strategies.  Fourth  is  a  selection  of 
an  algorithm  for  determining  the  feedback  structure  of  the  controller. 

It  is  import  .int  to  note  in  passing  that,  whereas  each  of  these  elements 
when  viewed  separately  is  quite  amenable  to  sophisticated  mathematical  analy¬ 
sis,  the  entire  synthesis  process  that  incorporates  the  interactions  between 
these  elements  is  nu  nearly  so  amenable,  if  at  all.  This  fact  represents  a 
substantial  obstacle  to  obtaining  a  satisfactory  synthesis. 

In  the  designs  attempted  to  date,  attention  has  focused  overwhelmingly 
on  the  fourth  element:  algorithms  for  determining  the  controller  feedback 
structure.  Selections  of  the  first  two  elements  (structural  design  and 
reduced-order  model  selection)  have  been  made  with  a  certain  arbitrariness, 
principally  for  simplicity  and  convenience,  and  the  third  element  (active 
devices)  has  been  selected  with  certain  somewhat  arbitrary  restrictions  (e.g. , 
no  nodal  actuators).  These  three  selections  have  remained  fixed  through  all 
of  the  attempts  at  determining  an  appropriate  controller  feedback  structure. 
The  results  suggest  that  in  so  doing,  the  overall  synthesis  problem  has  been 
overconstrained.  Either  a  solution  does  not  exist,  or  the  process  of  finding 
it  has  an  unacceptably  low  probability  of  success.  In  particular,  there 
appears  to  be  a  strong  possibility  that  the  reduced-order  models  (design  and 
evaluation)  and  the  active  device  selections  are  mutually  incompatible. 

We  propose  to  approach  this  impasse  by  conducting  a  "controlled  experi¬ 
ment"  in  which  each  principal  element  in  the  overall  synthesis  process  is  sys¬ 
tematically  examined.  Essential  features  of  the  experiment  are  that: 

(1)  Experiment  "variables”  are  adjusted  only  one  at  a  time. 

(2)  The  effect  of  adjusting  any  experiment  "variable”  is  assessed  by 
conducting  a  complete  end-to-end  synthesis  attempt. 

(3)  Results  of  each  synthesis  attempt  are  analyzed  and,  as  appropriate, 
taken  into  recount  in  subsequent  stages  of  the  experiment. 

5.3  Definition  of  a  "Controlled  Experiment" 

"Variables"  of  the  controlled  experiment  identify  important  features  of 
each  of  the  principal  elements  in  the  complete  synthesis  process.  They  are: 

a:  The  basic  structural  design. 

b  =  (bj,  b2»  bj):  Reduced-order  models. 

bj:  The  criterion  for  ranking  of  structural  modes. 
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b2 :  The  basis  for  retaining  modes  in  the  design  model. 

b3:  The  basis  for  retaining  modes  in  an  evaluation  model. 

c  =  (cj,  C2»  C3):  Active  devices  (sensors  and  actuators). 

cj:  The  types  of  active  devices  to  be  allowed  (e.g.f  axial 
only ) . 

C2:  The  criterion  for  selecting  locations  of  active  devices. 

C3:  The  class  of  reduced-order  models  to  be  considered  in 
sele  "lng  active  devices. 

d:  The  algorithm  for  determining  controller  feedback  structure. 

To  clarify  these  variable  definitions,  we  note  that  the  following  assignment 
of  initial  values  to  these  variables  describes  the  design  process  with  ACOSS 
Model  No.  2  as  it  has  developed  to  date 

a  *■  a(0)  ,  b  <-  b^O)  ,  c  «•  (5-1) 

where 


a  -  Revision  3  of  ACOSS  Model  No.  2,  incorporating  details  of  the 

interaction  between  certain  rigid  components  and  the  associated 
flexible  supports,  and  lightweighting  of  nodal  connecting  ele¬ 
ments 


(0)  A 

m 


the  root-mean-square  (RMS)  LOS-error  induced  by  the  deflection 
of  a  mode  (ignoring  interactions  with  other  inodes)  in  response 
to  the  broadband  disturbance  (Reference  5-1,  Section  5.2.1) 


(0)  A 


the  first  10  modes. 


ranked  according  to  criterion  b 


(0) 

1 


(0)  & 


the  first  20  modes,  ranked  according  to  criterion 


c<°>  A 


only  node-connecting  elements  (e.g. 
motion  sensors) 


rnerabe  r  -d  ampe  r  s , 


axial- 


»  a  best  approximation  (in  the  least-squares  sense)  to  each 

component  of  the  modal  representation  (or  modal  coefficient)  of 
the  LOS-rotation  vector  by  linear  combinations  of  node  shapes 
corresponding  to  a  candidate  selection  of  active  devices 
(Reference  5-1,  Section  3) 
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c<°)  -  the  modal  representation  of  the  LOS-rotatlon  vector  in  crite- 

-*  (0) 

rlon  c,  reflects  only  those  modes  Included  In  design  model 
h(°)  l 
b2 

Variable  “d"  has  assumed  a  number  of  values  in  the  various  design  attempts, 
Including 


,(0)  4 


basic  disturbance  rejection  control  (DRC)  design  in  which 
complete  knowledge  of  the  disturbance  statistics  is  assumed 
(Reference  5-1,  Section  4) 


,(D  A 


lines.  quadratic  optimal  control  design,  in  which  no  explicit 
knowledge  of  the  disturbance  statistics  is  assumed  (Reference 
5-1,  Section  5) 


,(2)  4 


algorithm  dv  augmented  by  the  incorporation  of  measures  to 
alleviate  spillover  (Reference  5-1,  Section  6) 


among  others. 

In  the  initial  phase  of  the  experiment,  we  propose  to  fix  the  values  of 
variables  "a”  and  "d",  and  focus  on  the  effects  of  changes  in  variables  “b" 
and  "c” — particularly  the  interaction  between  them.  The  fixed  values  chosen 
for  variables  "a”  and  “d"  are 


a  «-  a 


d*d' 


where 

a(00)  A  ggyigion  i  0f  ACOSS  Model  No.  2,  a  stiffer  structural  design 
which  does  not  contain  the  lightweighting  feature  of  design 
a(0) 

A  pictorial  representation  of  the  overall  synthesis  process  is  shown  in  Figure 
5-1.  A  most  significant  feature  of  this  process  is  the  interaction  of  vari¬ 
ables  "b"  and  "c",  which  has  not  been  given  sufficient  attention  in  previous 
design  attempts. 

Specific  classes  of  values  for  variables  "b"  and  "c"  that  appear  to  be 
of  interest  Include  those  listed  below.  These  listings  are  most  conveniently 
expressed  as  functions  of  other  experiment  variables  and  of  additional  vari¬ 
ables  nc,  n^  that  may  be  considered  parameters  of  the  experiment. 

b^  «  criterion  bj^  augmented  by  specific  a  priori  informa¬ 
tion  on  physical  or  modal  deformations 
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STRUCTURAL 
DESIGN  SELECTION 
(a) 


Figure  5-1.  Process  for  active  control  synthesis 


b21)(nC»bl) 


the  first  n^  modes,  ranked  according  to  criterion  b^ 
(i.e.,  the  value  being  defined  for  b£  is  a  function 
of  whatever  value  has  been  selected  for  variable  bj ) 


b22)(vV 


the  modes  selected  by  criterion  b 


(1) 


(n^  ,b^ )  augmented 


by  all  modes  whose  natural  frequencies  lie  in  the  range 
associated  with  those  selections  (i.e.,  not  allowing  an 
interlacing  of  controlled  and  residual  nodes) 


b3U<W 


b<2)(VV 


the  first  n  modes  (n  >  n  ),  ranked  according  to  cri¬ 
terion  bL 

the  modes  selected  by  criterion  b^^i^.b^)  augmented 
by  all  modes  whose  natural  frequencies  lie  in  the  range 
associated  with  those  selections 


(1)  A 


.(2)  A 


devices  c^  together  with  devices  located  at  nodes  and 
capable  of  translational  action  only 

devices  c  j  ^  augmented  with  a  rotational  action  capabil¬ 
ity  at  nodes 


(1)  A 


.(0) 

'2 


criterion  c^~'  augmented  by  specific  a  priori  informa¬ 


tion  on  physical  or  modal  deformations 


C32)(c2’b3} 


c3  (c2*b2^  *  tbe  mo<*a^  representation  of  the  LOS  rotation  vector  in 

criterion  c-2  reflects  all  modes  included  in  design 
model  b2 

the  modal  representation  of  the  LOS  rotation  vector  in 
criterion  C2  reflects  all  modes  included  in  evaluation 
model  b3 

A  pictorial  representation  of  typical  stages  of  the  experiment  in  terms 
of  experiment  variables  is  shown  in  Table  5-1.  Each  column  represents  an 
end-to-end  attempt  of  the  entire  synthesis  process.  The  succession  from  one 
column  to  the  next  (from  left  to  right)  involves  changing  only  one  variable. 

5.4  Expected  Outcome 

It  is  expected  that  the  conduct  of  this  experiment  will  lead  to  a  deeper 
understanding  of  the  complete  process  of  active  control  synthesis.  It  is 
especially  important  to  understand  more  clearly  the  essential  elements  of 
selecting  appropriate  reduced-order  models,  selecting  active  devices  for 
actuation  and  sensing,  and  the  nature  of  the  relationship  between  them.  In 


particular,  it  is  expected  that  the  results  of  the  experiment  will  generate  a 
process  by  which  a  reduced-order  design  model  and  a  set  of  active  devices 
which  are  mutually  compatible  can  be  found.  Such  mutual  compatibility  is  the 
key  to  successful  synthesis  for  active  control. 


SECTION  5 

LIST  OF  REFERENCES 


R.  Strunce,  et  al. ,  ACOSS  Eleven  Second  Semiannual  Technical  Report, 
Volume  2:  Active  Controller  Designs,  Report  CSDL-R-1583,  Charles  Stark 
Draper  Laboratory,  Cambridge,  MA,  August  1982. 

Ho,  Y.-C.,  "Is  It  Application  or  Is  It  Experimental  Science?”  (Edi¬ 
torial),  IEEE  Trans.  Automatic  Control,  Vol.  AC-27,  No.  6,  December 
1982,  p.  1142. 


P ADC /DAP 

OR  IFF  I SS  ARB  NY  13441 
-MOC/TSLD 

I.RIFFISS  AFB  NY  13441 


ADMINISTRATOR 
DDF  TECH  INF  CTR 
YrTN:  DTIC-DDA 

AMERON  STA  BG  5 
ALEXANDRIA  VA  22314 

'^arles  Stark  Draper  Lab 
rD5  Technology  Sq>  ^re 
Cambridge.  MA  02  J 39 


■  i.arles  Stark  Draper  Lab 
r+vtn:  Dr.  Keto  Soosaar 
5^/5  Technology  Square 
M  S.  -95 

Cambridge.  MA  02139 

NASA  Headquarters 
ATTN:  Mr.  J.  B.  Dahlgren 
Code  RTH-6 

Washington,  DC  20546 


Charles  Stark  Draper  Lab 
Attn:  Mr.  R.  Strunce 
-55  Technology  Square 
M  S.  -60 

Cambridge,  MA  02139 

Charles  Stark  Draper  Lab 
Attn:  Dr.  Daniel  R.  Hegg 
C 55  Technology  Square 
M.  3.  -60 

"ambridge.  MA  02139 
ARPA/STO 

Attn  Lt  Col  A.  Herzberg 
-  **00  Wilson  Blvd 
Arlington,  VA  22209 


ARPA/STO 

Attn:  Mi  j  E.  Dietz 
\ 400  Wilson  Blvd 
Arlington,  VA  22209 


DL- 


Riverside  Research  Institute 
Attn  Mr.  A.  DeVilliers 
1701  N.  Ft.  Myer  Drive.  Suite  711 
Arlington.  VA  22209 


Riverside  Research 

Attn  HALO  Library.  Mr.  Bob  Passut 
1701  N.  Ft.  Myer  Drive 
Arlington.  VA  22209 


'tek  Corp 

Optical  Systems  Division 
10  Maguire  Rd. 

,-exington.  MA  02173 


Perkin  Elmer  Corp 
Attn:  Mr.  H.  Levenstein 
Electro  Optical  Division 
nain  Avenue 
Norwalk.  CT  06856 

Hughes  Aircraft  Company 
Attn:  Mr.  George  Speak 
MS.  B_1 56 

Culver  City.  CA  09230 


Hughes  Aircraft  Company 
Attn:  Mr.  Ken  Beale 
Centinela  Teale  Sts 
Culver  City.  CA  90230 


Air  Force  Flight  Dynamics  Lab 

Attn+  Dr.  Lynn  Rogers 

Wright  Patterson  AFB.  0*  45433 


i-’FWL/FIBG 

Attn:  Mr.  Jerome  Pearson 
aright  Patterson  AFB.  OH  45433 


->ir  Force  Wright  Aero  Lab.  FIGC 

Actn  Siva  S  Banda 

Wright  Patterson  AFB.  OH  45433 


Air  Force  Institute  of  Technology 
Attn.  Prt>f.  R.  Calico/ENY 
Wright  Patterson  AFB,  OH  45433 


Aerospace  Corp. 

Attn:  Dr.  G.  T.  Tseng 
2350  E.  El  Segundo  Blvd 
El  Segundo*  CA  90245 


Aerospace  Corp. 

Attn:  Mr.  J.  Mosich 
2350  E.  El  Segundo  Blvd 
El  Segundo*  CA  90245 


Aerospace  Corp/Bldg  125/1054 
Attn:  Mr.  Steve  Burrin 
Advanced  Systems  Tech  Div 
2400  E  El  Segundo  Blvd 
£1  Segundo,  CA  90245 

SD/SD/YLVS 

Attn:  Mr.  Lawrence  Weeks 
P  0.  Box  92960 
Worldu/ay  Postal  Center 
Los  Angeles  CA  90009 

SD/YCD 

Attn.  YCPT/Capt  Gajewski 
P  0.  Box  92960 
Worlduiay  Postal  Center 
Los  Angelel*  CA  90009 

Grumman  Aerospace  Corp 
Attn:  Dr.  A.  Mendelson 
South  Oyster  Bay  Road 
Bethpage.  NY  1171*t 


OUSDR&E/DS 

Attn  Mr.  A.  Bertapelli 
Hoorn  3D 136 

Pentagon,  Washington*  DC  20301 


Jet  Propulsion  Laboratory 
Dr  S.  Szermay 
4800  Oak  Grove  Drive 
Pasadena*  CA  91103 


DL-3 


MIT/Lincoln  Laboratory 
Attn:  S.  Wright 
P  0.  Box  73 
Lexington.  MA  02173 


MIT/Lincoln  Laboratory 
Attn:  Dr.  D.  Hyland 
P.  □.  Box  73 
Lexington.  MA  02173 


MlT/Lincoln  Laboratory 
Attn:  Dr.  N.  Smith 
P  0.  Box  73 
Lexington.  MA  02173 


Control  Dynamics  Co. 

Attn:  Dr.  Sherman  Seltzer 
Suite  1414  Executive  Plaza 
555  Sparkman  Drive 
Huntsville.  AL  35805 

Lockheed  Space  Missile  Corp. 

Attn:  A  A.  Woods.  Jr.  ,  0/62-E6 
p  0.  Box  504 

Sunnyvale,  California  94088-3504 


Lockheed  Missiles  Space  Co. 
Attn:  Mr.  Paul  Williamson 
1251  Hanover  St. 

Palo  Alto,  CA  94304 


General  Dynamics 
Attn:  Ray  Halstenberg 
Convair  Division 
5001  Keary  Villa  Rd 
San  Diego.  CA  92123 

ST  I 

Aetn:  Mr.  R.  C.  Stroud 
20065  Stevens  Creek  Blvd. 
Cupertiono,  CA  95014 


NASA  Langley  Research  Ctr 
Attn:  Dr.  Earle  K.  Huckins  III 
Dr.  M.  F.  Card 

Langley  Station,  Bldg  1293B.  MS  230 
Hampton,  VA  23665 


NASA  Johnson  Space  Center 
Attn.  Robert  Piland 
Ms.  EA 

Houston/  TX  77058 


McDonald  Douglas  Corp 

Attn:  Mr.  Read  Johnson 

Douglas  Missile  Space  Systems  Div 

5301  Bulsa  Ave 

Huntington  Beach/  CA  92607 

Integrated  Systems  Inc. 

Attn:  Dr.  N  K.  Gupta  and  M.  G.  Lyons 
151  University  Avenue/  Suite  400 
Palo  AltOz  California  94301 


TRW  Defense  Space  Sys  Group  Inc. 

Attn:  Ralph  Iwens 

Bldg  82/2054 

One  Space  Park 

Redondo  Beach/  CA  90278 

TRW 

Attn  Mr.  Len  Pincus 
uldg  R-5/  Room  2031 
Redondo  Beach/  CA  90278 


Department  of  the  NAVY 
Attn  Dr  K.  T.  Alfriend 
Naval  Research  Laboratory 
Code  7920 

Washington,  DC  20375 

Airesearch  Manuf.  Co.  of  Calif. 
Attn:  Mr.  Oscar  Buchmann 
2525  West  190th  St. 

Torrance/  CA  90509 


Analytic  Decisions/  In.. 
Attn  Mr.  Gary  Glaser 
i401  Wilson  Blv. 
Arlington/  VA  22209 


Ford  Aerospace  &  Communications  Corp. 

Drs  I  P.  Leliakov  and  P.  Barba/  MS/GBO 

3939  Fabian  way 

Palo  Altoz  California  94304 


Center  for  Analysis 
Mr  James  Justice 
13  Corporate  Plaza 
Newport  Beach,  CA  92660 


W  J  Schafer  Associates 
Dr  R  Kappesser 
Suite  800 

1901  Fort  Meyer  Drive 
Arlington.  VA  22209 

General  Research  Corp 
Attn  Mr  Thomas  Zakriewski 
"Cj55  Old  Springhouse  Road 
McLean,  VA  22101 


Air  Force  Ueapons  Laboratory 
Attn:  Lt  Col  D.  Washburn 
ARAA 

nirt land  AFB,  NM  87117 


Harman  Sciences  Corp. 

Attn:  Dr.  Walter  E.  Ware 
1500  Garden  of  the  Gods  Road 
P  0  Box  7463 

Colorado  Springs,  CO  80933 
MR  J,  Inc. 

10400  Eaton  Place 
Suite  300 
Fairfax,  VA  22030 


Photon  Research  Associates 
Attn.  mr.  Jim  Myer 
P  0.  Box  1318 
La  Jolla,  CA  92039 


Rockwell  International 

Attn:  Russell  Loftman  (Space  Systems  Group) 
(Mail  Code  -  SL56) 

12214  Lakewood  Blvd. 

Downey,  CA  90241 

Science  Applications,  Inc. 

Attn:  Mr.  Richard  Ryan 
3  Preston  Court 
Bedford,  MA  01730 


U.  S.  Army  Missile  Command 
Attn:  DRSMI-RAS/Mr .  Fred  Haak 

Redstone  Arsenal,  AL 


1 


56 


Naval  Electronic  Systems  Command 
Attn:  Mr,  Charles  Good 
PME_106-4 
National  Center  I 
Washington*  DC  20360 

Lockheed  Palo  Alto  Research  Laboratory 

Attn:  Dr.  J.  N.  Aubrun.  0/52-56 

3251  Hanover  Street 

Palo  Alto,  California  94304-1187 


U  S.  Army/DARCOM 

Attn:  Mr  Bernie  Chasnov 

AMC  Bldg 

5001  Eisenhower  Ave 
Alexandria,  VA  22333 

R  Carman 
RADC/OCSE 


Honeywell  Inc. 

Attn:  Dr.  Thomas  B.  Cunningham 
Attn  Dr.  Michael  F.  Barrett 
2600  Ridgway  Parkway  MN  17-2375 
Minneapolis,  MN  55413 

NASA  Marshal  Space  Flight  Center 
Attn.  Dr.  J.  C.  Blair,  ED01 
Henry  B.  Waites 

Marshal  Space  Flight  Center,  AL  35812 

q 

TRW 

Attn:  Robert  Benhabib 
Bidg  82/2024 
One  Space  Park 
Redondo  Beach,  CA  90278 

NASA  Langley  Research  Center 
Attn:  Dr.  L.  Pinson 
MS  -  230 

Hampton,  VA  23665 


H.  R.  Textron 

Attn.  Mr.  Richard  Quartararo 
2485  McCabe  Way 
Irvine,  CA  92714 


Naval  Research  Lab 
Attn.  W.  Bennett 
Mail  Code:  7926 
Washington,  DC  20375 


DL-7 


i  « 

\  MISSION  5 

\  of  \ 

^  Rome  Air  Development  Center  ? 

jj  RA.VC  pta.nl>  and  executu  research,  development,  tut  and  \ 

*  selected  acquisition  programs  in  s upport  o{,  Command,  Control  Q. 
5  Comunccationi  and  Intelligence.  {Ch}  activities .  Technical  f 

*  and  engineering  Support  within  areas  oh  technical  competence  C 

S  AS  Provided  to  ESP  Program  Ohhicu  IPOs]  and  other.  BSD  § 

^  elements.  The  principal  technical  mission  areas  are  \ 

k  communications ,  electromag netic  guidance  and  control,  sur-  a 
\  veillance  oh  ground  and  aerospace  objects,  intelligence  data  \ 

*  collection  and  handling,  inhormation  system  technology,  V 

r  ionospheric  propagation,  solid  state  sciencu ,  microwave 
{  physics  and  electronic  reliability,  maintainability  and 
>  compatibility . 


